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Preface 


I  began  this  research  as  a  repeat  participant  in  the  “barf  chair-’  last  fall  be¬ 
fore  funds  had  been  completely  exhausted  for  data  collection.  It  certainly  wasn't 
glamorous,  but  I  had  become  immediately  interested  upon  my  entering  AFIT  in 
the  leading-edge  research  of  the  motion  sickness  research  group.  I  was  saddened  L* 
Dr.  Bill  Chelen’s  departure  as  a  result  of  the  budget  reductions,  but  was  glad  to 
have  *  thesis  topic  in  this  fascinating  area  and  the  support  of  the  group,  who  took 
me  onboard  despite  the  throngs  of  other  students  at  their  office  doors. 

I  recently  read  that  one  person  can  learn  when  left  on  his  own.  But  that  person 
must  be  in  the  company  of  others  in  order  to  learn  well.  This  statement  summarizes 
the  thesis  experience  for  me.  I  am  glad  this  page  provides  a  forum  to  thank  those 
who  accompanied  me. 

My  deepest  thanks  go  to  my  wife,  Kathy,  and  two  children,  Andrew  and 
Patrick.  I  have  appreciated  their  endless  understanding  and  support,  but  especially 
the  way  that  their  presence  and  love  were  frequent  reminders  of  the  immeasur¬ 
able  value  of  my  family.  Maj  Steve  Rogers,  LtCol  T.S.  Kelso,  and  Dr.  Matthew 
Kr’  isky,  my  committee  members,  provided  the  right  mixture  of  guidance  and  ex¬ 
pectation.  They  were  mentors,  not  taskmasters,  for  which  I  am  most  grateful.  Their 
enthusiasm  was  catching  and  sustaining.  Capt  Chuck  Wright  started  me  on  the  way 
by  graciously  offering  his  program  and  expertise  in  answer  to  my  initial  question  of 
"How?’’  Capts  Dave  Swanson  and  Art  Earl  were  dependable  allies  and  continuous 
sources  of  workarounds  and  shortcuts  which  kept  bearable  the  tedium  of  program¬ 
ming.  Capt  Ken  Langford  and  Mr.  Ron  Milliron  were  mv  “friends  in  high  places" 
that  kept  the  computer  workstations  willing  to  work  with  me. 

I  thoroughly  enjoyed  this  epic  journey  into  the  world  of  motion  sickness,  EEGs, 
.J  computer  programming.  Furthermore,  I  will  never  experience  motion  sickness 
in  the  same  way  again. 
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I  hope  that  in  this  search  for  pattern,  Grey  Walter  would  have  seen  the  essence 
of  science.  But  moreover,  I  hope  that  in  the  pursuit  of  science,  more  people  will 
recognize  the  work  of  a  loving  and  living  God. 


Dwight  Andrew  Roblyer 
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Abstract 

A  toposcope  was  constructed  as  a  new  tool  to  study  signal  frequency  and  phase 
relationships  in  electroencephalogram  (EEG)  records  collected  while  subjects  were 
experiencing  motion  sickness.  This  new  tool,  named  TOPOS,  is  a  software- based, 
multi-featured  version  of  Grey  Walter  and  Harold  Shipton’s  device  which  they  first 
produced  in  the  late  1940s.  The  TOPOS  graphical  display  permits  the  study  of 
instantaneous  frequency  relationships  between  the  input  channels  and  a  reference 
signal  of  fixed  or  varying  frequency.  TOPOS  adds  a  correlation  grid  to  aid  observers 
in  detecting  channel-to-reference  correlation  levels.  Users  can  also  vary  several  dis¬ 
play  parameters  via  menus  to  optimize  the  analysis  environment.  Sinusoidal  test 
inputs  of  known  frequency  produced  recognizable  and  predictable  patterns  on  the 
TOPOS  display,  depending  on  the  existing  channel-to-reference  frequency  relation¬ 
ships.  Motion-sickness-afFected  EEG  was  input  to  TOPOS  in  order  to  study  the 
correlation  between  the  displays  of  each  channel  and  four  separate  references:  a  1.5- 
Hz  sinusoid,  and  three  channels  of  the  EEG  itself.  Rapidly  varying  correlations  were 
observed  in  each  case. 
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A  TOPOSCOPIC  INVESTIGATION 
OF  BRAIN  ELECTRICAL  ACTIVITY 
INDUCED  BY  MOTION  SICKNESS 


I.  Introduction 

1.1  Background:  Motion  Sickness  on  Earth  and  in  Space 

Motion  sickness  is  a  common  companion  of  people  traveling  by  land,  sea,  or 
air  (33:29-30).  These  sufferers  find  that  the  slowly  rolling  road,  choppy  seas,  or  air 
turbulence  can  quickly  biing  on  headaches,  sweating,  stomach  discomfort,  nausea, 
and  vomiting — ail  symptoms  of  motion  sickness  ( 4 2 : A 1 ) . 1  Victims  often  either  take 
medication  which  makes  them  drowsy  and  the  symptoms  less  pronounced,  or  stal- 
wartly  suffer  until  the  motion  stops,  they  adapt,  or  emesis  brings  some  relief.  No 
cure  exists,  and  the  exact  etiology  (cause)  of  the  symptoms  is  poorly  understood 
(28:218).  While  wretchedly  uncomfortable  and  annoying,  terrestrial  motion  sick¬ 
ness  is  relatively  harmless,  causing  delays  in  travel  or  task  schedules  but  posing  few 
threats  to  human  life  (33:73). 2 

But  what  can  be  done  when  a  person  experiences  motion  sickness  symptoms 
in  space,  where  task  delays  carry  very  expensive  price  tags  and  compromise  of  a 
crewmember’s  attentiveness  and  fitness  could  risk  his  health,  or  even  his  life?  Even 
after  30  years  of  research,  this  question  remains  without  a  wholly  satisfactory  an¬ 
swer.  Moreover,  it  is  an  important  question,  because  upon  entry  into  orbit  the  first 

'Motion  sickness  w  really  not  a  “sickneus”  at  nil,  but  was  so  labeled  because  of  the  symptoms 
it  produces  (33:28).  Dhenin  calls  it  a  “normal  response  to  an  abnormal  environment"  (9:472). 

2One  significant  exception  is  the  case  of  military  forces.  Motion  sickness  is  a  significant  problem 
for  soldiers  being  transported  into  battle,  explaining  why  •  *c  military  has  sponsoi^u  much  of  the 
search  for  the  causes  and  cures  of  motion  sickness  ( 42: A 1 ) ( 33:31  -  33). 
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physiologic  change  experienced  by  up  to  67  percent  of  US  astronauts  is  space  motion 
sickness  (SMS),  also  known  as  space  adaptation  syndrome  (SAS)  (15:1). 

Space  adaptation  syndrome  has  been  called  by  Vanderploeg,  et  al.  the  most 
clinically  significant  medical  phenomenon  during  the  first  several  days  of  spaceflight 
(8:1185).  During  that  period,  space  travelers  experience  some  symptoms  common  to 
terrestrial  motion  sickness  and  some  unique  to  SMS-  An  affected  person  may  develop 
a  hypersensitivity  to  motion  or  experience  illusions  of  being  constantly  upside  down 
or  pitched  forward  (42:A2).  The  astronaut  will  almost  always  have  no  appetite  and 
may  have  jepeated  episodes  of  brief,  but  violent  ,  emesis  (42:  A2).  Nearly  one-third  of 
US  astronauts  have  reported  at  least  one  episode  of  vomiting  while  in  orbit  (25:101) 
While  the  sufferer  rarely  sweats  or  shows  pallor,  he  will  almost  certainly  experience 
a  headache,  malaise,  iack  of  initiative,  and  irritability  (42:A3).  These  symptoms  not 
only  vary  between  individuals,  but  between  separate  trips  into  orbit  by  the  same 
persou  (25:103).  However  it  presents  itself,  the  SMS  syndrome  lasts  for  an  average 
of  3  to  5  days,  until  the  victim  adapts  to  the  weightless  environment  (25:20). 

in  1982,  the  Air  Force  Institute  of  Technology  began  research  into  causes  and 
treatments  of  SMS.  Recently,  this  has  focused  on  the  brain  electrical  activity  which 
appears  as  motion  sickness  symptoms  develop.  The  AFIT  research  team  found  this 
activity  formed  a  pattern  which  was  repeated  among  different  subjects,  and  has 
followed  this  lead  in  search  of  insights  into  the  causes  and  effective  treatment  of 
motion  sickness.  Banducci  and  Vogeu  began  the  process  of  analyzing  these  patterns 
by  separating  them  into  their  component  frequencies  using  spectral  analysis  and 
then  graphically  mapping  the  results  (2;ix)(44:vii).  Their  methods  averaged  data 
across  time  and  based  results  on  interpolated  data.  Because  of  this  methodology, 
there  still  remains  questions  concerning  the  origin  and  instantaneous  propagation  of 
this  abno.iiial  activity  in  the  brain.  These  issues  can  be  further  investigated  using  a 
different  tool — the  toposcope. 
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1.2  Problem  and  Scope 

The  problem  is  to  construct  a  toposcope  as  a  new  tool,  test  its  use  and  flex¬ 
ibility,  and  then  interface  it  with  motion-sickness-affected  EEG.  The  device  will  be 
computer-based  and  coded  in  the  C  programming  language  for  implementation  on 
a  Silicon  Graphic  workstation.  This  brain  mapping  approach  should  provide  a  con¬ 
tinuous,  graphical  display  which  can  focus  on  brain  activity  in  very-low-frequency 
regions  without  the  loss  of  time  resolution  resulting  from  spectral  analysis.  It  will 
be  able  to  derive  its  frequency  selection  from  a  fixed  frequency  source  or  the  signal 
in  a  single  EEG  channel.  The  tool's  performance  will  be  characterized  by  inputting 
known-frequency  signals.  The  toposcope  will  then  be  interfaced  to  a  motion-sickness- 
affected  EEG  record. 

1.3  Role  of  Pattern  Recognition 

The  pursuit  of  information  about,  and  classification  of,  spatial  and  temporal 
patterns  is  the  focus  of  the  field  of  pattern  recognition  (43:5).  In  this  thesis,  the 
principles  of  this  field  will  be  applied  to  a  motion-sickness-affected  EEG  record  to 
search  for  features  or  invariant  attributes  of  SMS  brain  activity  patterns.  This  re¬ 
search  effort  will  involve  a  partnership  between  human  and  machine  in  performing 
this  task.  The  machine — the  toposcope — will  possess  no  automatic  recognition  ca¬ 
pability,  as  is  commonly  associated  with  machines  in  the  field  of  pattern  recognition. 
Instead,  the  toposcope  will  exhibit  the  features  of  the  pattern  which  are  thought  to 
be  the  best  for  recognition:  highlighting  the  frequency  and  spatial  relationships  of 
the  EEG  signals.  Furthermore,  the  toposcope  will  further  aid  recognition  by  tuning 
its  entire  display  to  the  predominant  frequency  of  a  selected  channel.  Therefore, 
by  extracting  significant  features  of  the  EEG  data  from  a  background  of  irrelevant 
detail,  a  human  can  approach  a  much  simpler  task  of  recognizing  any  pattern  present 
and  further  defining  the  paitem  class  encompassing  motion  sickness.  Grey  Walter, 
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one  of  the  toposcope’s  creators  and  a  pioneer  in  brain  activity  research,  wrote  that 
in  the  seeking  of  pattern  lies  the  essence  of  science  (47:69). 

1,4  Assumptions 

The  are  several  key  assumptions  for  this  research  regarding  data  collection 
and  analysis.  Because  of  the  constraint  on  data  availability  as  explained  in  the  first 
assumption,  several  of  those  that  follow  are  similar  to  the  assumptions  proposed  by 
Vogen. 

1.  The  only  EEG  records  available  for  analysis  are  those  collected  during  1991. 
Justification:  Air  Force  sources  could  no  longer  fund  AFIT’s  motion  sickness 
research  after  September  1991  due  to  budgetary  cutbacks.  Limiting  Effects: 
No  additional  data  can  be  collected  for  the  purpose  of  pattern  validation  or  in 
the  case  of  data  loss  or  corruption. 

2.  All  the  equipment  used  in  the  collection  of  EEG  data  was  properly  calibrated 
and  functioning  as  designed.  Justification:  Data  were  collected  prior  to  this 
research  effort  and  Vogen  documented  his  calibration  and  functional  checks 
of  equipment  (44:31).  Limiting  Effects:  If  calibration  or  other  problems  are 
found  in  the  data,  efforts  are  restricted  to  finding  workarounds  or  eliminating 
the  affected  data  set(s). 

3.  The  EEG  data  subject  to  analysis  were  relatively  artifact  free  and  any  artifact 
which  was  undetected  was  insignificant  to  the  results.  Justification:  Data  were 
collected  prior  to  this  research  effort  and  artifacts  were  assessed  as  minimal 
(44:4).  Limiting  Effects:  If  excessive  artifacts  are  discovered  in  the  data,  efforts 
are  restricted  to  finding  workarounds  or  eliminating  the  affected  data  set(s). 

4.  The  research  problem  is  one  of  pattern  recognition.  A  pattern  is  postulated  to 
exist,  but  finding  and  displaying  the  best  feature  of  that  pattern  will  reveal  the 
most  information  about  its  behavior.  Justification:  1)  Brain  activity  patterns 
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are  the  basis  of  the  entire  field  of  brain  topographic  mapping.  2)  Vogen  reported 
a  standard  pattern  during  motion  sickness.  Limiting  Effects:  Focuses  the 
research  on  the  qualitative  detection  of  pattern. 

5.  The  physiologic  changes,  including  EEG,  which  occur  as  motion  sickness  devel¬ 
ops,  correlate  to  the  subjective  level  of  sickness  reported  by  the  subject.  Jus¬ 
tification:  In  1987,  Drylie,  Fix,  and  Gaudreault  developed  an  equation  which 
correlates  these  two  groups  of  data  and  trained  a  neural  network  to  predict 
subjective  levels  of  sickness  using  physiological  data,  according  to  Banducci 
(2:49).  Limiting  Effects:  None.  Rather,  it  permits  the  accurate  labeling  of 
EEG  states  with  the  current  level  of  sickness. 

6.  The  results  will  be  representative  of  those  which  would  be  obtained  from  the 
healthy,  normal  population  experiencing  motion  sickness  symptoms.  Justifica¬ 
tion:  Samples  were  drawn  from  10  healthy,  normal  volunteer  Air  Force  officers 
ranging  from  26  to  41  years  of  age  (18:1153).  Limiting  Effects:  None.  Rather, 
it  forms  a  portion  of  the  basis  for  further  extrapolating  the  results  to  cases  of 
space  motion  sickness  among  astronauts. 

1.5  Format  of  Thesis 

The  next  two  chapters  provide  background  information  for  this  research.  Chap¬ 
ter  II  frames  the  space  motion  sickness  problem  and  emphasizes  the  role  of  the  EEG 
in  potential  solutions.  Chapter  III  then  overviews  three  brain  mapping  techniques 
used  in  deciphering  and  understanding  the  complexities  of  the  EEG.  Chapter  IV 
overviews  the  approach  used  in  each  stage  of  the  thesis  and  explains  the  relation¬ 
ships  between  the  AFIT  display  and  that  of  Walter  and  Shipton  (45:283).  Chapter  V 
reports  on  the  design  and  coding  of  the  AFIT  toposcope.  The  results  from  analyz¬ 
ing  test  signals  and  motion-sickness-affected  EEG  records  using  this  tool  appear  in 
Chapter  VI,  followed  by  conclusions  in  Chapter  VII. 


II.  SMS:  A  Summary  of  Current  Knowledge 


2. 1  Why  SMS  is  Significant 

The  impact  of  SMS  on  space  operations  is  disputed  by  some.  Although  ac¬ 
knowledging  the  syndrome  s  physiological  effects,  Thornton,  et  al.  state, 


*Vhile  it  is  obvious  that  a  person  is  hors  de  combat1  during  vomiting,  this 
is  brief.  Conversely,  trained  astronauts  have  in  every  case,  performed 
assigned  tasks.  There  have  been  two  precautionary  delays  of  scheduled 
EVAs2.  It  isn't  necessarily  easy  . . .  there  is  a  lack  of  initiative,  but  tasks 
trained  for  and  scheduled  are  done  and  done  well  (42:A4). 

However,  Thornton,  et  al.  also  report  that  individuals  with  less  extensive 
training  prior  to  the  flight  have  sometimes  been  unable  to  complete  all  assigned 
tasks  (42:A4).  An  example  was  the  first  journalist  in  space,  a  Tokyo  Broadcasting 
correspondent  who  traveled  with  a  Soviet  resupply  mission  to  the  Mir  space  station 
in  December  1990.  He  had  participated  in  18  months  of  rigorous  training  with  the 
cosmonauts,  but  still  his  reports  back  to  Earth  were  filled  with  comments  about  his 
high  level  of  discomfort  and  lack  of  normal  function  due  to  space  motion  sickness 
and  other  factors  (36:1).  In  addition,  Hettinger  maintains  that  any  period  of  time  in 
which  crew  members  are  disabled  by  SMS  symptoms,  no  matter  how  severe  or  pro¬ 
tracted,  is  “extremely  critical  in  terms  of  the  risks  and  expense  involved”  (14:424). 
A  1988  study  by  Ratino,  et  al.  found  a  correlation  between  increase  in  reaction 
time  and  magnitude  of  motion  sickness  symptoms  (32:223).  Furthermore,  activities 
outside  the  Shuttle  are  not  scheduled  until  after  the  third  day  of  a  flight  and  are  lim¬ 
ited  to  immune  or  adapted  individuals  because  of  the  dangers  of  asp. rating  vomitus 
inside  the  closed  environment  of  a  space  suit  (8:1 185)(30:2S).  This  is  particularly 

'In  a  disabled  condition. 

2 Extravehicular  activities  (space  walks). 
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important  because  emesis  in  SMS  usually  comes  suddenly,  without  the  prior  nausea 
or  other  gastrointestinal  symptoms  that  would  normally  precede  it  on  Earth  (42:A2). 
So.  SMS  poses  not  only  a  nuisance  but  also  a  threat  to  astronauts  and  their  mission. 

2.2  Why  Docs  SMS  Occur ? 

Many  researchers  consider  SMS  just  another  form  of  the  terrestrial  motion 
sickness.  They  also  believe  the  two  syndromes  have  the  roots  of  their  existence  in 
the  neural  mismatch3  hypothesis  first  formed  by  Reason  (33:166).  This  hypothesis 
explains  motion  sickness  symptoms  as  the  body’s  natural  response  when  the  infor¬ 
mation  transmitted  by  the  eyes,  the  vestibular  system,  and  other  receptors  i«  not  the 
information  the  person  expects  to  receive,  based  on  past  experience  (9:474).  This 
memory  of  sensory  inputs  is  accessed  by  a  comparator  center  within  the  central 
nervous  system  (CNS),  which  correlates  them  with  the  current  experience  and  then 
generates  a  response  dependent  upon  the  duration  and  the  intensity  of  the  exist¬ 
ing  mismatch.  Frank  motion  sickness  would  be  the  result  of  an  ‘‘overload"  in  this 
comparator  center. 

For  astronauts,  discordant  sensory'  inputs  are  quick  to  arise  upon  entry  into  or¬ 
bit.  Their  otoliths,  the  vestibular  receptors  dependent  upon  gravity  for  determining 
the  head’s  position,  are  left  without  any  reference  at  all.  Visual  and  tactile  orien¬ 
tation  cues  also  become  ambiguous:  “up"  and  “down”  are  now  redefined  constantly 
by  the  changing  visual  and  tactile  experience  of  a  floating  astronaut  (28:218). 

Although  it  is  widely  accepted  that  SMS  is  a  form  of  motion  sickness,  there  are 
other  theories  which  also  implicate  CNS  involvement.  The  nullification  of  gravity  in 
orbit  affects  more  than  sensory  inputs.  Some  hypothesize  that  the  large  fluid  shift 
within  the  bodies  of  crewmembers  in  orbit  might  directly  increase  cerebrospinal  fluid 
pressure  or  chemically  change  its  composition,  triggering  a  response  from  the  CNS 


3A1so  known  as  sensory  conflict. 


(28:219).  Also,  similar  changes  in  the  pressure  or  composition  of  the  fluid  within  the 
vestibular  organs  could  cause  a  similar  result  (28:219). 

2.3  Motion  Sickness  Treatments 

Recognizing  the  probable  role  of  the  CMS  in  motion  sickness,  researchers  have 
developed  three  main  types  of  treatments:  biofeedback,  desensitization,  and  drug 
therapy.  To  date,  no  method  has  been  completely  satisfactory  and  each  has  specific 
drawbacks.  Biofeedback  trains  a  person  to  use  relaxation  techniques  to  override  the 
autonomic  (involuntary)  nervous  system  and  control  the  physiologic  symptoms  of 
motion  sickness  as  he  or  she  experiences  them.  This  approach  has  had  positive  ef¬ 
fects  for  some  airsick  flight  trainees  but  can  be  very  lengthy  and  requires  the  person’s 
complete  cooperation  and  concentration  (25: 106 )( 1 8: 1 1 53, 1 157)  .  To  desensitize  or 
protectively  adapt  someone  to  SMS,  they  are  subjected  to  graded  levels  of  the  visual 
and  acceleration  stimuli  believed  to  induce  SMS.  The  Space  Science  Board  of  the 
National  Research  Council  reported  in  1987  that  preadaptation  hail  not  been  effec¬ 
tive  against  SMS,  although  it  did  improve  the  effectiveness  of  biofeedbacl;  training 
(25:106). 

Researchers  have  also  recognized  that  motion  sickness  symptoms  can  some¬ 
times  be  controlled  by  drugs.  Astronauts  and  cosmonauts  have  taken  several  differ¬ 
ent  drugs,  most  of  which  are  used  to  treat  terrestrial  motion  sickness,  to  attempt 
to  prevent  or  reduce  SMS  symptoms.  These  drugs  include  scopolamine,  an  anti¬ 
cholinergic  agent4  which  reduces  certain  inter-neuron  communications  within  the 
autonomic  nervous  system,  and  meclizine,  an  antihistamine  which  has  inherent  an¬ 
ticholinergic  effects  (25:106)(30:29).  However,  each  of  the  drugs  has  undesirable  side 
effects  varying  from  sedation  to  blurred  vision.  These  side  effects  pose  their  own 
operational  hazards  which  lead  researchers  to  combine  drugs  in  an  attempt  to  offset 

4Class  of  drugs  which  oppose  the  transmission  of  aoetycholine,  a  chemical  which  carries  nerve 
impulses  across  intercellular  gap9. 


8 


these  undesired  effects.  A  typical  example  is  scopdex,  a  mixture  of  scopolamine  and 
amphetamine  (30:29).  No  drug  therapy  has  yet  been  found  which  shuts  down  SMS 
symptoms  and,  unfortunately,  astronauts  continue  to  experience  episodes  of  emesis 
while  on  these  medications  (19:2). 

In  an  effort  to  locate  a  more  effective  means  of  preventing  or  treating  SMS. 
AFIT  in  1982  first  developed  a  means  of  measuring  the  physiologic  symptoms  of 
motion  sickness  as  they  progressed.  One  of  the  measurements  used  was  the  eloc- 
troen  phalogram  (FEG),  w-hich  detects  and  displays  the  electrical  activity  of  the 
brain.  Despite  some  earlier  discounting  of  EEC  changes  during  motion  sickness 
(33:63),  the  AFIT  research  team  found  a  consistent  pattern  of  high-power,  very-low- 
frequency  electrical  activity  which  grew,  waxing  and  waning,  as  symptoms  developed 
(5:1022, 1024)(44:39-40.43-44).  The  pattern  is  shown  in  Figure  1. 

The  activity  was  located  in  a  band  of  frequencies  lying  below  4  hertz  which 
is  commonly  called  the  delta  and  sub-delta  frequencies  (6:111).  The  research  team 
recognized  that  although  the  subjects’  external  symptoms  were  those  common  to 
motion  sickness,  their  brain  activity  patterns  resembled  that  of  a  partial  (focal) 
epileptic-type  seizure.  (5:1024) 

Because  the  pattern  was  so  distinctive,  AFIT  began  testing  the  ability  of  the 
classic  anticonvulsant  drug  phenytoin  to  suppress  the  activity  and  perhaps  reduce 
the  subject’s  motion  sickness  symptoms.  They  discovered  that  phenytoin  therapy 
provided  a  mean  symptom-free  time  that  was  11.9  times  longer5  than  that  experi¬ 
enced  under  the  placebo  (5:1074).  In  addition,  none  of  the  experiment  participants 
experienced  any  side  effects  other  than  slight  "light-headedness”  or  a  sensation  of 
enhanced  alertness  (5:1023).  Phenytoin  treatment  also  delayed  the  appearance  of 
the  characteristic  EEG  pattern  in  most  of  the  subjects,  with  one  of  the  participant  s 
EEGs  never  registering  the  onset  of  the  powerful  delta  pattern  (5:1024).  While 

4  With  a  standard  deviation  of  6.2. 
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Frontal 


Figure  1.  Brain  Activity  Pattern  During  Motion  Sickness.  Vogen  reported  this 
signal  propagation  pattern  which  increased  in  signal  power  with  time: 
1)  Initial  symptoms  brought  a  left  parietal  focus;  2)  Almost  simultane¬ 
ously.  or  shortly  thereafter,  there  was  an  ipsilateral  spread  (on  the  same 
side)  to  the  front o-temporal  region;  3)  A  contralateral  spread  (to  the 
opposite  side)  follows;  4)  More  pronounced  symptoms  brought  a  right- 
temporal  focus,  sometimes  unlocalized.  (44:48) 
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the  results  of  the  phenytoin  trials  were  very  encouraging,  they  only  emphasized  the 
question  of  how  motion  sickness  could  be  related  to  the  energetic  delta- waves. 

To  further  study  the  pattern  of  delta-wave  propagation,  AFIT  has  now  turned 
to  a  different  form  of  brain  topographic  mapping,  the  toposcope.  The  following 
chapter  introduces  the  field  of  brain  mapping  and  three  major  mapping  techniques, 
as  well  as  explains  the  choice  of  the  toposcope  for  this  next  phase  of  AFIT  motion 
sickness  research. 


11 


Ill.  Brain  Mapping:  A  Summary  of  Current  Knowledge 


The  field  of  brain  topographic  mapping,  also  railed  brain  mapping,  exists  in  or¬ 
der  to  expand  the  knowledge  of  how  the  human  brain  functions— and  malfunctions. 
It  does  so  by  pioviding  information  about  the  location  of  brain  electrical  activity 
and  the  spatial  and  temporal  relationships  between  separate  locations.  Thus,  brain 
topographic  mapping  is  true  to  its  name  both  the  prefix  fopoand  the  word  map  con¬ 
vey  the  idea  of  “location."  More  importantly,  brain  mapping  can  offer  a  significant 
advance  over  the  usual  methods  used  to  describe  certain  brain  phenomena. 

Biain  mapping  is  a  term  describing  a  group  of  display  methods  for  which  EEG 
signals  are  often  used  as  input.  The  EEG  is  generated  by  brain  cells  as  they  alter 
their  states  of  electrical  charge  by  communicating  with  each  other  via  signals  of 
constantly  shifting  amplitude  and  frequency  (41:1).  These  electrical  changes  extend 
to  the  scalp,  where  electrodes  cau  detect  the  very  minute  voltages  they  create.  As 
these  voltages,  or  potentials,  continuously  change,  they  can  be  traced  onto  strips 
of  paper  to  form  extremely  complex  graphs  similar  to  a  polygraph  (27:168)  (see 
Figure  2).  In  addition,  the  “brain  print’’  of  these  constantly  changing  signals  is 
different  in  each  individual,  as  distinctive  as  the  person's  signature  (46:55). 

These  rapidly  changing  lines  hold  information  about  more  than  the  surface  of 
the  brain.  The  EEG  is  the  result  of  both  spontaneous  changes  in  brain  state  and 
function,  and  the  elicited  changes  resulting  from  a  stimulus,  such  as  hearing  a  noise 
(2l:309)(6:i).  The  tracings  of  an  EEG  can  even  be  used  to  infer  the  location  of  signal 
sources  in  a  third  dimension,  inside  the  volume  thv>-  the  brain  occupies  (26:704).  Yet, 
interpreting  the  graphs  of  an  EEG  with  the  human  eye  is  often  considered  to  be  more 
ofe'"'  n't  than  a  science,  for  the  EEG  contains  too  much  data  in  an  unsuitable  form 
5  .’.s  ' a’,  analysis  (48:1 )  and,  in  all  likelihood,  a  great  deal  of  noise.  This  is  especially 
true  ui  identifying  brain-related  diseases  and  conditions  caused  by  subtle  alterations 
of  '^flekgv  and  activity,  such  as  mental  retardation  ( 10.455)  or  motion  sickness.  It  was 
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Figure  2.  Example  of  Eight  Channels  of  EEG  Tracings  (35:73). 

because  of  this  data  presentation  problem  that  researchers  developed  the  techniques 
of  brain  topographic  mapping. 

Brain  mapping  is  an  ever  broadening  field  as  new  technologies  become  avail¬ 
able.  Upon  encountering  the  words  "brain  mapping,"  the  techniques  that  might 
first  come  to  mind  are  those  that  produce  three-dimensional  images:  computed  to¬ 
mography,  nuclear  magnetic  resonance  imaging,  and  positron  emission  tomography.1 
However,  these  technologies  are  considered  brain- imaging  methods,  as  opposed  to 
brain- mapping.  The  field  of  brain  topographic  mapping  is  limited  to  those  techniques 
applied  to  data  measured  using  the  electroencephalogram  (EEG)  or  the  magnetoen¬ 
cephalogram  (MEG).  This  chapter  will  discuss  only  topographic  mapping  methods, 
and  specifically  only  those  applied  to  interpreting  spontaneous  EEGs,  which  AFIT 
used  to  record  all  of  its  research  data.  Although  brain  mapping  techniques  are  also 
used  to  study  short  bursts  of  brain  activity  which  are  induced  by  specific  visual, 
auditory,  or  sensory  stimuli2  (16:52),  these  applications  will  not  be  addressed. 

'These  methods  are  probably  better  known  by  their  respective  acronyms:  CT-scan,  NMR  or 
MRI,  and  PET-scan. 

2This  method  of  recording  is  known  as  evoked-potential,  or  event-related  potential  recordings. 
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This  chapter  will  overview  several  areas.  It,  will  begin  with  an  explanation  of 
the  role  of  brain  mapping  as  a  research  tool.  This  will  be  followed  by  a  survey  of 
three  groups  of  brain  mapping  techniques  which  appear  in  the  current  literature: 
the  Walter-Shipton  class  of  toposcopes  (29:14),  contour  mapping,  and  statistical 
mapping.  Each  discussion  will  include  a  description  of  the  technique  and  its  strengths 
and  limitations.  Finally,  the  conclusion  will  analyze  the  various  mapping  methods 
with  respect  to  the  research  problem  and  the  goal  of  continuing  the  investigation  of 
brain  nctivir  during  motion  sickness. 

3.1  The  Role  of  Brain  Topographic  Mapping 

A  mapping  pioneer,  Dietrich  Lehmann,  compared  brain  mapping  to  using  a 
subway  map  (48:1).  He  commented  that  just  as  a  verbal  description  of  the  subway 
system  would  be  greatly  aided  by  a  graphical  representation  of  the  various  routes, 
so  too  would  raw  EEG  data  be  complemented  by  a  brain  map.3  Topograp  ic  maps 
do  not  add  new  information,  but  are  intended  to  make  readily  available  the  spatial 
data  locked  inside  the  EEG  tracings  by  showing  the  data  in  a  space-oriented  form 
(22:29)(24:548).  This  is  accomplished  by  regarding  the  electrodes  as  distributed  in 
two  dimensions  over  the  surface  of  the  skull.  The  way  these  signals  are  distributed 
over  this  surface  is  the  key  to  gaining  more  information  as  to  their  origins  (6:108). 

The  field  of  brain  topographic  mapping  continues  to  grow  with  technology  and 
has  found  a  definite  place  in  brain  research.  H.  Petsche  said,  “The  overwhelming 
invasion  of  the  EEG  by  mapping  methods  demonstrates  that  scientists  dealing  with 
the  EEG  have  become  aware  of  the  fact  that  traditional  electroencephalography 
has  neglected  an  essential  aspect,  namely  location''  (29:15).  The  earliest  efforts 
in  the  brain  mapping  area,  dating  back  to  the  1940s  and  early  1950s,  have  given 
way  to  computer-based  methods  which  offer  unique  characteristics  to  the  search 

3Lopes  da  Silva,  along  with  several  other  experts,  stresses  the  importance  of  not  forsaking  the 
raw  EEG  as  the  primary  record  (24'549). 
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for  more  insight  into  brain  functions.  However,  the  first  device  to  be  discusseci, 
the  Toposcope,4  is  especially  unique  due  to  its  foundational  role  in  the  field,  the 
perspective  it  offers,  and  its  relatively  recent  reappearance  in  the  literature. 

3.1.1  The  Toposcope.  The  first  brain  mapping  device  used  exclusively  on  hu¬ 
mans  was  built  by  Grey  Walter  and  Harold  Shipton  in  the  late  1940s  and  called  the 
Toposcope  (45:283).  Its  purpose  was  to  study  the  intrinsic  rhythms,  or  frequencies, 
of  the  brain's  interna!  communications  created  by  the  often  rhythmic  waveforms  of 
electrical,  intercellular  activity  (37:659).  Walter's  theory  was  that  these  rhythms 
were  produced  by  groups  of  millions  of  neurons  firing  in  harmony  and  could  be  tied 
to  brain  function.  In  addition,  if  properly  displayed,  they  would  show  interdepen¬ 
dence  between  separate  areas  of  the  brain  as  messages  were  propagated  from  spatially 
separated  groups  of  neurons  (46:62)(29:14).  He  shared  the  vision  of  Sir  Charles  Sher¬ 
rington  that  the  brain  was  ‘“an  enchanted  loom  where  millions  of  flashing  shuttles 
weave  a  dissolving  pattern,  always  a  meaningful  pattern  though  never  an  abiding 
one’"  (47:14).  Moreover,  he  hoped  to  design  the  toposcopic  display  to  eliminate  the 
interference  of  signals  in  which  he  had  no  interest,  while  highlighting  the  parts  of  the 
brain  where  activity  was  rela^d  in  frequency,  as  if  distinguishing  “from  the  gossip 
backchat  of  bystanders  and  the  welter  of  routine  traffic”  (46:62). 

3.1.1. 1  How  the  Toposcope  Works.  The  Toposcope  was  built  in  several 
ronfiguiations,  but  the  third  model  was  implemented  using  the  then  current  vacuum 
tube  technology’  with  22  miniature  cathode  ray  tubes  arranged  in  the  same  pattern 
as  the  scalp  electrodes  (39:219)  to  display  information  tying  together  the  amplitude 
and  frequency  of  the  brain's  activity  at  each  electrode’s  location.  Within  each  display 
tube,  a  vector  would  rotate  at  a  controlled  rate  like  the  hand  of  a  clock  about  the 

4  "Toposcope"  will  be  capitalized  in  this  report  only  when  referring  to  a  device  built  by  Walter 
and/or  Shipton.  The  lower  case  usage  will  refer  to  similar  tools  created  by  others  or  to  the  class  of 
devices. 
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tube’s  center  (29:14),  looking  much  like  a  World  War  II  radar  display.  A  photograph 
of  the  display,  taken  in  the  1950s,  appears  in  Figure  3. 


Figure  3.  The  Toposcope  Display  of  Walter  and  Shipton.  The  faint  lines  between 
the  tubes  indicate  the  linkage  of  each  tube  to  the  electrodes  on  the  head 
of  the  subject.  The  electrode  locations  appear  as  small  circles  between 
the  displays  (46:63). 

The  brightness  of  th*'  luminous  point  located  at  the  edge  of  the  sweeping  vector 
was  controlled  by  the  signal’s  amplitude,  making  the  display’s  brightness  instanta¬ 
neously  responsive  to  the  amplitude  of  electrical  signal  present  at  that  location  by 
turning  off  or  on  the  electron  beam  tracing  the  display  (46:62).  The  vector  rotation 
rate  was  manually  adjusted,  tuning  the  display  to  a  particular  frequency  component 
of  the  EEG  (46:62). 
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The  vectors  on  all  22  screens  rotated  in  lock-step,  but  the  display  regions  of 
light  and  dark  in  each  screen  varied  because  brightness  was  controlled  in  each  by  a 
different  EEG  channel.  These  displays  could  not  only  differ  from  each  other,  indi¬ 
cating  different  dominant  frequencies  at  spatially  separate  brain  locations,  but  each 
would  also  change  rapidly  over  time,  indicating  the  constantly  changing  frequency 
and  amplitude  modulation  of  the  EEG  signals.  As  quoted  by  Walter  and  Shipton, 
the  cumulative,  mesmerizing  effect  was  summarized  by  Sir  Charles  Sherrington  as 
‘"...a  sparkling  field  of  rhythmic  flashing  points  with  trains  of  travelling  sparks 
hurrying  hither  and  thither”'  (45:282). 

From  this  flashing  field,  Walter  and  Shipton  sought  to  display  relationships 
between  EEG  channels.  The  key  to  these  relationships  the  vector  sweep  rate 
in  the  displays,  Because  the  vector  in  each  screen  rotated  at  the  same  rate,  the 
displays  were  to  present  continuous  snaoshots  of  the  relationship  between  the  22 
EEG  channels  and  the  reference  frequency  set  by  the  sweep  rate. 

3.J.J.2  Strengths  and  Limitations.  The  creators  of  the  Toposcope  found 
that  when  the  sweep  rate  was  set  to  match  the  frequency  of  the  dominant  EEG 
rhythm,  one  or  more  sections  of  the  tubes  would  become  constantly  bright  whenever 
they  reflected  some  harmonic  of  that  ever-changing  rhythm  (29:14).  Examples  of 
the  wedge,  hour-glass,  and  semicircular  patterns  that  resulted  are  shown  in  Figure  4 
under  the  column  “Rotational  Scan”  .  Phase  relations  could  also  be  easily  recognized 
by  comparing  the  relative  positions  of  these  areas  of  brightness  on  different  tubes 
(12:54).  Signals  of  identical  frequency,  but  different  phase,  presented  similar  patterns 
but  rotated  with  respect  to  each  other. 

For  the  first  time,  the  Toposcope  made  information  about  spatial  and  tempo¬ 
ral  variations  in  the  EEG  instantly  and  simultaneously  available  (39:217).  Screens 
displaying  similar  stationary  patterns  indicated  similar  frequency  signals  at  those 
locations.  That  signal  could  then  travel  to  neighboring  regions,  causing  their  screens 


17 


Fundamental 


1st  Harmonic 


2nd  Harmonic 


3rd  Harmonic 


Ralaranca  Signal  Channel  Signal  Rotational  Circular 

(Swaap)  (inianaKy)  Scan  Scan 


Figure  4.  Examples  of  Toposcope  Displays  of  Simple  Signal  Relationships.  In  each, 
the  1-Hz  reference  signal  is  used  to  correspond  to  one  sweep  of  the  rotating 
vector.  Positive  amplitudes  of  channel  signals  create  bright  regions  on 
the  screen  at  locations  corresponding  to  their  position  in  time.  If  these 
frequencies  were  constant,  they  would  result  in  stationary’  patterns  of  the 
same  design. 
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to  show  the  characteristic  pattern.  In  addition,  a  signal  could  be  considered  to  orig¬ 
inate  between  two  adjacent  locations  whose  screens  displayed  the  same  frequency 
pattern,  but  90  degrees  out  of  phase  (45:289).  However,  a  consistently  phase-shifted, 
but  similar,  stationary  pattern  occurring  between  two  separate  areas  of  the  brain 
could  result  from  a  propagation  delay  of  the  same  signal  and  point  to  an  origin  lying 
other  than  between  the  two  electrodes.  The  Toposcope  highlighted  those  frequencies 
harmonically  related  to  the  sweep  frequency  and  filtered  out  others  because  related 
frequcries  created  stationary  patterns  and  unrelated  ones  created  patterns  that 
drifted  around  the  si  reen,  never  remaining  in  one  place  long  enough  to  “summate" 
for  the  eyes  or  film  (38:484)(45:287). 

However,  the  Walter-Shipton  toposcope  did  have  some  limitations. 

•  It  was  difficult  to  interpret.  Walter  wrote,  “When  we  began  to  use  this  machine, 
we  found  the  time  maps  hard  to  understand.  But  gradually  the  new  code  has 
begun  to  penetrate  our  thick  heads  and  much  of  what  was  quite  bewildering 
in  the  ordinary  brain  prints  now  seems  to  be  taking  on  new  form  and  luster" 
(46:62). 

•  It  displayed  only  the  signals  at  the  electrode  locations  since  it  performed  no 
interpolation  between  electrodes,  thus  limiting  its  spatial  resolution  (39:220). 

•  It  was  incapable  of  exact  amplitude  measurements,  but  such  measurements 
were  not  among  its  intended  purposes  (45:285,  287)  (39:218,  220).  Shipton 
considered  this  a  necessary  sacrifice  of  any  area-display  seeking  spatial  and 
temporal  resolution  (38:485). 

•  It  required  photography  as  an  integrating  mechanism  over  several  revolutions 
to  extract,  and  capture  the  stationary  patterns  (46:62)  (45:289). 

•  It  was  difficult  to  set  the  speed  of  the  servomechanism  controlling  vector  ro¬ 
tation  because  of  the  operator's  tendency  to  focus  on  moving  rather  than  sta¬ 
tionary  patterns  and  because  the  frequency  of  the  EEG  is  not  stable  (38:484). 
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3.1. 1.3  The  Updated.  Version  of  the  Toposcope.  Additional  versions  of 
toposcopes  were  created  later  by  Shipton,  with  the  most  recent  being  announced  in 
1981.  In  building  this  latest  toposcope,  Shipton  updated  it  with  a  mixture  of  analog 
and  digital  technology,  but  retained  the  concept  and  type  of  display  and  features 
from  an  earlier,  interim  toposcope  announced  in  1963  (37:659)(38:483).  This  display, 
christened  the  "circular  scan,”0  altered  the  original  rotational  sweep  to  cause  it  to 
inscribe  a  spiral  on  the  screen  by  using  a  luminous  point  positioned  at  the  end 
of  a  rotating  vector  of  increasing  length.  This  added  a  fixed-length  time  axis  to 
the  display,  making  it  possible  for  an  observer  to  see  a  brief  history  of  the  most 
recent  signals  for  comparison.  Figure  4  provides  examples  of  this  display  under 
the  column  "Circular  Scan.”  Signals  harmonically  related  to  the  rotational  speed 
appeared  on  the  new  display  at  the  same  angular  position  but  at  an  increased  radius 
on  each  successive  sweep  (38:484).  Those  signals  differing  slightly  from  a  harmonic 
relationship  would  result  in  a  drifting  pattern  now  more  easily  recognized.  Figure  5 
illustrates  these  features  as  it  compares  the  temporal  and  frequency  resolution  of  the 
initial  and  revised  displays. 

Now  time  could  be  read  over  short  periods  by  measuring  the  distance  of  any 
portion  of  the  spiral  to  the  screen's  center  and  photography  was  only  necessary 
as  a  recording  medium  and  not  an  integration  mechanism  (38:484).  In  addition, 
the  circular  scan  system  could  more  easily  differentiate  between  signals  of  slightly 
differing  frequency  (39:220).  The  time  period  of  the  display  was  a  maximum  of  60 
sweeps  (38:485). 

An  additional  feature  of  the  new  Toposcope  was  the  ability  to  use  the  signal 
in  a  selected  channel  as  the  reference  so  that  the  activity  of  the  entire  brain  could 
be  correlated  and  compared  with  it  (38:484).  By  measuring  the  frequency  of  the  ref¬ 
erence  signal  from  the  intervals  between  time-axis  crossings,  that  signal  would  drive 
the  rate  of  vector  rotation  in  the  displays  of  all  channels  (38:490),  thus  eliminating 

5  Also  known  as  the  “helical"  scan  (39:14). 
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Figure  5.  Comparison  of  Resolution  of  Original  and  Updated  Walter-Shipton  Dis¬ 
plays.  A.  The  photographically  integrated  original  display  is  ambiguous 
when  presented  with  short  duration  signals  or  those  slightly  offset  from 
the  frequency  of  interest.  B.  The  circular  display  clarifies  these  ambigu¬ 
ities,  showing  brief  signals  in  the  10  and  12  o’clock  directions,  a  steady, 
harmonic  signal  at  8  o’clock,  and  a  steady,  but  non-harmonic  signal  at  4 
o’clock,  resulting  in  a  constant  shift  in  subsequent  traces  (39:220). 


the  need  to  manually  control  the  rotational  frequency.  The  rotational  speed  could 
be  related  either  to  the  fundamental  frequency  or  to  some  sub-multiple  of  the  signal 
by  use  of  an  attenuator  (38:490).  Although  Waiter  did  mention  that  the  sweep  rate 
on  the  original  Toposcope  could  be  controlled  by  signals  from  the  subject’s  brain 
(46:62),  no  indication  of  how  this  was  done  could  be  found  in  the  literature. 

This  1981  version  of  the  Toposcope  made  several  positive  contributions  but  had 
its  limitations  as  well.  By  virtue  of  the  new  circular  display,  it  did  remove  the  need 
for  photographic  integration  in  order  to  focus  on  stationary  patterns  in  the  display. 
The  addition  of  an  automatic  mode  which  drove  the  sweep  rotation  according  to 
one  EEG  channel  eliminated  the  restrictions  and  difficulties  of  manual  control.  Also, 
Shipton  reported  that  this  system  was  better  suited  than  the  original  Toposcope  for 
studying  spontaneous  EEG,  as  opposed  to  evoked  potentials  (38:485).  However,  its 
limitations  included: 
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»  It  did  not  address  the  limitations  of  its  predecessor  in  spatial  resolution  and 
amplitude  measurement  accuracy. 

•  It  was  unable  to  display  evoked  responses  to  slow,  repetitive  stimuli  (rates 
below  3  per  second)  without  losing  adequate  frequency  resolution  (38:485). 

The  Toposcope  continues  to  offer  a  unique  perspective  on  brain  activity  re¬ 
search,  displaying  instantaneous  relationships  in  frequency  and  phase  between  sep¬ 
arate  areas  of  the  human  brain  (39:221-222).  But  its  Shipton-defined  role  as  “an 
extremely  sensitive  indicating  device,  rather  than  a  precise  measuring  instrument” 
is  an  important  distinction  (38:485). 

3.1.2  Contour  Brain  Mapping.  Ever  since  Antoine  Remind  created  the  first 
chronotopogram  in  the  early  1960s  (35:73-93),  contour  mapping  has  been  an  ef¬ 
ficient  and  popular  way  of  overcoming  the  problem  of  discrete  EEG  sampling  in 
space  (23:145).  Contour  mapping,  unlike  the  Walter-Shipton  class  of  toposcopes, 
uses  interpolation  methods  in  order  to  smooth  out  the  informational  discontinuities 
between  electrodes  (23:143).  The  output  is  an  image  very  similar  to  a  contour  map 
of  a  geographic  area,  with  free-form  curves  indicating  areas  of  similar  elevation  and 
the  spacing  between  curves  showing  an  incremental  elevation  change. 

3. 1.2.1  How  Contour  Mapping  ITorfcs.  There  are  two  basic  groupings 
of  contour  mapping  techniques  (26:705).  One  consists  of  recording  scalp  potentials 
from  a  two-dimensional  array  of  EEG  electrodes  at  a  single  moment  in  time  and 
creating  a  map  of  the  potential  field  using  interpolation  methods.  Thus  a  series  of 
maps  can  be  produced  showing  changes  over  time.  The  second  group  records  only 
from  a  line  of  electrodes  which  measure  the  EEG  potential  at  specific  points  along 
one  dimension.  However,  the  recording  is  continuous,  producing  a  two-dimensional 
plot  with  the  time  axis  running  the  length  of  the  output.  Here,  too,  an  interpolation 
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scheme  is  used  to  estimate  those  points  not  measured.  The  chronotopogram  belongs 
to  this  class  of  contour  mapping  (26:705).  An  example  appears  in  Figure  6. 
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Figure  6.  Rernond’s  Chronotopogram.  An  example  of  the  first  contour  brain  map 
(35:81). 

There  are  several  different  means  of  interpolation  in  use,  both  linear  and  non¬ 
linear.  Although  interpolation  is  a  common  practice,  its  use  should  not  be  considered 
trivial.  Instead,  careful  consideration  should  be  given  to  selecting  the  best  interpo¬ 
lation  method  since  99  percent  of  the  resulting  brain  map  will  result  from  the  al¬ 
gorithm’s  calculations,  rather  than  from  the  measured,  known  points  (27:168).  The 
most  common  linear  methods  used  in  brain  mapping  applications  are  the  nearest - 
three  -  neighbor  and  the  nearest-four-neighbor  algorithms,  according  to  Wong  (48:14). 
They  cause  the  interpolated  values  to  reflect  the  trend  set  by  the  values  of  the  bound¬ 
ing  three  or  four  points.  Variations  of  these  linear  methods  differentially  weight  the 
bounding  or  spatially  distant  points  to  affect  the  calculated  value  (11:21).  In  non¬ 
linear  interpolation,  the  calculation  fits  a  complex  curve  to  the  data  points  and  uses 
the  resulting  equation  to  determine  the  needed  values. 

Sometimes,  averaging  procedures  are  also  used.  Measurements  of  potentials  at 
each  electrode  are  averaged  together  over  some  predetermined  period  of  time  and 
then  these  average  values  are  mapped  as  before  (26:705).  This  results  in  the  loss 
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of  information  on  high  frequency  changes  in  brain  activity,  but  works  to  highlight 
the  most  significant  activity  patterns  since  they  were  present  during  all  or  most  of 
the  averaged  period  (48:12).  A  secondary  benefit  is  the  resulting  data  reduction  in 
generating  only  a  single  average  map  (48:12). 

Contour  maps  are  not  limited  to  displaying  instantaneous  voltages,  but  can 
also  display  voltages  from  portions  of  a  signal's  frequency  spectrum.  The  EEG 
data  can  be  processed  by  a  computer-based  mathematical  technique  called  a  fast 
Fouri'  •  transform  (FFT)  which  transforms  a  signal  stream,  often  1 -to- 10-seconds 
long,  from  the  time  domain  into  the  frequency  domain  (13:348).  The  resulting 
individual  frequency  components  could  then  be  summed  back  together  to  reproduce 
the  original  signal.  These  components  can  be  displayed  for  a  single  electrode  as  traces 
in  which  the  vertical  axis  is  voltage  or  power  and  the  horizontal  axis  is  frequency. 
However,  the  spectrum  region  for  the  entire  array  of  electrodes  can  also  be  displayed 
as  a  contour  map,  affording  an  image  of  the  amplitude  of  overall  brain  activity  in 
that  frequency  band  (48:24). 

Artificially  coloring  the  regions  of  a  contour  map  can  highlight  the  areas  of 
similarity  and  dissimilarity,  and  therefore  is  a  very  popular  option  (27:168).  Shades  of 
gray  or  variations  in  the  density  of  dots  or  other  symbols  are  comparable  highlighting 
methods. 


3. 1.2. 2  Use  of  Contour  Drain  Maps  at  AFIT.  Both  Banducci  and  Yo- 
gen  used  contour  maps  in  analyzing  motion  sickness  EEG  patterns  (2:86)(44:20-21). 
Both  researchers  used  the  contour  maps  created  by  a  software  package  called  the 
Brain  Atlas6  to  display  voltage  maps  from  the  delta  region  of  the  brain’s  frequency 
spectrum.  Examples  of  these  color  maps  are  reproduced  in  black  and  white  in  Fig¬ 
ure  7. 

8  A  product  of  Bio-logic  Systems  Corporation. 
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8.5-13.0 


13.5-23.5 


Figure  7.  Examples  of  Contour  Brain  Maps  Used  in  Prior  Research  at  AFIT.  These 
four  maps  display  activity  in  four  different  frequency  bands,  as  labeled 
by  the  numerical  range  below  each  map. 
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They  processed  14-channel  EEG  data  using  FFTs  to  transform  2-to- 10-second 
epochs  (blocks)  of  data  into  the  frequency  domain  (44:32).  The  resulting  maps  were 
produced  in  color  and  clearly  displayed  the  spread  of  an  electrical  activity  pattern 
as  motion  sickness  symptoms  grew  in  severity. 

3. 1.2.3  Strengths  and  Limitations.  Contour  mapping  offers  two  distinct 
advantages  to  the  field  of  brain  mapping.  First,  it  presents  data  in  a  manner  easily 
assimilated,  even  by  people  who  know  relatively  little  about  EEGs  (27:168).  Second, 
when  used  in  conjunction  with  the  FFT,  contour  mapping  can  present  this  data  in 
either  the  time  or  frequency  domain,  permitting  researchers  to  look  at  brain  activity 
across  all  frequencies  or  by  frequency  region. 

Despite  the  power  of  the  F*-  "_  function  and  the  appeal  of  the  smooth,  colorful 
display,  there  are  several  acknowledged  drawbacks  to  the  contour  mapping  approach. 
The  following  are  precautions  to  consider  when  interpreting  a  contour  brain  map. 

•  The  vast  majority  of  the  mapped  values  are  interpolated,  and  the  type  of  in¬ 
terpolation  determines  how  the  peak  voltages  should  be  considered.  In  linear 
methods,  any  peak  voltages  are  forced  to  lie  at  the  locations  of  the  elect rodps 
rather  than  at  a  location  of  a  calculated  value  (48:14).  In  non -linear  interpola¬ 
tion,  voltage  peaks  can  occur  at  any  location,  but  this  method  can  also  assign 
maxima  where  none  actually  exist  (48:16).  Thus,  it  is  important  to  remembci 
that  all  interpolated  points  are  approximations  (11:21). 

•  High  frequency  changes  in  brain  activity  will  not  appear  if  average  values  were 
used  in  constructing  the  display. 

•  Fast  Fourier  Transform  maps  display  integrated,  not  instantaneous,  values 
(7:342).  Specifically,  the  frequency  band  power  will  correspond  to  the  stan¬ 
dard  deviation  of  the  average  of  all  momentary  maps  over  the  time  period 
(20:58-59).  Shipton  also  states  that  additional  difficulties  and  uncertainty  are 
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introduced  in  the  low  frequency  regions  when  attempting  to  splice  together  a 
reliable  picture  of  the  average  activity  over  time  (39:218-219). 

•  Both  the  number  of  electrodes  and  the  reference  voltage  used  can  drastically 
affect  the  appearance  of  the  contour  map  (11 :23)( 2G:705). 

•  It  is  possible  to  misuse  color  in  the  contour  map.  If  contrasting  colors  are  used 
to  represent,  neighboring  amplitude  ranges,  they  can  lead  to  the  interpretation 
that  an  important  threshold  exists,  when  it  really  does  not  ( 1 1:23)(39:223). 

These  limitations  of  contour  mapping  apply  to  the  previous  research  at  AFIT 
based  on  brain  mapping.  While  Yogen  was  specifically  careful  in  selecting  the 
time  periods  (epochs)  he  used  in  averaging  his  samples  (44:32-33),  the  mapping 
method  used  in  his  and  Banducci’s  efforts  could  not  display  the  frequency-dependent, 
momentary  changes  in  the  EEG  pattern,  but  only  averaged  changes  which  were 
“blurred'’  in  the  time  domain  In  addition,  the  majority  of  the  hundreds  of  dis 
played  pixels  in  any  given  map  were  generated  by  interpolating  only  14  data  points. 
In  order  to  produce  a  time-sequenced  approximation  of  the  propagating  pattern, 
Vogen  produced  an  animation  using  many  separate  maps  averaged  over  overlapping 
periods  of  time  recorded  on  a  video  camera  ^  44:34). 

3.1.3  Statistical  Brain  Mapping.  In  statistical  mapping,  the  image  created 
does  not  represent  the  voltage  or  power  of  brain  activity.  Instead,  it  maps  statistics 
comparing  each  point  of  the  EEG  spatial  field  with  some  reference,  with  the  goal 
of  indicating  the  degree  of  difference  among  images.  The  introduction  of  statistical 
techniques  has  greatly  changed  the  outlook  for  brain  topographic  methods  (39:222). 
Gevins  further  emphasized  the  role  of  this  approach: 


Modern  information  processing  of  EEG  recordings  brings  3  essential  ad¬ 
vantages:  'precision  (and  accuracy)  of  measurement, ’  'speed  of  process¬ 
ing,’  and  ‘significance,’  that  is,  a  new  understanding  of  the  results  with 
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the  help  of  statistical  tools.  Any  statement  which  possesses  great  signifi¬ 
cance  must  derive  from  a  large  enough  number  of  accurate  measurements 
and  be  made  after  a  statistical  operation  accompanied  by  a  strong  enough 
statistical  ‘weight.’  (34:vin) 

3. 1.3.1  How  Statistical  Mapping  Works.  Statistical  mapping  usually 
begins  with  a  plot  of  the  brain  potential  field  in  some  or  all  frequencies  produced  by 
a  contour  mapping  method.  It  then  performs  a  pixel-bv-pixel  statistical  comparison 
with  a  reference  map  which  can  be  another  individual  image  or,  mort  commonly,  a 
map  of  mean  values  from  some  population.  The  results  are  then  transfoimed  into  a 
new  image  via  a  statistical  to  optical  transformation  (10:455) 

In  the  case  of  significance  probability  mapping  (SPM),  Duffy,  et  al.  use  two 
different  type  of  statistics:  the  Student’s  t-test  and  the  2-transform  to  delineate 
regional  topographic  differences  in  brain  activity  (10:456-457).  In  the  ^-statistic 
form  of  SPM,  the  individual  points  are  summed  between  multiple  maps  within  two 
separate  sets:  the  control  and  experimental  sets.  The  results  are  two  maps  consisting 
of  the  summed  pixel  values  of  their  respective  sets.  Each  of  these  images  are  then 
converted  to  a  map  of  the  mean  values  and  a  map  of  the  variances.  Finally  the 
t-test  is  applied  between  the  control  and  experimental  sets  and  the  result  is  the  t- 
statistic  map,  which  reveals  the  regions  in  which  the  differences  between  the  two 
populations  are  statistically  significant  (10:457).  In  2-statistic  SPM,  an  individual 
map  is  compared  against  the  mean  and  variance  maps  of  reference  population.  The 
2-transform  map  results,  where  the  values  are  the  number  cf  standard  deviations  the 
individual  map  lies  from  the  mean  of  the  reference  group  (10:457). 

Another  approach  to  statistical  mapping  is  the  neurometries  method  devel¬ 
oped  by  E.R.  John,  et  al.  which  estimates  the  probability  that  the  quantitative 
measurements  of  brain  activity  reflect  dysfunction  (31:153).  Neurometries  first  ex¬ 
tracts  particular  features  from  the  EEG  data.  This  method  then  produces  brain 
maps  using  the  2-transform  with  the  scale  representing  the  probability  that  the  val¬ 
ues  shown  might  be  found  in  a  healthy,  normally  functioning  person  of  the  same  age. 
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John  has  also  demonstrated  the  capability  of  neurometries  to  accurately  classify  the 
transformed  individual  maps  of  psychiatric  cases  according  to  the  disorder  present 
(17:116). 


3. 1.3. 2  Strengths  and  Limitations.  As  Remond's  quote  pointed  out  ear¬ 
lier,  the  strength  of  the  statistical  methods  lies  in  their  power  to  highlight  the  sig¬ 
nificant  differences  in  brain  activity.  Another  advantage  is  the  capability  to  then 
dependably  classify  maps  into  categories  of  dysfunction.  All  this  increases  the  use¬ 
fulness  of  brain  mapping  as  a  research  tool. 

The  limitations  of  the  method  are: 

•  “Normal”  is  a  relative  label.  There  is  significant  debate  over  defining  the 
statistical  distribution  of  the  normal  EEG.  Different  techniques  sometimes  use 
substantially  dissimilar  limits  of  normal  (27:171). 

•  Medications,  age,  gender,  and  handedness7  affect  the  statistical  comparison. 
These  factors  must  be  considered  because  they  affect  the  EEG  and  could  in¬ 
validate  any  comparison  to  a  normal  population  (27:171). 

•  Abnormal  results  are  not  necessarily  significant.  A  computer  can  easily  gener¬ 
ate  5000  or  more  statistical  results.  Thus,  a  few  results  outside  even  a  three 
standard  deviation  criterion  for  abnormality  could  be  expected  due  to  random 
chance  (27:171). 

•  Statistical  mapping  requires  a  sizable  normative  database  of  maps,  a  signif¬ 
icant  investment  of  resources  (11:25).  Location-by-location  comparisons  be¬ 
tween  individual  maps  cannot  provide  information  about  the  significance  of 
any  differences  because  no  variances  are  available.  Instead,  individual  maps 
should  be  compared  only  against  reference  populations  (21:345). 

7  Whether  the  person  is  right-  or  left-handed. 
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3.2  Further  Application  of  Brain  Mapping  at  AFIT 

Each  group  of  brain  mapping  methods  has  distinctive  strengths  and  limitations. 
Because  of  this,  they  would  seem  to  complement  each  other:  a  toposcope  displays 
instantaneous  timing  and  phase  relationships,  contour  mapping  “fills  in  the  gaps” 
between  electrodes  and  provides  spectral  maps,  and  statistical  mapping  calculates 
the  significance  of  differences  between  maps. 

However,  with  respect  to  AFIT's  quest  to  learn  more  about  brain  activity  pat¬ 
tern'  luring  motion  sickness,  each  technique  may  or  may  not  be  practical.  Contour 
mapping  of  FFT-generated  spoctrums  was  the  method  used  in  the  initial  discovery 
of  the  pattern  of  interest.  Further  application  of  this  technique  does  not  appear  to 
offer  the  most  hope  of  additional  insights.  Statistical  mapping,  while  a  very  powerful 
tool,  requires  a  sizeable  reference  database  which  AFIT  neither  has  nor  will  be  able 
to  afford  in  the  near  future.  A  toposcope  like  that  of  Walter  and  Shipton  .11  reveal 
more  of  the  brain  activity  pattern  which  motion  sickness  “weaves.”  It  is  especially 
suited  to  investigate  instantaneous  spatio-temporal  images  and  should  be  tunable 
to  the  delta/sub-delta  frequency  region.  Although  less  complex  in  implementation 
than  contour  mapping  with  its  interpolation  and  FFT  operations,  the  toposcopic 
technique  will  nevertheless  offer  a  continuous,  alternative  display  of  brain  activity- 
patterns  as  they  spread  over  the  brain's  surface  moment  by  moment. 

Both  Walter  and  Shipton  saw  the  toposcope  as  more  of  a  qualitative  rather  than 
a  highly  quantitative  tool,  stating  that  encephalographers  are  usually  more  concerned 
with  the  frequency  and  spatial  location  of  signals  than  their  exact  amplitude  (38:483). 
Shipton  wrote  in  1986, 


Numbers  are  not  (Lord  Kelvin  notwithstanding)  the  only  route  to  knowl¬ 
edge.  An  analogy  can  be  found  in  the  analysis  of  handwriting.  It  is  ex¬ 
ceedingly  difficult  to  represent  a  signature  by  a  polynomial;  a  bank  clerk 
can  quickly  recognize  a  vaiid  or  invalid  signature  even  in  the  presence  of 
considerable  artifact  (39:222). 
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Given  th«  dual  qualitative/quantitative  nature  of  AFIT’s  search  for  motion 
sickness  brain  activity  patterns,  the  desire  to  use  a  different  topographic  method 
to  compare  against  the  pievious  applications  of  contour  mapping,  and  the  lack  of 
the  sizeable  database  required  to  implement  statistical  mapping,  the  toposcope  was 
chosen  as  the  tool  for  this  latest  phase  .of  motion  sickness  research.  The  AFIT 
toposcope  will  function  similarly  to  that  of  Walter  and  Shipton.  Furthermore,  while 
the  limitations  of  spatial  and  amplitude  resolution  are  characteristic  of  the  toposcope 
class  of  devices  and  will  not  be  improved  upon  by  AFIT's  new  toposcope,  the  use 
of  computer  code  instead  of  electronics  as  the  display  s  building  blocks  offers  new 
flexibilities,  features,  and  interpretation  aids.  How  this  tool  was  implemented  is 
described  in  the  next  chapter. 
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IV.  Research  Approach 


There  were  three  major  efforts  which  comprised  this  research.  The  first  was  to 
access  and  convert  the  binary  EEG  data  in  the  Bio-logic  computer  files  in  the  Motion 
Sickness  Lab.  The  next  task  was  to  design  and  construct  a  toposcope  by  writing 
software  using  C  to  drive  a  computer  graphics  terminal.  The  last  task  was  to  run  test 
signals  and  digitized  EEG  data  through  the  toposcope,  making  adjustments  to  the 
display  and  frequency  selection,  in  order  to  characterize  the  toposcopes  performance 
and  study  the  brain  activity  pattern  in  the  delta  frequency  band.  This  chapter  will 
describe  the  approach  used  in  each  of  these  areas  as  well  as  list  the  required  material 
and  equipment 

4. 1  Preprocessing  the  Data 

The  first  task  was  to  ensure  digitized  EEG  records  would  be  available  to  drive 
the  computer-based  toposcope.  No  additional  recordings  could  be  made  due  to  the 
loss  of  funding  for  AFIT  motion  sickness  research  at  the  end  of  FY91.  Records 
from  1991  were  available  on  14-channel  beta  tapes,  with  the  data  recorded  in  analog 
form.  Each  of  these  records  had  also  been  previously  digitized  using  the  Brain  Atlas 
software  and  stored  on  diskettes  or  a  hard  disk,  but  the  format  of  the  data  within 
each  file  was  unknown. 

The  data  preparation  task  is  documented  in  Appendix  A  and  consisted  of  the 
following  steps: 

•  STEP  1.1  -  Determine  format  of  data  in  Bio-logic  files. 

•  STEP  1.2  -  Validate  file  interpretation  and  voltage  values  obtained. 

•  STEP  1.3  -  Move  data  files  to  Silicon  Graphics  machines  and  convert  to  l' NIX 
format. 
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•  STEP  1.4  -  Read  pertinent  file  header  data,  reorganize  data  by  channels,  and 
read  in  data  points. 

«  STEP  1.5  -  Validate  output. 

4-2  Designing  and  Coding  the  Display 

The  next  task  was  to  design  and  construct  a  toposcope  on  a  computer  graphics 
terminal.  The  steps  in  this  process  follow: 

•  STEP  2.1  -  Learn  basic  C  syntax  and  grammar,  as  well  as  the  fundamentals 
of  4Sight,  the  Silicon  Graphics  windowing  system. 

•  STEP  2.2  -  Obtain  public-domain  graphics  code  as  similar  as  possible  to  the 
desired  product  and  learn  how  it  functions. 

•  STEP  2.3  -  Construct  display  prototype  program  from  sections  of  original  and 
borrowed  code. 

•  STEP  2.4  -  Build  test  data  files  using  known  signals  (sinusoids  in  the  delta 
frequency  range). 

•  STEP  2.5  -  Test  Phase  I  display  using  test  data  files.  Test  display's  ability  to 
highlight  harmonics  of  sweep  frequency  to  a  human  observer. 

•  STEP  2.6  -  Construct  Phase  II  toposcope  based  on  lessons  learned. 

Instead  of  the  circular  displays  of  the  Walter-Shipton  toposcope,  square  dis¬ 
plays  for  each  channel  were  chosen  for  the  AFIT  version  because  of  the  relative  ease 
of  creating  and  updating  moving  lines  on  the  Silicon  Graphics  workstations.  The 
display  for  a  given  channel  consists  of  a  horizontal,  gray-scale  “bar"  of  fixed  dimen¬ 
sions  which  moves  up  and  down  wit  hin  the  region  of  a  box.  The  boxes  appear  on  the 
display  in  the  same  relative  orientation  as  their  sources  on  the  subjects’  scalps.  The 
intensity  of  all  bars  is  determined  by  a  single  signal  or  channel,  known  as  the  refer¬ 
ence,  while  the  vertical  displacement  of  each  bar  within  each  of  the  14  channel  boxes 
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is  driven  independently  by  the  signal  of  that  channel.  This  results  in  the  brightest 
bands  appearing  in  all  channels  whenever  the  amplitude  of  the  reference  signal  is 
peaking,  but  at  locations  varying  across  channels  due  to  the  different  instantaneous 
amplitudes  of  the  EEG  in  each  channel.  Figure  8  shows  the  original  concept  of  the 
AFIT  toposcope  while  Figure  9  gives  simplified  examples  of  expected  displays  of 
sinusoidal  test  signals. 


Figure  8.  Original  Concept  of  the  AFIT  Toposcope.  Note  that  horizontal  traces 
replace  the  circular  sweeps  of  the  Walter-Shipton  version. 


Harmonic  relationships  between  reference  and  channel  signals  would  be  in¬ 
dicated  by  one  or  more  stationary  bands  of  maximum  intensity  in  that  channel’s 
display  (see  Figure  9).  However,  it  should  be  noted  that  because  of  the  difference  in 
the  role  of  the  reference  signal  between  the  Walter/Shipton  and  original  AFIT  topo- 
scopeS,  different  harmonic  relationships  are  displayed.  The  original  toposcope  used 
the  reference  signal  as  an  angular  displacement  driver  to  set  the  speed  of  the  rotating 
vector  and  driving  the  display  to  key  on  harmonics  of  that  reference  frequency.  For 
example,  a  stationary  hourglass  pattern  (two  opposing  wedges)  indicated  that  the 
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Figure  9.  Expected  Test  Displays  on  the  AFIT  Toposcope.  Input  signals  are  si¬ 
nusoids.  The  reference  signal  determines  the  brightness  of  the  bar  (only 
the  peak  brightness  is  shown  in  the  simplified  version)  and  the  channel 
signal  determines  the  vertical  position  of  the  bar  within  the  box.  Note 
the  intensity  gradients  occurring  in  the  actual  display.  Displays  shown 
are  not  instantaneous  but  represent  the  composite  display  as  seen  by  an 
observer. 
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channel  signal  was  a  second  harmonic  of  the  reference  signal.  But  in  the  original 
AFIT  version,  the  reference  signal  drives  the  brightness  of  the  displayed  bar  and 
the  channel  signal  drives  the  displacement,  vertically  in  this  case.  If  a  double- bar 
stationary  pattern  results,  it  indicates  that  the  channel  signal  is  potentially  a  second 
sub-harmonic  of  the  reference  signal. 

Thus,  while  both  toposcopes  indicate  channels  whose  signals  are  at  the  refer¬ 
ence’s  fundamental  frequency,  they  also  highlight  signals  on  opposite  sides  of  the  ref¬ 
erence  frequency.  This  role  change  for  the  reference  signal  was  required  by  the  design 
chosen  for  the  AFIT  toposcope  and  hardware  restrictions,  all  of  which  are  explained 
in  Chapter  V.  Also,  a  means  of  automatic  frequency  selection  was  implemented  in 
a  manner  similar  to  Shipton’s  1963  toposcope,  as  discussed  in  Section  3. 1.1. 3. 

4-3  Analyzing  Test  and  EEG  Signals  with  the  AFIT  Toposcope 

The  final  task  was  to  use  the  motion-sickness-affected  EEG  records  as  inputs  to 
the  completed  toposcope  and  analyze  the  resulting  displays  for  evidence  of  frequency- 
based  relationships  between  separate  areas  of  the  brain  as  reported  by  Vogen.  The 
steps  in  this  process  follow: 

•  STEP  3.1  -  Use  test  signals  to  determine  the  toposcope’s  response  to  signals 
sharing  fundamental,  harmonic,  sub-harmonic,  and  phase-shifted  relationships. 
Adjust  the  display  as  necessary'  based  on  lessons  learned. 

•  STEP  3.2  -  Determine  from  Vogen’s  notes  the  time-spans  in  the  EEG  record 
where  the  subjects  began  reporting  motion  sickness  symptoms  and  in  which 
channel(s)  the  high-power,  low-frequency  EEG  pattern  characteristically  ap¬ 
peared. 

•  STEP  3.3  -  Transfer  the  EEG  data  record  to  the  Silicon  Graphics  workstation 
where  the  AFIT  toposcope  was  resident. 
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•  STEP  3.4  -  Display  the  record  on  the  toposcope,  referencing  the  signal  of  all  the 
channels  to  a  fixed-frequency  reference,  as  well  as  to  EEG  channels  determined 
to  be  active  in  the  pattern. 

The  expected  result  was  that  the  display  of  the  reference  channel  would  initially 
show  a  stationary,  single  bar  while  other  channels  showed  no  meaningful  pattern. 
However,  as  the  high  energy,  low  frequency  signal  began  to  propagate  across  the 
brain,  the  corresponding  patterns  would  fade  into  and  out  of  the  reference  channel 
pattern,  perhaps  harmonically  related  and/or  phase-shifted,  as  the  activity  pattern 
waxed,  waned,  and  spread.  The  pattern  was  expected  to  follow  that  described  by 
Yogen  and  shown  earlier  in  Figure  1,  Section  2.3. 

4-4  Chapter  Summary 

This  phased  approach,  accessing  the  data,  constructing  the  display,  and  then 
studying  the  toposcope’s  presentations  of  the  data,  was  intended  to  provide  a  sys¬ 
tematic  means  of  using  a  new  tool  to  further  study  the  brain  activity  patterns  during 
motion  sickness.  While  similar  in  many  ways  to  its  predecessors,  the  AFIT  toposcope 
uses  a  unique  display  methodology.  The  first  step,  EEG  data  record  conversion,  is 
discussed  in  Appendix  A.  How  the  new  toposcope  was  designed  and  coded  to  display 
the  converted  signals  is  the  subject  of  the  next  chapter. 
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V.  Building  the  Display 


The  program  fr  r  the  graphics  display  and  user  interface  of  the  AFIT  toposcope 
grew  from  afounda  ion  formed  by  a  “host”  program  Profile.1  Because  of  the  absence 
of  C  and  graphics  programming  experience,  this  method  of  using  another  program 
as  a  learning  tool  and  foundation  was  adopted  by  the  researcher.  The  program 
chosen  was  the  result  of  an  earlier  AFIT  computer  graphics  project  to  modify  the 
demonstration  program  Curve  Demo  provided  bv  Silicon  Graphics  for  users  as  public- 
domain  software.  Profile  was  chosen  because  it  used  windows,  drew  lines,  ran  on  the 
Silicon  Graphics  workstations,  and  at  least  one  programmer  at  AFIT  knew  how  it 
operated.  The  display  portions  were  eventually  gutted  and  rebuilt,  but  the  event 
handling2  portions  were  left  largely  untouched. 

5.1  From  Code  to  Display 

There  were  several  tasks  in  this  portion  of  the  research:  boxes  and  bars  were  to 
be  drawn  on  the  screen;  the  bars  were  to  move  in  sequence  with  an  incoming  signal; 
and  the  brightness  of  the  bars  was  to  vary  in  accordance  with  a  separate  signal. 

5.1.  J  Drarnng  to  the  Screen.  The  graphics  library  functions  of  the  Silicon 
Graphics  IRIS-4D  workstation  used  for  the  project  made  the  display  of  the  simple 
objects  required  relatively  straightforward,  Functions  existed  to  draw  boxes  as  well 
as  lines  of  varying  widths.  The  bars  used  in  the  display  were  actually  line  segments 
which  were  2  to  30  pixels  wide.  The  vertical  displacement  of  the  bars  resulted  from 
tying  the  y-coordinates  of  the  line  segment’s  vertices  to  the  value  of  data  in  that 
channel,  which  was  clipped  at  three  standard  deviations  from  zero  and  then  scaled 
to  use  the  entire  height  of  the  box.  Initially,  the  intensity  of  the  bar  was  controlled  by 

1  Written  by  Capt  Chuck  Wright. 

2Refers  to  means  of  incorporating  operator  actions  into  program  execution. 
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similarly  clipping  the  reference  channel's  input  and  scaling  it  to  the  256- value  gray¬ 
scale.  Thus  the  bar  would  appear  very  dark  at  a  minima  of  the  intensity  driving 
signal,  very  white  at  a  maxima,  and  some  shade  of  gray  at  points  in-between. 

5.1.2  Animation  Using  Double  Buffering.  Smooth  animation  was  possible 
because  of  the  technique  called  double  buffering.  The  method  involves  displaying 
a  front  buffer  (bit  plane)  while  updating  a  hidden,  back  buffer.  These  buffers  are 
then  switched  and  the  buffer  which  was  earlier  on  display  is  now  updated  in  the 
background.  This  process  continues  as  long  as  required.  Because  the  buffers  can 
only  be  switched  as  fast  as  a  screen  can  be  drawn  or  refreshed,  animation  frequency 
is  limited  to  72  Hz.  Attempting  to  animate  an  object  on  a  single  bit  plane,  even 
a  simple  object  such  as  a  bar,  results  in  a  jerky  and  confusing  display.  Using  this 
double-buffer  technique,  a  previously  drawn  object,  such  as  the  display  bar  in  a 
channel  box,  could  be  erased  from  its  old  position  and  redrawn  at.  a  new  position 
while  in  the  back  buffer  and  then  the  two  buffers  switched,  displaying  the  newly 
drawn  screen.  When  this  operation  was  repeated  rapidly  it  resulted  in  the  bars’ 
smooth  motion  on  the  screen. 

5.2  The  Phase  I  Display 

A  prototype  display  was  coded  and  produced  a  single  window  which  contained 
two  boxes,  corresponding  to  two  EEG  channels.  Within  each  of  these  boxes  was  a 
single,  horizontal  bar  whose  deflection  would  be  controlled  by  the  frequency  of  one 
EEG  channel  while  the  brightness  of  the  bar  would  be  determined  by  the  amplitude 
of  a  second  channel.  If  the  two  channels  were  of  equal  frequency,  a  single  band  of 
brightness  should  be  stationary  in  the  box.  If  the  signal  driving  the  brightness  was 
a  harmonic  of  the  signal  driving  the  deflection,  then  more  than  one  stationary  band 
would  be  displayed.  In  addition,  if  the  signals  differed  only  slightly  in  frequency,  the 
brightness  bands  would  slowly  drift  up  or  down  within  the  box. 
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The  initial  display  was  difficult  to  interpret  because  too  much  information  was 
being  displayed.  The  smoothly,  but  rapidly  changing  shade  of  the  bar  as  it  traveled 
the  height  of  the  box  was  interesting,  but  confusing,  making  it  hard  to  track  the 
locations  of  the  brightest  bands  on  the  screen.  A  wider  bar  (20  vs.  5  pixels)  and  the 
decision  to  only  display  the  bar  when  the  intensity  driver  was  at  an  amplitude  peak 
helped  significantly.  Since  the  number,  location,  and  drift  of  the  brightest  bars  was 
the  key  to  interpreting  the  display,  both  of  these  changes  worked  well  in  focusing 
the  v’-'wer’s  attention  on  the  important  aspects  of  the  display.  However,  it  was  this 
change  which  also  determined  the  role  of  the  reference  signal  in  the  display.  Since 
the  display  was  only  drawn  as  indicated  by  the  occurrences  of  peaks  in  the  intensity 
driver,  the  display  would  be  faster  if  all  14  channels  were  updated  with  each  screen 
redraw  event.  This  mandated  that  the  reference  signal,  since  it  was  common  to  all 
channels,  should  be  the  intensity  driver.  To  do  otherwise  would  significantly  slow 
the  display  update  rate. 

With  these  modifications,  the  Phase  I  display  provided  the  displays  expected 
when  driven  with  sinusoidal  test  signals  between  0.5  and  2  Hz.  The  display  suc¬ 
cessfully  identified  first  and  second  harmonic  relationships  between  brightness  and 
displacement  driving  signals.  These  test  signals  were  generated  on  a  function  gen¬ 
erator,  verified  on  an  oscilloscope,  digitized  using  the  Brain  Atlas,  and  converted 
using  the  prep.dataO  function  in  the  toposcope  code.  While  test  inputs  could  have 
been  generated  by  a  C  program  and  directly  fed  into  the  display,  merit  was  seen  in 
using  as  many  as  possible  of  the  same  steps  to  create  the  test  data  as  were  used  to 
preprocess  the  EEG  data  files. 

5.3  The  Phase  II  Display 

While  the  preceding  display  had  established  the  basic  principles  used  in  the 
display  of  the  AFIT  toposcope,  the  Phase  II  toposcope  created  the  tool  which  was 
better  suited  for  the  EEG  data  analysis  which  was  to  follow.  It  was  during  this 
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period  that  the  AFIT  toposcope  was  given  the  name  TOPOS ,  the  Greek  word  for 
“place."  Appendix  I  contains  the  code  for  the  TOPOS  program  with  associated 
comments. 

5.3.1  Appearance  Modifications.  The  two  boxes  in  the  Phase  i  display  were 
replaced  with  14  smaller  ones,  arranged  in  the  pattern  in  which  EEG  electrodes 
where  placed  during  data  collection.  Text  was  added  to  clearly  label  the  channels 
and  the  standard  electiodc  locations  for  reference  purposes.  Status  and  configuration 
information  was  written  to  the  screen  to  keep  the  user  appraised  of  current  data  file, 
elapsed  time,  reference  channel,  and  selected  display  options. 

5.3.2  Display  and  Processing  Modifications.  Vertical  displacement  in  the  dis¬ 
plays  was  no  longer  scaled  to  each  channel’s  maximum  and  minimum  values  because 
extraneous  spikes  in  the  data  caused  the  majority  of  data  points  to  be  compressed 
into  a  narrow  region.  Instead,  the  upper  and  lower  limits  of  bar  travel  were  set  to 
three  standard  deviations  of  each  channel,  using  values  calculated  during  data  file 
conversion. 

A  significant  change  was  the  introduction  of  a  time  axis  in  the  display,  similar 
in  principle  to  Shipton’s  circular  display.  The  display  of  TOPOS  changed  so  rapidly 
that  the  advantage  which  Shipton  found  in  his  circular  scan  over  the  original  rota¬ 
tional  scan  became  very  apparent;  TOPOS  also  needed  a  means  of  visually  tracking 
the  display's  history.  The  time  axis  of  the  circular  scan  was  adapted  to  TOPOS  in 
the  same  linear  fashion  as  its  display:  time  “bins"  were  created  in  the  software  and 
each  call  to  draw  the  bar  wrote  only  a  vertical  slice  of  the  bar  as  wide  as  one  bin  and 
right-shifted  but  adjacent  to  the  slice  last  drawn.  An  example  is  shown  Figure  10. 

Thus  a  time  record  of  the  bar  s  path  in  each  channel  display  appeared  on  the 
screen,  being  erased  only  as  the  display  wrapped  around  in  continuation  of  its  task. 
All  this  permitted  the  observer  to  clearly  note  the  bar’s  path  over  thr-  last  69  draw- 
events  at  any  time  during  the  display’s  operation.  This  time-shifted  display  method 
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Figure  10.  Time-Shifted  TOPOS  Display.  Note  the  visual  history'  of  the  display 
bar’s  path  not  provided  by  the  single-bar  display. 


resulted  in  the  first  clear  indication  of  a  fourth-harmonic  relationship  between  signals 
during  TOPOS  testing. 

Five  other  modifications  followed,  all  targeted  at  improving  the  display’s  abil¬ 
ity  to  focus  on  lower  frequency  components  of  the  EEG.  The  first  was  the  option  to 
trigger  draw  events  at  axis-crossings  (sign  reversals)  of  the  reference  signal,  rather 
than  at  its  amplitude  peaks.  Figure  11  provides  an  example  of  signal-display  rela¬ 
tionships  in  both  the  single-bar  and  time-shifted  displays  using  this  event  trigger. 
This  option  was  required  by  the  need  to  study  the  lower  frequency  components  of 
the  EEG  signals.  While  peak-driven  drawing  was  satisfactory  for  the  smooth,  si¬ 
nusoidal  test  signals,  the  technique  keys  to  the  highest  frequency  components  of  an 
EEG  signal.  This  led  to  the  second  modification:  using  “smoothed’’  EEG  inputs. 
Axis  crossings,  especially  those  of  an  EEG  signal  that  has  been  smoothed,  occurs 
at  the  signal’s  lower  frequencies.  Smoothing  was  necessary  because  energetic,  high 
frequency  signal  components  were  often  large  enough  in  amplitude  to  drive  the  signal 
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Figure  11.  TOPOS  Test  Displays  Resulting  from  Reference  Signal  Axis  Cross¬ 
ings.  Display  intensity  is  muted  in  comparison  to  reference  signal 
peak-generated  displays  and  three  bands  now  indicate  a  first  harmonic 
relationship. 
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across  the  axis  while  the  mean  local  amplitude  of  the  signal  was  still  well  above  or 
below  the  x-axis.  The  relationships  between  the  peaks  and  axis-crossings  of  an  EEG 
signal  and  the  benefits  of  signal  smoothing  are  shown  in  Figure  12. 


ORIGINAL  EEG  SIGNAL 


Figure  12.  Advantage  of  a  Smoothed  EEG  Signal.  Note  the  multiple  axis  crossings 
in  the  original  signal  resulting  from  high  energy,  high  frequency  compo¬ 
nents.  These  crossings  are  absent  in  the  smoothed  version,  leaving  only 
the  axis-crossings  at  the  lower  frequencies. 
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The  smoothing  algorithm  used  was: 


Vi 
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where:  y,  =  the  amplitude  at  time  (sample)  i 

n  =  the  radius  of  the  smoothing  window 


(1) 


The  radius  of  the  smoothing  window  was  determined  heuristically.  An  n- value  of  14 
was  found  to  reduce  the  energy  of  high-frequency  components  sufficiently  to  uncover 
some  of  the  lower  frequency  axis  crossings  of  sample  EEG  signals.  These  smoothed 
signals  were  first  substituted  only  for  the  reference  signal  in  order  to  set  the  timing 
of  the  draw  events  in  the  display.  However,  this  sacrified  amplitude  resolution  when 
using  the  scaled  intensity  option  since  the  reference  signal  determined  the  intensity 
of  the  bars  in  the  display.  The  channel  signals,  which  drove  the  displacements  of  the 
bars,  continued  to  use  the  original  EEG  inputs. 

The  third  modification  was  an  option  to  skip  over  a  selected  number  of  reference 
signal  draw  events  (peaks  or  axis-crossings).  This  lowered  the  effective  frequency  of 
the  screen  updates  while  maintaining  some  control  link  between  the  brain's  actual 
activity  and  the  display's  output. 

To  provide  another  option  in  exploring  the  low  frequency  components  of  the 
EEG  signals,  the  fourth  modification  caused  the  Phase  II  display  to  calculate  the 
mean  sample  number  between  every  two  axis  crossings  of  the  reference  signal.  When 
these  calculated  sample  numbers  were  treated  as  actual  axis  crossings,  the  effect 
was  the  same  as  doubling  the  axis-crossing  frequency.  During  the  testing,  it  was 
found  that  reference/channel  relationships  were  more  easily  seen  in  a  peak-driven 
display  when  the  reference  signal  was  a  first-order  harmonic  of  the  channel  signal, 
producing  a  2-band  display.  It  was  hypothesized  that  by  effectively  halving  the 
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constantly  modulating  frequency  of  the  reference  signal  in  this  manner,  that  a  2- 
band  display  might  result  which  would  be  easier  to  detect  than  the  1-band  display 
of  2  signals  fundamentally  related  to  each  other.  Tests  with  simple  and  summed 
sinusoids  supported  this  approach,  but  it  was  not  as  helpful  in  EEG  signals  analysis. 
Because  of  this,  and  due  to  the  time  delays  caused  by  the  preprocessing  to  calculate 
the  additional  data  points,  this  option  was  later  abandoned. 

The  fifth  modification  was  the  addition  of  computer-generated  low-frequency 
sinusoids  as  reference  signal  options.  Complete  data  records  could  then  be  displayed 
and  their  relationships  to  0.5.  1,  1.5,  and  2-H7  signals  studied. 

Other  options  were  also  added  for  controlling  the  intensity  of  the  plots  on 
the  display.  Intensity  scaling  according  to  reference  signal  amplitude  often  seemed 
more  distracting  than  helpful  and  so  the  user  was  given  the  alternatives  of  clamping 
the  intensity  to  a  maximum  value  or  using  a  bi- level  intensity  scheme  based  on  a 
user-selected  threshold.  This  threshold  option  was  intended  to  address  the  display 
distortions  arising  from  the  limited  resolution  of  the  digital  input  signals:  signals 
crossing  the  axis  will  probably  not  be  sampled  at  the  zero-amplitude  point  but 
at  some  point  very  near  the  axis  crossing  with  some  non-zero  amplitude.  If  the 
signal  crosses  the  axis  at  some  steep  slope,  as  EEG  signals  do,  the  amplitude  of 
the  sample  may  be  significantly  displaced  from  zero.  As  a  result,  digitized  test 
signals  were  producing  ragged  patterns  using  axis-crossing  events  which  were  not  as 
easily  recognized  as  the  smooth  patterns  produced  by  computer-generated  signals, 
which  contained  all  of  the  exact  axis-crossing  values.  Using  this  option,  those  bars 
resulting  from  a  reference  signal  value  which  had  changed  sign  from  its  predecessor, 
but  which  was  greater  than  the  selected  threshold,  was  drawn  in  a  subdued  shade  to 
deemphasize  its  presence.  This  acted  to  focus  the  viewer’s  attention  on  the  bars  in 
the  overall  pattern  which  resulted  from  reference  signal  values  under  the  displacement 
threshold. 
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5.3.3  User  Interface.  Linking  control  over  these  many  display  options  within 
the  AFIT  toposcope  to  user  menus  greatly  simplified  the  task  of  viewing  various  data 
records  repeatedly  while  varying  display  parameters.  Descriptions  of  these  menus 
and  their  associated  functions  appear  in  the  TOPOS  User’s  Guide,  Appendix  G.3.4. 

5.4  Chapter  Summary 

A  computer  program  was  written  in  two  phases  to  implement  the  display 
methodology  selected  early  in  the  research.  TOPOS,  the  resulting  AFIT  toposcope, 
offers  various  display  and  processing  options  which  were  added  as  the  tool  matured 
during  testing.  These  options  allow  the  user  to  vary  bar  intensity  and  width  in  order 
to  enhance  the  patterns  present.  Draw  events  can  be  selected  in  order  to  control  the 
effective  frequency  of  the  reference  signal  and  both  fixed-frequency  signals  and  EEG 
channels  can  be  used  as  reference  signals.  The  results  of  using  TOPOS  to  analyze 
test  signals  and  an  EEG  record  of  motion-sick  individuals  are  the  subject  of  the  next 
chapter. 
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VI.  Results:  The  Toposcope’s  Point  of  View 


Finally,  the  last  task  of  this  thesis  was  using  TOPOS  to  study  test  signals 
and  digitized  EEG  data.  The  test  signal  analysis  was  an  important  first  step  to 
attempt  to  characterize  the  toposcope’s  display  for  various  signal  frequency  and 
phase  relationships  by  using  inputs  with  known  characteristics.  This  would  provide 
the  basis  for  interpreting  the  more  complex  EEG  displays  and  indicate  possible 
problems  in  analyzing  them. 

6.1  Test  Signal  Analysis 

Test  files  were  built  in  which  each  channel  was  a  different  signal;  both  periodic 
signals  and  EEG  signals  combined  with  a  sinusoid  were  used.  Appendix  G.3.4  lists 
the  test  files’  contents  by  channel.  Using  TOPOS,  any  channel  could  be  referenced 
to  any  other  channel,  and  the  relationship  displayed.  Low-frequency  sinusoids,  both 
of  constant  and  varying  frequency,  were  generated  by  the  digitization  procedure  used 
for  the  EEG  signals,  producing  an  off-frequency  signal.  They  were  also  generated 
by  computer  program,  producing  true-frequency  signals.  Later,  display  performance 
was  checked  with  EEG  signals  which  had  been  summed  with  1-Hz  sinusoids  of  various 
amplitudes. 

6.1.1  TOPOS'  Initial  Configuration.  The  time-shifted  display  was  used  be¬ 
cause  of  its  addition  of  a  time  axis  to  the  display.  Seventy  “time  bins”  were  used 
to  provide  the  viewer  with  as  much  history  of  the  patterns  presented  as  the  display 
would  allow.  Screen  updates  were  driven  by  reference  signal  axis  crossings.  \V  ile 
peak-driven  drawing  was  also  a  valid  method  for  displaying  test  signal  relationships, 
it  was  confusing  with  EEG  signals  because  of  their  amplitude  modulation  and  so 
was  not  used.  Bar  width  was  varied  as  necessary  for  clarity.  Narrow  bars  were 
helpful  in  detecting  patterns  in  complex  displays,  as  were  broad  bars  in  simpler  dis- 
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plays.  Bar  intensities  were  not  varied  using  either  of  the  developed  options  because 
the  peak  reference  signal  amplitudes  they  represented  in  one  case  were  valid  only  in 
peak-  driven  TOPOS  displays,  and  in  the  other  case,  their  representation  of  sampling 
resolution  error  was  deemed  more  confusing  than  helpful.  Instead,  bars  were  drawn 
using  maximum  intensity. 

6.1.2  Display  Interpretation.  The  displacement  of  the  displayed  bars  from 
the  central  horizontal  axis  in  each  display  was  afforded  the  most  importance  and 
used  dua  the  sole  feature  in  determining  the  signal's  relationship  to  the  reference  and 
in  search  of  any  pattern  presented  across  multiple  channels.  This  was  justified  be¬ 
cause  the  bars  represented  both  the  reference  signal’s  frequency  and  the  observed 
channel’s  relationship  to  it.  The  reference  signal  frequency  was  represented  on  the 
screen  by  the  timing  of  each  new  bar’s  appearance  simultaneously  in  every  chan¬ 
nel’s  display,  and  the  displacement  of  each  bar  represented  the  scaled  amplitude  of 
that  channel’s  signal  at  that  instant  in  time.  It  was  hypothesized  that  similar  rela¬ 
tionships  between  channel  signals  and  the  reference  would  produce  similar  displays. 
Likewise,  dissimilar  relationships  would  result  in  distinguishable  patterns.  The  re¬ 
sulting  patterns  from  the  low-frequency  sinusoids  fell  into  four  categories,  based  on 
the  signals’  relationship. 

6.1.2. 1  Patterns  of  Fundamental  Relationship.  When  the  reference  chan¬ 
nel  and  the  observed  channel  were  of  matching  frequencies,  the  TOPOS  display 
showed  a  single,  horizontal  line  in  the  observed  channel.  This  pattern  was  a  ragged 
approximation  of  the  horizontal  line  if  the  input  signals  were  not  sampled  very  near 
their  axis  crossings,  as  was  often  the  case  with  digitized  signals.  However,  when  the 
signals  differed  slightly  in  frequency,  this  single-line  pattern  split  into  two  interlaced 
and  phase- reversed  sinusoids.  Examples  of  these  three  patterns  appear  in  Figure  13. 
In  addition,  Figure  14  shows  how  two  signals  of  slightly  differing  frequencies  result  in 
the  interlaced  sinusoids  using  draw  events  triggered  by  reference  signal  axis  crossings. 
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Figure  13.  A  TOPOS  Display  of  Fundamental  Relationships.  Pattern  A  resulted 
from  two  signals  of  matching  frequency.  Pattern  B  was  a  ragged  ap¬ 
proximation  of  Pattern  A  due  to  the  signal’s  matching  frequency  but 
insufficient  sampling  resolution.  Pattern  C  arose  when  the  two  signals 
were  nearly,  but  not  exactly,  the  same  frequency. 


Figure  14.  A  TOPOS  Display  Point-by- Point.  This  figure  provides  insight  into 
how  and  why  two  signals  of  approximately,  but  not  exactly,  the  same 
frequency  produce  an  interlaced,  phase-reversed  pattern  on  the  screen. 
The  reference  signal  is  drawn  with  a  dashes.  Vertical  lines  mark  the 
refererence  signal  axis  crossings  and  boxes  mark  the  observed  channel 
signal's  amplitude  at  the  time  of  those  crossings.  These  boxes,  which 
trace  the  diverging  portion  of  the  overall  pattern,  are  what  appear  on 
TOPOS’  display. 
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6. 1.2. 2  Patterns  of  Reference  Signal  Harmonics.  Three  different  pat¬ 
terns  also  emerged  whenever  the  observed  channel’s  frequency  was  an  integer  mul¬ 
tiple  (a  harmonic)  of  the  reference  frequency.  If  the  channel  signal  was  a  perfect 
harmonic  of  any  order,  the  display  was  a  single  horizontal  line,  indistiguishable  from 
that  of  the  true  fundamental  frequency  relationship.  However,  if  the  channel  signal’s 
frequency  was  only  approximately  an  odd  multiple  of  the  reference,  the  interlaced 
pattern  of  near-fundamental  relationship  re-emerged.  When,  instead,  the  channel 
signal  represented  only  an  approximate,  even  multiple  of  the  reference  frequency, 
a  single  sinusoid  appeared  in  that  display,  with  its  period  inversely  related  to  the 
frequency  difference  between  channel  and  reference  signals.  These  harmonic/near- 
harmonic  patterns  are  shown  in  Fig  15. 
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Figure  15.  A  TOPOS  Display  of  Harmonic  Relationships.  Pattern  A  resulted  from 
a  signal  channel  which  was  an  exact  harmonic  of  the  reference  signal.  If 
the  channel  signal  was  a  near,  but  not  exact,  harmonic  and  was  most 
nearly  an  odd  multiple  of  the  reference  frequency,  Pattern  B  was  dis¬ 
played.  Pattern  C  arose  when  the  channel  signal  was  again  off-frequency, 
but  nearly  an  even  multiple  of  the  reference  frequency. 
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6. 1.2. 3  Patterns  of  Reference  Signal  Subharmonics.  Figure  16  presents 
another  trio  of  patterns  produced  when  the  observed  channel’s  frequency  was  exactly 
or  approximately  an  integer  factor  (a  subharmonic)  of  the  reference  frequency.  A 
new  pattern  consisting  of  three  horizontal  lines  was  created  by  channel  signals  which 
were  perfect  first-order  subharmonics  of  the  reference  signal.  The  sinusoidal  pattern 
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Figure  16.  A  TOPOS  Display  of  Subharmonic  Relationships.  Pattern  A  resulted 
from  a  signal  channel  which  was  an  exact  first-order  subharmonic  of  the 
reference  signal.  If  the  channel  signal  was  approximately  a  second  or 
higher-order  subhannonic,  Pattern  B  was  displayed.  Pattern  C  arose 
when  the  channel  signal  was  again  off-frequency,  but  nearly  a  first-order 
subharmonic  of  the  reference  frequency. 


appeared  again  in  the  cases  when  channel  signals  were  approximate  second  or  higher- 
order  subharmonics  of  the  reference  signal.  Another  new  pattern  consisting  of  four 
interlaced  sinusoids  grouped  in  two  phase-reversed  pairs  resulted  from  signals  which 
were  approximate  first-order  subharmonics. 

6. 1.2. 4  Patterns  of  Phase  Shift.  Several  test  signals  in  one  data  set 
were  of  matching  frequencies  but  incorporated  varying  amounts  of  phase  shift.  The 
patterns  they  displayed  when  referenced  to  a  common  signal  indicated  three  cate¬ 
gories  of  phase  shift  effects.  First,  no  change  in  display  would  result  if  the  original 
and  shifted  waveform’s  axis  crossings  were  indistiguishable.1  Second,  if  the  origi¬ 
nal  signals  produced  one  or  more  sinusoidal  patterns  on  the  screen,  the  phase  shift 
caused  a  corresponding  horizontal  shift  in  the  displayed  pattern.  Third,  when  the 
original  pattern  was  one  or  more  horizontal  lines,  the  phase-shifted  signal  caused  a 
vertical  “spreading”  of  the  pattern,  initially  increasing  the  number  of  horizontal  lines 

'For  example,  the  amplitude  at  a  specific  time  sample  of  a  4-Hz  sinusoid  after  a  90-degree  phase 
shift  does  not  change  ftom  the  sample  of  the  original  waveform  at  the  same  time.  This  i9  not  the 
case,  however,  with  a  3-Hz  sinusoid. 
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and  then  increasing  the  distance  between  them.  These  lest  two  types  of  phase-shift 
effects  are  shown  in  Figure  17. 
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Figure  17.  A  TOPOS  Display  of  Phase  Shift.  Series  A  shows  the  leftward  shift  in 
a  sinusoidal  pattern  caused  by  the  negative  phase  shifting  of  one  signal. 
Series  B  show's  the  “spread"  that  occurs  in  a  linear  pattern. 


6.1.3  The  Ambiguity  Problem.  In  the  tests,  nine  categories  of  signal  fre¬ 
quency  relationships  presented  only  five  unique  pattern  types,  and  no  one-to-one 
mapping  existed  by  which  each  pattern  could  point  to  only  one  possible  signal  rela¬ 
tionship.  However,  useful  information  was  still  available  in  the  sense  that  a  funda¬ 
mental  or  harmonic/subharmonic  relationship  was  signaled,  even  if  the  exact  nature 
of  that  relationship  was  unknown.  For  example,  when  TOPOS  displays  a  single- 
line  pattern,  it  indicates  that  the  input  signals  are  of  exactly  the  same  frequency 
or  that  the  observed  channel  signal  is  a  perfect  harmonic  of  the  reference  signal. 
Also,  if  the  pattern  is  the  two  interlaced,  phase-reversed  sinusoids,  the  signals  are 
near-fundamentals  or  the  observed  channel  signal  is  an  approximate,  odd  multiple 
of  the  reference  frequency.  This  ambiguity  had  to  be  considered  when  interpreting 
the  patterns  resulting  from  TOPOS’  analysis  of  EEG  signals.  To  evaluate  TOPOS' 
usefulness  with  EEG  inputs,  the  test  signals  composed  of  EEG  summed  with  a  1-Hz 
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sinusoid  were  referenced  to  signals  fundamentally,  sub-harmonically,  and  harmoni¬ 
cally  related.  The  significant  result  was  when  the  composite  signal  was  referenced 
to  a  1-Hz  sinusoid.  The  pattern  presented  followed  the  same  random  pattern  as 
that  of  the  unaltered  EEG  when  referenced  to  the  1-Hz  signal,  but  was  noticeably 
compressed  towards  the  central  horizontal  axis  in  the  display  box.  It  was  hypoth¬ 
esized  that  the  complex  sinusoidal  patterns  that  were  associated  with  subharmonic 
and  harmonic  relationships  in  the  earlier  tests  were  functions  of  the  periodicity  of 
the  signals  used  and  might  not  be  presented  during  EEG  analysis. 

6.2  EEG  Analysis 

6.2.1  Preparations.  Prior  to  beginning  the  EEG  analysis,  an  initial  criteria 
for  what  constituted  a  significant  pattern  was  needed  and  several  more  changes  were 
required  to  TOPOS  to  better  adapt  it  to  the  differences  between  the  test  signals  and 
EEG  inputs  it  would  now  process. 

6. 2. 1.1  Target  Pattern.  It  was  decided  to  primarily  look  for  compres¬ 
sions  of  the  EEG-produced  patterns  as  signals  of  a  significant  relationship  with  the 
reference  signal.  This  would  be  most  easily  recognized  as  an  attempt  for  the  display 
to  approximate  the  pattern  presented  if  the  reference  signal  was  referenced  to  itself.2 
However,  despite  the  outcome  of  the  tests  using  the  EEG/ 1-Hz  composite  signal, 
spreading  of  the  pattern  sinusoidally  (interlaced  sine  waves)  or  linearly  (multiple 
horizontal  lines)  might  also  occur  and  would  be  considered  important  indications  as 
well. 

The  TOPOS  display  would  be  monitored  for  simultaneous  occurrence  of  either 
pattern:  the  compression  or  spread.  The  resulting  cross-channel  pattern  could  ap- 

2This  pattern  is  available  on  the  TOPOS  display  if  the  selected  reference  signal  is  one  of  the 
14  channels.  However,  if  the  reference  signed  is  one  of  the  four  sinusoidal  signals  which  are  also 
options,  no  display  is  available  for  comparison.  Nonetheless,  the  user  should  remember  that  the 
pattern  would  be  that  of  a  true  fundamental  relationship  with  no  phase-shift:  a  vertically-centered 
straight  line. 
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proximate  the  development  described  by  Yogen  and  shown  earlier  in  Figure  1.  In 
addition,  it  was  hypothesized  that  one  channel  would  drive  the  cross-channel  syn¬ 
chronization  and  when  that  channel  was  used  as  the  reference  signal,  the  activity 
pattern  would  be  most  easily  visible. 

6. 2. 1.2  Recognition  Problems.  The  EEG  records  differed  in  three  sig¬ 
nificant  ways  from  the  majority  of  the  test  signals  used  earlier.  First,  they  were 
constantly  and  rapidly  changing  frequency  under  normal  conditions.  Second,  the 
brain  activity  pattern  which  Yogen  reported  they  held,  was  also  reported  to  contin¬ 
ually  wax  and  wane.  This  meant  that  without  time-averaging  the  data,  as  Yogen 's 
FFT/contour-mapping  method  did,  the  sought-after  pattern  could  be  intermittent 
as  well  as  different  in  shape  from  screen  update  to  screen  update.  The  third  signif¬ 
icant  difference  was  that  the  EEG  signals  were  also  rapidly  amplitude  modulated, 
meaning  that  axis  crossings  would  often  be  made  at  very  high  slopes  and  the  x-axis- 
displacement  resulting  from  the  inadequate  sampling  resolution  could  be  amplified. 
This  would  contribute  to  a  lack  of  smoothness  in  any  resulting  pattern. 

6.2. 1.3  Final  Configuration.  Four  additional  modifications  were  made 
to  the  TOPOS  configuration.  The  first  was  necessary  because  it  quickly  became 
apparent  that  the  time-shifted  display  was  presenting  far  too  much  information  to 
support  the  detection  of  a  pattern  which  might  exist  for  only  a  few  screen  updates 
at  a  time  due  to  an  EEG's  rapid  frequency  modulation.  Therefore,  the  number  of 
time  bins  displayed  on  the  screen  was  made  user  selectable,  permitting  as  few  as  6 
or  as  many  as  74.  The  lower  settings  allowed  the  user  to  focus  his  or  her  attention 
on  a  region  small  enough  to  detect  cross-channel  patterns  which  were  short-lived. 
Figure  18  shows  the  second  change,  the  addition  of  a  correlation  grid  in  the  center 
of  the  TOPOS  screen.  This  was  necessary  to  perform  one  additional  processing  step 
and  better  integrate  for  the  user  the  information  of  TOPOS’  14  channel  displays. 
The  integration  was  required  because  it  was  difficult  for  the  user  to  visually  track  the 


55 


TIME-SHIFT  1  TIME-SHIFT  2  TIME-SHIFTS  TIME-SHIFT  4  TIME-SHIFT  6  TIME-SHIFTS 


TOROS  Correlation  Grid 

Figure  18.  TOPOS’  Correlation  Grid.  Each  of  the  square  regions  corresponds  by 
location  to  an  EEG  channel  in  the  larger  display.  The  six  colored  bars  in 
each  region  represent  three  levels  of  similarity  in  that  channel's  displace¬ 
ment  from  the  reference  signal’s  amplitude  in  the  six  samples.  Black 
indicates  a  very  close  relationship  and  white  indicates  large  differences. 
The  top  row  of  smaller  grids  are  how  the  display  would  look  over  time 
without  the  time-shifting. 
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patterns  in  all  14  channel  displays  simultaneously  because  of  their  spatial  separation. 
The  additional  processing  came  in  the  form  of  the  correlation  grid  showing  each 
channel’s  degree  of  displacement,  not  from  the  x-axis  as  the  rest  of  the  TOPOS 
display  was  showing,  but  from  the  reference  signal's  amplitude  at  that  sample.  This 
was  shown  instantaneously  in  each  region  of  the  correlation  grid  in  the  same  time- 
shifted  format  used  in  the  channel  displays  and  in  three  colors  to  indicate  degrees  of 
agreement  between  the  signals. 

The  third  modification  resulted  from  the  need  to  slow  the  display  to  permit 
the  user  to  analyze  the  constantly  changing  presentation.  Display  update  speed 
was  inversely  scaled  to  the  number  of  time  bins  in  use  and  options  were  added  to 
momentarily  pause  the  display,  as  well  as  to  step  it  one  screen  update  at  a  time,  or 
one  '‘batch’’  (a  complete  update  of  all  time  bins)  at  a  time.  The  fourth  modification 
did  not  involve  software  changes,  but  required  the  use  of  an  option  not  found  helpful 
before.  It  was  necessary  to  use  only  every  seventh  axis  crossing  of  the  reference  signal 
to  trigger  screen  updates  in  order  to  lower  the  effective  frequency  into  the  region  of 
interest  (below  4  Hz)  when  using  an  EEG  channel  as  reference.3  Using  only  every 
third  axis  crossing  was  required  even  with  a  smoothed  EEG  signal  as  reference. 

6.2.2  TOPOS’  Output.  Because  of  the  complexity  of  the  display  resulting 
from  EEG  inputs  and  the  length  of  the  EEG  record,  an  11-second  block  of  data  was 
chosen  for  analysis.  The  data  was  from  Vogen’s  Subject  5  and  began  at  4  minutes, 
53  seconds  into  the  record.  This  same  block  was  singled  out  by  Yogen  as  a  prime 

imple  of  peak  motion-sickness  brain  activity  (44:48-49.70).  His  thesis  included  12 
color  contour  maps  of  6-second  FFTs  made  at  0. 5-second  intervals  (44:70). 

’Because  of  the  EEG's  rapidly  varying  frequency,  the  instantaneous  frequency  calculated  from 
inter-event  times  constantly  varied  and  not  ali  draw-events  were  accompanied  by  frequency  readouts 
of  less  than  4  Hz. 


6.2.2. 1  1.5-Hz  Reference.  The  1.5- Hz  sinusoid  was  selected  as  the  ini¬ 

tial  reference  signal  because  the  Vogen  contour  maps  of  interest  were  displayed  at 
that  frequency.  This  TOPOS  run  was  a  check  for  how  the  EEG  data  block  related 
to  a  fixed  frequency  which  Vogen  had  considered  optimal.  Because  the  reference’s 
frequency  was  in  the  delta  region,  no  skipping  of  draw  events  was  required  to  adjust 
the  frequency. 

The  channel  displays,  provided  in  Figures  19-25,  were  only  sometimes  highly 
correlated4  to  the  reference  signal  as  the  correlation  grid  changed  rapidly.  In 
addition,  each  channel  showed  frequent  episodes  of  activity  that  was  highly  related 
to  the  reference  signal  for  periods  spanning  2  samples.  However,  at  intervals  during 
the  data  block,  certain  channels  were  very  closely  related  to  the  reference  signal  for 
a  larger  number  of  samples.  Figure  2C  shows  the  location  and  frequency  of  highly 
related  channel-reference  activity  as  shown  by  3  or  more  consecutive  black  bars  on  the 

4  As  measured  by  the  correlation  grid. 


58 


TOPOS 

The  AITT 
Interactive 
Topotcope 


'jS? _ , 


S3 


UATCH-DfAW 
CO  S'  FIG I*  RATION 


D:—«^  *-arX., 
'Sl«<%rv«  I 
rr*v*-y  h*.imi 


I^TIh  UWk 

r«M»a*  tw  kj 

i  on* 


Figure  21.  TOPOS’  Display  of  EEG  with  1.5-Hz  Reference  (Part  3). 
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Figure  22.  TOPOS’  Display  of  EEG  with  1  5-Hz  Reference  (Part  4). 
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Figure  23.  TOPOS’  Display  of  EEG  with  1.5-Hz  Reference  (Part  5). 
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Figure  25.  10P0S’  Display  of  EEG  with  1.5-Hz  Reference  (Part  7). 
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Figure  26.  Pattern  of  High  Correlations  on  the  TOPOS  Correlation  Grid.  Each 
black  dot  signifies  three  or  more  consecutive  highly  correlated  ampli¬ 
tudes  between  that  channel  and  the  1.5- Hz  reference.  Parts  1-7  are 
summaries  of  consecutive  displays  spanning  the  same  data  segment  as 
Figures  19-25.  The  intermediate  patterns,  Parts  la-7a,  display  high 
correlations  which  overlap  between  the  original  displays. 


correlation  grid.  The  diagram  also  accounts  for  the  artificial  breaks  in  the  correlation 
grid  due  to  the  still-captures  of  the  updating  display. 

The  correlation  mapping  confirmed  the  constantly-shifting  nature  of  high, 
channel-reference  relativity  with  a  1.5-Hz  reference  signal.  Also,  the  number  of 
highly  correlated  channels  waxed  and  waned  until  Displays  4  and  4a  in  Figure  26. 
At  this  point,  the  number  dropped  sharply,  with  the  low  occuring  during  samples  at 
average  data  record  times  between  5:04.5  (minisec)  and  5:05,  and  then  built  back  to 
a  higher  level  in  Displays  6  -7.  Although  Vogen’s  maps  were  time-averaged,  the  se¬ 
quence  he  provided  shows  the  development  of  the  brain  activity  pattern  he  reported. 
His  pattern  showed  an  average  decline  until  a  low  was  reached  at  approximately 
5:05  before  activity  increased  again.  Thus  the  TOPOS  display  appeared  to  closely 
approximate  the  Vogen’s  pattern.  However,  Vcgen’s  pattern  did  not  then  continue 
to  increase,  as  did  the  correlation  grids  on  TOPOS.  Instead,  the  pattern  repeated 
twice  again  the  build-drop-build  pattern  in  a  rhythmic  manner  prior  to  the  end  of  the 


data  block.  The  TOPOS  correlation  grid  did  not  support  these  final  two  drops,  even 
when  the  criteria  for  meaningful  correlation  was  raised  to  four  or  five  consecutive 
black  bars. 

However,  the  single-sinusoid  and  “spread”  patterns  seen  earlier  in  the  test  dis¬ 
plays  did  occur  in  the  TOPOS  displays  of  the  EEG  record.  Channel  14  in  Figure  22 
and  Channel  5  in  Figure  23  show  the  sinusoidal  pattern  and  and  Channel  11  in  Fig¬ 
ure  19  show's  the  convergence  of  a  sinusoidal  spread  and  Channel  3  in  Figures  22 
and  23  show  both  the  divergence  and  convergence.  Both  of  these  patterns  indicated 
a  near-fundamental,  subharmonic,  or  harmonic  relationship  between  channel  and 
reference  signals.  Furthermore,  these  more  complex  relationships  are  not  considered 
in  the  simple  comparisons  display  of  the  correlation  grid. 

6. 2. 2. 2  EEG  Channel  References.  Changing  the  reference  to  one  of  the 
EEG  signals  now  changed  the  focus  from  the  EEG’s  relationship  to  a  constant  fre¬ 
quency.  to  its  correlation  to  one  channel  within  the  record.  Also,  the  displays  would 
also  be  compared  with  Vogen’s  map  sequence.  Three  channels  were  selected,  based 
on  the  locations  of  earliest  and  highest  activity  in  Vogen’s  contour  maps.  Channel 
6  was  chosen  as  an  example  of  the  initial  left-parietal  activity,  Channel  3  repre¬ 
sented  the  activity  in  the  left  fronto- temporal  region,  and  Channel  11  included  the 
right  temporal  focus  Vogen  noted  (44:48).  Because  the  original,  digitized  EEG  sig¬ 
nals  were  used,  only  every  seventh  axis-crossing  was  displayed  in  order  to  lower  the 
instantaneous  reference  frequencies  into  the  1-5  Hz  range. 

The  resulting  displays  were  analyzed  for  the  three  activity  drops  shown  on 
Yogen’s  maps,  as  well  as  for  any  pattern  of  synchronization  between  the  remaining 
channels  and  the  reference.  All  three  reflected  the  first  drop  in  low-frequency  brain 
activity,  as  did  the  TOPOS  display  using  the  1.5-Hz  reference.  However,  none  of  the 
EEG-referenced  displays  reflected  the  two  later  decreases. 
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Also,  none  of  the  displays  appreciably  increased  their  correlation  levels  over 
that  of  the  1.5-Hz  reference.  None  of  the  three  reference  channels  “captured  the 
attention”  of  the  other  channels  any  more  than  another.  Thus,  no  channel  could  be 
hypothesized  to  be  driving  the  activity,  much  less  have  originated  it. 

6. 2. 2. 3  Smoothed  References.  Smoothed  EEG  signals  were  used  in  two 
scenarios:  1)  only  the  reference  channel  was  smoothed,  and  2}  all  the  EEG  signals 
were  smoothed.  The  first  scenario  was  to  use  smoothing  as  a  means  to  obtain  the 
lower  iiequency  component  of  the  reference  signal  without  distorting  the  remaining 
channels.  The  second  scenario  was  to  obtain  the  lower  frequency  component  while 
distorting  all  channels  somewhat  equally.  These  displays  were  evaluated  similarly  to 
those  in  the  previous  section  and  found  to  perform  approximately  the  same.  One 
difference  was  overall  higher  displacement  correlation  levels  in  the  display  using  all 
smoothed  inputs,  probably  due  to  the  slow  rise  and  fall  of  smooth  curves. 

6.2. 2-4  Other  Observations.  The  most  disturbing  observation  was  the 
lack  of  any  increasing  trend  from  low  to  high  levels  of  correlation  in  the  TOPOS  dis¬ 
play  of  EEG  as  Subject  5  progressed  from  sitting  stationary  without  motion  stimuli 
to  severe  motion  sickness  and  emesis.  While  fluctuations  were  constant,  the  channel- 
reference  correlation  was  higher  at  the  start  of  the  record  when  baseline  EEG  was 
being  ollected.  This  could  be  due  to  the  lower  amplitudes  of  the  resting  EEG  in 
each  channel  being  displayed  on  a  scale  determined  by  the  standard  deviation  of  the 
entire  record. 

It  should  also  be  noted  that  high  correlation  levels  did  not  imply  fundamental 
or  near-fundamental  relationships  between  channel  and  reference  signals.  Such  re¬ 
lationships  could  have  existed,  hut  harmonics  of  the  reference  frequency  would  also 
produce  the  same  pattern. 
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6.3  Chapter  Summary 

The  software  base  of  TOPOS  continued  to  play  a  key  role  in  its  usefulness  as 
an  analysis  tool.  As  a  result  of  studying  the  test  signal  displays  and  in  anticipation  of 
the  different  characteristics  of  EEG  data,  several  modifications  to  TOPOS  then  pre¬ 
ceded  its  application  to  a  motion-sickness-'affected  EEG  data  block.  These  included 
the  ability  to  vary  the  number  of  time  bins  in  the  time-shifted  display  and  to  control 
the  rate  of  screen  updates  in  order  to  better  focus  on  short-lived  patterns.  Also,  the 
correlation  grid  aided  the  observer  by  integrating  in  one  location  the  instantaneous 
correlation  levels  of  axis-crossings  of  channel  and  reference  signals.  Test  signal  anal¬ 
ysis  using  TOPOS  showed  several  categories  of  patterns  that  could  be  produced  by 
various  relationships  between  two  signals.  However,  none  of  them  could  uniquely 
point  to  one  of  the  relationships.  A  block  of  motion-sickness-affected  EEG  data 
studied  earlier  by  Vogen  was  then  studied  using  TOPOS  and  both  fixed- frequency 
and  EEG  channel  references.  The  TOPOS  display  did  not  show  evidence  of  Vo- 
gen’s  pattern,  did  not  identify  any  of  three  key  channels  as  the  pattern's  driver,  and 
did  not  show  a  trend  of  increasing  channel-reference  correlation  as  motion  sickness 
symptoms  began  and  worsened. 
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VII.  Conclusion 


The  last  eight  months  of  work  have  translated  the  Toposcope  of  Walter  and 
Shipton  into  a  multi-featured,  software-based  tool  capable  of  translating  frequency 
and  phase  relationships  between  two  signals  into  a  graphical  display.  All  research 
objectives  were  met:  the  AFIT  toposcope  was  designed  and  built,  test  signals  were 
displayed  and  characterized,  and  motion-sickness-affected  EEG  data  were  analyzed. 

The  toposcope  was  chosen  for  this  research  because  of  the  instantaneous,  un¬ 
embellished  nature  of  its  display.  It  offers  signal  analysis  which  does  not  depend  on 
time-averaging  and  does  not  create  an  entire  landscape  of  electrical  activity  from 
only  a  handful  of  data  points.  The  most  recent  tool  last  used  in  AFIT's  motion  sick¬ 
ness  research,  the  FFT-based  contour  brain  map,  had  both  of  these  characteristics. 
Thus,  the  toposcope  offered  a  different  point  of  view. 

The  design  of  TOPOS  was  evolutionary.  A  vertically  travelling  bar  was  se¬ 
lected  over  the  rotating  vector  of  Walter  and  Shipton  to  better  facilitate  coding,  but 
the  principle  of  their  earlier  devices  was  retained.  That  bar  was  chopped,  shifted, 
narrowed,  and  enlarged  as  the  display  was  changed  to  increase  its  usefulness.  In  its 
final  configuration,  TOPOS  first  notes  the  data  samples  in  which  one  signal,  the  ref¬ 
erence,  crosses  the  x-axis.  The  tool  then  shows  the  user  the  degree  of  displacement 
from  the  x-axis  of  14  different  signals  at  those  samples.  There  are  many  possible 
display  patterns  that  result,  depending  on  the  offset  between  the  axis  crossings  of  the 
refeience  and  displayed  signals,  and  these  offsets  vary  with  the  signals  instantaneous 
frequency  and  phase  shift.  The  patterns  appear  on  the  screen  in  segments  as  the 
display  updates  are  right-shifted  until  the  display’s  boundary  is  reached.  The  user 
can  visually  compare  the  14  channel  displays  for  similarities  to  the  reference  signal’s 
pattern  or  to  other  patterns  which  are  indicative  of  sub-harmonic  or  harmonic  rela¬ 
tionships.  Or,  the  user  can  watch  the  correlation  grid  at  the  center  of  the  display  for 
an  indication  of  the  degree  of  similarity  between  the  channel  and  reference  signals’ 
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displacement  at  each  sample.  A  user’s  first  impression  of  TOPOS’  display  is  well 
characterized  by  Walter’s  description  of  his  own  toposcope:  bewildering  and  difficult 
to  interpret  (46:62).  Vet  watching  the  display  over  a  period  of  time,  coupled  with 
an  understanding  of  how  the  displays  results,  leads  quickly  to  more  familiarity  with 
the  tool  and  the  meaning  of  its  output. 

TOPOS  offers  users  distinct  advantages  over  its  hardware-based  predecessors. 
Because  it  is  constructed  in  computer  code,  TOPOS  gives  users  real-time  control  of 
display  parameters  in  order  to  enhance  and  detect  displayed  patterns.  Varying  the 
number  of  time  bins,  the  width  and  intensity  of  bars,  and  the  draw  events  which  are 
displayed  can  highlight  the  presence  of  a  significant  pattern  for  an  observer.  TOPOS 
also  makes  the  computer  an  ally  in  helping  the  user  to  interpret  the  displays.  The 
correlation  grid  iri  the  TOPOS  display  is  such  an  application  of  computer  processing 
to  bring  together  information  for  better  assimilation. 

Because  it  is  a  computer  program,  TOPOS  works  only  with  digital  inputs. 
Some  test  signals  and  the  available,  analog  EEG  data  were  interfaced  to  the  topo¬ 
scope  using  the  results  of  investigations  which  discovered  the  Brain  Atlas  digitization 
and  file  formats.  Sinusoidal  test  inputs  were  used  to  explore  the  tool's  indications  of 
frequency  and  phase  relationships  between  two  signals.  A  single  EEG  file  recorded 
during  a  severe  motion-sickness  episode  was  successfully  interfaced  to  TOPOS.  It  was 
then  studied  for  evidence  of  fundarm  al,  sub-harmonic,  and  harmonic  relationships 
with  an  external,  low-frequency  signal,  as  well  as  between  channels  in  the  record. 

As  the  test  and  analysis  phase  progressed  and  the  lesulting  TOPOS  displays 
were  collected  and  studied,  it  became  increasingly  evident  that  TOPOS  was  not  dis¬ 
playing  the  anticipated  brain  activity  pattern  in  a  form  recognizable  to  observers. 
This  raised  two  issues  of  significance.  The  first  issue  was  the  display's  inability  to  un¬ 
ambiguously  differentiate  between  any  two  significant  frequency  relationships.  The 
common  patterns  shared  by  fundamentally  related  signals  and  certain  harmonics, 
for  example,  reduced  the  usefulness  of  the  tool.  The  second  issue  was  the  lack  of 
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any  overall  increase  of  channel-reference  correlations  in  the  motion-siekness-affected 
EEG  record.  These  two  issues  begged  the  question  “Why?” 

In  attempting  to  answer  this  inquiry,  the  tool's  data-handling  integrity  was 
the  first  area  considered.  Soon  after  the  first  test  signals  were  displayed  on  TOPOS, 
quality  control  became  a  regular  concern.  Internal  computations  and  data  manipu¬ 
lations,  as  well  as  the  graphics  they  produced  were  thoroughly  cross-checked. 

The  next  area  of  investigation  was  TOPOS’  design.  An  initial  decision  to  adapt 
Walter  and  Shipton's  circular  format  to  a  line-based  display  began  the  evolution  of 
how  TOPOS  functioned.  Design  restrictions  came  and  went,  such  as  the  peak 
drawing  mode,  as  the  display  changed  to  increase  readability  and  usability.  This 
usability  came  in  both  the  form  of  user-control  and  the  ability  to  focus  on  frequencies 
below  4  Hz.  Central  to  the  final  display  was  the  use  of  reference  signal  axis-crossings 
for  frequency  derivation.  This  method  was  also  used  by  Shipton,  but  TOPOS  did 
not  offer  a  corresponding  display  of  channel  amplitudes  as  Shipton’s  display  did,  and 
the  opportunity  was  created  for  confusion  between  signals  of  differing  frequencies, 
but  similar  axis  crossings.  While  the  Waiter-Shipton  device  effectively  displayed 
the  amplitude  of  channel  activity  between  reference  signal  axis  ciossings,  TOPOS 
“strobes”  the  channel  signals  and  measures  their  amplitude  only  at  the  sample  where 
the  reference  signal  crosses  the  axis. 

The  third  area  of  consideration  was  the  tool's  application.  All  of  the  Walter 
and  Shipton  references  reviewed  during  this  research  addressed  the  utility  of  their 
Toposcope  in  the  analysis  of  evoked  potential  EEQs.  This  application  confined  their 
attentions  to  extremely  short  periods  of  time  following  each  stimulus  provided  to 
the  subject,  and  when  these  stimuli  were  repeated  regularly  and  rapidly,  Walter  and 
Shipton  could  observe  the  brain's  response  over  and  over  again.  In  this  research, 
however,  the  arrival  of  the  various  motion  stimuli  at  the  brain  could  not.  be  timed, 
much  less  controlled,  leaving  the  observer  to  look  for  “suspicious”  patterns  developing 
over  long  stretches  of  an  EEG  record.  This  led  to  the  use  of  Vogen’s  contour  brain 
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maps  as  the  focus  of  the  EEG  analysis.  His  maps  spanned  a  small  segment  of  EEG 
data  and  offered  a  position  from  which  to  begin  the  pattern  search. 

The  final  area  of  investigation  was  the  validity  of  the  brain  activity  pattern 
mapped  by  Vogen.  This  was  never  really  questioned,  and  especially  given  the  design 
and  performance  problems  of  TOPOS,  should  not  be.  A  goal  of  this  research  had 
been  to  study  motion-sickness-affectrd  EEG  data  with  the  hope  of  observing  the 
instantaneous  development  of  the  brain  activity  patterns  and  increasing  the  under¬ 
standing  of  the  phenomenon,  not  to  disprove  it. 

While  any  user  could  benefit  from  more  time  to  study  the  output  of  TOPOS, 
the  design  issues  outlined  above  are  considered  to  contribute  the  most  to  its  limited 
success  as  a  signal  analysis  tool.  Its  task  of  simplifying  the  confusing  relationships  in 
the  human  EEG  is  formidable,  but  achievable.  In  1953,  Grey  Walter  was  addressing 
the  complexity  of  the  brain’s  electrical  activity  and  he  admonished,  “...it  would 
be  foolish  to  expect  the  enchanted  loom  to  weave  a  fustian  fabric”  (45:292).  It 
would  also  be  short-sighted  to  assume  TOPOS  to  be  incapable  of  displaying  a  more 
significant  portion  of  that  weave. 
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Appendix  A.  Preprocessing  the  Data 


This  appendix  highlights  the  knowledge  gained  from  the  EEG  data  records 
used  as  ir  to  the  toposcope.  Fourteen  data  records  existed  in  analog  form  on 
tape  and  in  digital  form,  but  in  an  unknown  Brain  Atlas  format,  on  AFIT's  Bio¬ 
logic  computer.  The  goal  was  to  either  re-digitize  the  analog  tapes  or  disco%-er  the 
format  of  the  digital  records,  and  then  convert  their  contents  into  their  original 
voltage  values. 

A.l  Re- Digitizing  Options:  NeXT  and  CODAS 

Given  the  unidentified  format  of  the  Brain  Atlas  data  files  and  past  difficulties 
in  obtaining  technical  information  on  the  program,  the  f  :st  point  of  investigation  was 
the  possibility  of  re-digitizing  the  records  from  the  analog  tapes.  While  the  NeXT 
workstations  in  the  AFIT  signal  processing  lab  could  perform  A-D  conversions,  they 
could  only  convert  one  channel  at  a  time.  This  was  a  significant  problem  since  the 
EEG  records  were  14  time-synched  channels.  If  a  NeXT  was  used,  a  series  of  synch 
pulses  would  have  to  have  been  recorded  over  each  channel,  each  channel  individually 
digitized  into  a  separate  file,  and  then  the  files  merged  together  while  simultaneously 
recovering  the  channel  synchronization. 

Another  alternative  was  CODAS,1  a  signal  processing  package  available  on  a 
personal  computer  (PC)  in  the  motion  sickness  lab.  This  program  could  simulta¬ 
neously  digitize  three  to  five  channels  using  a  special  boaid  installed  in  the  CPU. 
However,  this  option  would  have  still  required  the  creation  of  multiple  files  and  the 
resynchronization  of  the  channels  when  recombined. 

It  was  very  important  that  all  14  EEG  channels  retain  tluir  time  synchroniza¬ 
tion  in  order  to  ensure  the  validity  of  any  brain  activity  pattern  displayed  by  the 

'A  software  pro.luc".  <jf  Data/;  Instruments,  Inc . 
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toposcope.  Therefore,  before  expending  large  amounts  of  time  and  effort  to  learn 
the  synchronization  techniques  and  gather  the  equipment  required  by  the  NeXT  and 
CODAS  options,  it  was  decided  to  examine  the  Brain  Atlas  data  files  in  search  of 
insight  into  their  storage  format. 

A.  2  Brain  Atlas 

The  task  of  decoding  the  method  which  Bio  logic  Systems  Corporation  used 
for  data  storage  was  made  much  easier  when  they  mistakenly  sent  a  copy  of  that 
format  to  AFIT.  Because  of  past  difficulties  in  obtaining  information  from  the  com¬ 
pany,  obtaining  a  copy  of  the  file  structure  had  previously  been  thought  highly 
unlikely.  However,  typographical  errors  were  discovered  in  the  Bio-lc\  ■•ocument. 
The  following  section  will  discuss  some  of  the  specifics  of  the  Bio-logic  LEG  data  files 
resulting  from  the  study  of  Brain  Atlas  documentation  and  experimentation.  The 
procedures  referenced  to  the  appendices  are  based  on  the  Brain  Atlas  Users  Guide 
(4)  but  were  developed  during  the  experimentation  phase  to  improve  on  the  user 
interface  offered  by  the  user  guide. 

A. 2.1  Bio-logic  EEG  Data  File  Structure  The  Brain  Atlas  software  located 
on  the  Bio-logic  PC  in  the  Motion  Sickness  Lab  is  capable  of  simultaneously  dig¬ 
itizing  H  channels  of  analog  EEG  data.  Howear,  the  binary  files  it  creates  in¬ 
clude  six  additional  cl  ■  which  it  internally  interpolates  prioi  to  data  storage 
from  i;  use  channels  actually  recorded.  The  names  of  these  files  follow  the  format 
"Exxxxxxx.DAT  to  differentiate  them  from  evoked  potential  (EP)  data  files,  which 
are  prefixed  with  the  letter  “F."  Due  to  AFIT’s  nominal  collection  rate  of  G4  samples 
per  second  and  average  recording  time  of  20  minutes,  the  EEG  data  files  range  in 
si/e  from  uOO  to  800  kilobytes.  These  are  biliary  files  using  an  MS-DOS  format.  The 
first  1920  bytes2  are  header  information  in  a  pattern  common  to  both  EEG  and  EP 


files.  Therefore,  much  of  the  header  content  is  invalid  for  EEG  applications,  lhose 
values  which  are  important  appear  in  Equation  2  on  page  73  and  are  discussed  in 
that  section.  Also,  the  header  size  of  1920  bytes  is  different  than  that  shown  in  the 
Bio-logic  documentation.  The  documentation  errors  found  during  the  research  are 
listed  below. 

•  Header  length  —  Header  is  15  blocks3  long,  not  13  as  printed  (3:3). 

•  EEG  Data  Points  —  Begin  with  byie  1920. 4  Stored  as  1-byte  integer  numbers 
indicating  analog-to-digital  (A-D)  counts,  not  as  “sums  saved  as  three-byte 
integers”  (3:10). 

•  Voltage  Conversion  Equation  —  Produces  values  with  offsets  by  channel  of  0 
to  4  pvolts  from  those  used  by  the  Brain  Atlas. 

•  Header  location  of  SENS  and  (JALUV  —  Both  are  1-byte  (vs.  2-byte)  integers 
beginning  at  the  location  shown  in  the  Biologic  format  documentation. 

The  location  of  key  EEG  record  parameters  within  the  Bio-logic  data  files  are 
provided  in  Table  1. 


Table  1.  Key  to  Brain  Atlas  Data  File  Header.  All  values  are  in  hexadecimal. 

Notes:  1)  Invert  for  approximate  sample  rate  (either  64  or  128  Hz);  2)  See 
Table  4  for  location  of  each  channel. 


FIELD 

BYTES 

DESCRIPTION 

Sample  Period 

416 

in  msecs  (Note  l) 

Calibration  Values  (C.V.) 

425-438 

offset  in  A-D  counts  from  +80h  (Note  2) 

D.C.  Offsets  (D  C.) 

446-459 

offset  in  A-D  counts  from  +80h  (Note  2) 

Sensitivity  (SENS) 

467 

std  EEG  parameter  in  p  V/mm 

Calibration  Voltage  (CALUV) 

469 

in  p  V  at  std  SENS 

31  block  —  128  bytes. 

'’Header  contains  1920  bytes  numbered  0  to  1919. 


A. 2. 1.1  E EG  Data  Points.  Within  the  body  of  the  data  file  are  a  series 
of  1-byte,  integer  A-D  counts  representing  instantaneous  voltages  recorded  over  the 
EEG  channels.  They  range  from  Oh5  to  +0ffh.  However,  080h  is  considered  zero. 
Thus,  A-D  data  values  less  than  080h  represent  negative  numbers  and  the  actual 
number  of  collected  counts  for  any  data  point  is  found  by  subtracting  the  artificial 
zero  level  080h  from  the  binary  value  found  in  the  data  file.  This  offset  is  included 
in  the  voltage  calculation  shown  in  Equation  2. 

Prior  to  display  in  a  Brain  Atlas  plot  or  topographic  map,  Bio-logic  states  that 
these  A-D  counts  are  converted  into  voltages  using  the  following  equation  (3:10): 


Voltage^,. 


[{ADV- 80h)  -  (DC-  80h)j  x  (2  x 


CALUV  K  SENS 

-ce-)xTo- 


(2) 


where:  ADV  =  number  of  A-D  counts 

80h  =  the  assigned  zero- level 
DC  =  DC  offset  for  each  channel  (from  header) 

CALUV  =  size  of  calibration  signal  in  microvolts  (from  header) 
CV  —  CALUV  in  A-D  counts 


The  results  of  this  equation  should  !>«•  the  values  displayed  anytime  the  "Volt¬ 
age"  option  is  selected  within  t'.e  Cursoi  Analysis  funct  ior  option  of  the  Brain  Atlas 
(see  Appendix  F).  or  when  the  program  is  tasked  to  provide  an  ASCII  output  of 
a  subset  of  the  collected  data  (see  Appendix  E).  However,  this  equation  produces 
values  which  are  0  to  4  pvolts  less  than  t;ie  vdues  calculated  md  displayed  by  the 

6 The  format  Oxh  drnotes  hexadecimal. 
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Table  2.  Offsets  of  Bio-logic  Voltage  Conversion  Equation.  Differences  are  between 
calculated  values  and  values  produced  by  the  Brain  Atlas  and  are  in  /ivolts. 


|Ch 

Offset 

Ch 

Offset 

sim 

-2 

8 

-4" 

2 

-2 

9 

-1.98 

3 

+0 

10 

-4 

in 

-2 

11 

-1.98 

n 

_2 

12 

-2 

6 

-4 

13 

-2 

7 

mm 

14 

-2 

Brain  Atlas.  These  offsets  for  the  recorded  channels  are  constant  across  time  but 
vary  by  channel  as  shown  in  Table  2. 

The  offsets  for  interpolated  channels  were  not  tracked.  The  reason  for  these 
anomalies  was  not  found,  but  a  documentation  error  is  possible,  given  the  incidence 
of  errors  as  previously  noted. 

A. 2. 1.2  Uncovering  the  EEG  Data  Format  Using  Norton  Utilities,  the 
MS-DOS-format  Bio-logic  data  files  were  output  in  their  hexadecimal  notation  and 
used  to  verify  the  location  and  validity  of  key  header  data  as  shown  on  the  Bio-logic 
format  sheet.  The  documentation  was  correct  and  showed  the  21  1-byte  values6  of 
CV  contained  in  bytes  425-445  of  the  file  header.  Values  of  DC  offsets  by  channel 
were  located  in  bytes  446-466.  SENS  and  CALUV  are  actually  l-byte  integers  (vs. 
2-byte,  as  published)  beginning  at  bytes  437  and  469,  respectively. 

There  is  an  error  in  Bio-logic’s  description  of  how  the  EEG  data  points  were 
stored  in  the  body  of  the  data  file  (1).  The  points  are  not  recorded  as  3-byte, 
integer  sums,  but  as  lbyte  integers  representing  the  number  of  A-D  counts  from 

^Although  the  Drain  Allas  only  stores  and  displays  20  channels,  Ihe  format  for  header  parameters 
contains  locations  for  21  values. 


74 


each  channel  and  for  each  sample  obtained.  The  number  of  samples  depends  on  the 
update  rate  selected  and  the  length  of  the  recording. 

Port'  ns  of  the  data  files  were  output  using  the  ASCII  Conversion  option  in  the 
Bank  Math  function  of  the  Brain  Atlas  (1).  This  tool  only  converts  256  data  points, 
but  using  the  procedure  detailed  in  Appendix  E,  ASCII  values  of  voltages  could  be 
produced.  Vet  the  question  of  the  location  of  each  channel’s  value  within  the  block 
of  20  or  21  data  points  (see  Figure  27)  for  a  single  sample  was  still  unanswered.  For 
example,  was  Channel  1  the  first  data  point  or  located  elsewhere? 


H 

E 


.  Cluot«r 
.FPFFFFPF 
|r[’p  PFXPF 
.  FPFFFFFF 
FPFFFFPF 

ffffffff 

. FFFFFFFP 
. QCFFCOFF 
i  .OOFFCOFF 
;~n;  .42482331 
.43403035 
•ii  .43482039 
.43482031 
.49465445 
* .46372000 
•  .54332000 
. Ji  .54352000 


19  Socto 

f  FFFFFFF 

FPFFFFFF 

FFFFFFFP 

F  FFFFFFF 

FPFFt  FFK 

FFFFFFFP 

0FPF10FF 

OOFFOOFr 

20204318 

20204348 

20204348 

13204348 

5250494E 

52654633 

52654333 

52655031 


r  46 

FFFFFFFP 

FFFPPFFP 

PFFFFFFP 

FFPFFFFF 

FFFFFPFF 

02FF03FP 

11FF12PF 

oorroorr 

20322C2Q 
20363020 
20313030 
2031 1420 
54455250 
20005265 
20005265 
20005265 


FFPFFFFF 
F  PFPFFFP 
PFFFFPFP 
FFFFPFFF 
FFFFFPFF 
05FF06FF 
O0FF01FF 
OOFFOOFr 
43482033 
43482037 
43482031 
494B5445 
494E5445 
16342000 
43312000 
50343000 


■  8E858B83 
.80866084 

.858A847E 

•S3727F8U 


BFB88.489  S18C88BA  91946F6B 


8A858 181 
88SB8987 
837A8890 
9001 9A97 


86868483 

90899680 

967P9697 

3D847S89 


39859380 
977B9797 
8C817EB9 
898A8C8  8 


FFPFFFFF 
FPFFFFPF 
FFFFFPFF 
FPFFFFPF 
FPFFFFPF 
07FPC8FP 
04  FFC9FF 

oorroorr 

20204348 

20204348 

3120434B 

5250494E 

52504946 

52654638 

52655434 

52655436 

jaajB?cad 


9F83 A897 
887A7D7D 
88898A87 
87837F7D 


FFPFFFFF 

FFFFFFFP 

FPFFFFFF 

FFPFFFFF 

FPFFFFFF 

OAFFOBFF 

OEPF13PF 

OOFFOOFF 

20342020 

20382020 

20313220 

54455250 

54455250 

20005365 

20005265 

20005265 

9F8BAB97 

8C8285B2 

82807E81 

8985807C 

85737E0U 


Figure  27.  Hexadecimal  Representation  of  a  Brain  Atlas  EEG  Data  File  Segment. 

The  boxed  numbers  are  one  sample  of  20  channels  of  EEG  data  (6 
channels  are  interpolated). 


This  was  especially  difficult  since  the  Brain  Atlas  signal  junction  box  labeled 
the  channel  inputs  numerically,  but  the  Brain  Atlas  screen  used  alphanumeric  des¬ 
ignations.  In  addition,  the  Brain  Atlas  assumed  monopolar  recordings,  which  use  a 
common  reference  for  all  electrodes,  while  AFIT  held  made  bipolar  recordings  which 
used  neighboring  electrodes  as  recording  references.  This  difference  further  confused 
the  channel  assignment  search. 


Table  3.  Correlation  Between  Numeric,  AFIT,  and  Bio-logic  Channel  Labels. 


Channel 

AFIT 

Brain  Atlas 

I 

F7 

2 

K||pll 

F3 

3 

F4 

4 

C3-T5 

F8 

5 

T5-P3 

T3 

6 

P3-01 

C3 

7 

01-C2 

C4 

8 

02-P4 

T4 

9 

P4-T6 

T5 

10 

T6-C4 

P3 

11 

C4-T4 

P4 

12 

T4-F8 

T6 

13 

F8-F4 

01 

14 

F4-F3 

02 

To  determine  the  channel  locations  within  the  data  block,  the  Brain  Atlas 
was  calibrated  using  the  procedure  in  Appendix  C  and  a  200  mV  square  wave  was 
collected  following  the  procedure  in  Appendix  D.  EEG  data  was  played  over  all  but 
the  one  channel  containing  the  square  wave  and  14  separate  files  were  recorded,  each 
with  the  square  wave  located  on  a  different  channel.  The  initial  values  of  each  file 
were  then  converted  to  ASCII  and  the  values  correlated  with  those  displayed  using 
the  “Voltage”  option  in  the  Brain  Atlas.  This  was  possible  because  the  square  wave 
maintained  a  constant  value  across  enough  samples  to  clearly  indicate  the  match 
between  input  channel  and  Brain  Atlas  electrode  position.  The  resulting  assignments 
of  input  channel  number  and  AFIT  and  Brain  Atlas  electrode  locations  are  shown 
in  Table  3. 

Determining  these  assignments  made  it  possible  to  define  the  format  of  data 
storage  and  locate  all  the  necessary  values  to  input  into  Equation  2  in  order  to 
convert  the  A-D  counts  into  voltages.  The  channel  assignment  within  the  DC  offset 
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Table  4. 


Channel  Locations  in  a  Brain  Atlas  Header  Block  and  ASCII  File  (***  in¬ 
dicates  locations  of  interpolated  values). 


*** 

*** 

*  *  * 

Ch  01 

Ch  02 

*** 

Ch  03 

Ch  04 

Ch  05 

Ch  06 

*** 

Ch  07 

Ch  08 

Ch  09 

Ch  10 

*** 

Ch  11 

Ch  12 

Ch  13 

*  *  * 

Ch  14 

Table  5.  Channel  Locations  in  a  Single  Sample  From  a  Brain  Atlas  Binary  File 
(***  indicates  locations  of  interpolated  values). 


■mi 

*** 

Ch  02 

Ch  03 

Ch  04 

Ch  05 

Ch  06 

HB 

Ch  07 

Ch  08 

fauna 

Ch  10 

Ch  11 

Ch  12 

Ch  13 

Ch  14 

*** 

and  CV  blocks  in  the  header  is  the  same  as  in  the  Brain- Atlas-generated  ASCII  file 
and  is  shown  in  Table  4.  The  channel  assignments  for  the  EEG  data  samples  are 
contained  in  Table  5.  Using  these  data  locations,  the  Equation  2  could  be  used  to 
calculate  voltages  and  determine  the  offsets  in  Table  2. 

A. 3  Reading  EEG  Data  Files 

This  knowledge  of  data  locations  and  conversion  calculations  provided  all  the 
required  input  to  process  the  EEG  data  within  a  dio-logic  file.  The  process  involved 
transferring  the  data  to  the  computer  to  be  used  during  the  programming,  writing 
the  code,  and  testing  and  validating  the  program  and  its  output. 

A -3.1  Transferring  the  Data.  It  was  anticipated  that  a  format  conversion 
might  be  required  before  the  MS-DOS  binary  data  could  be  correctly  read  by  the 
Unix-based  Silicon  Graphics  Iris/4 D  workstation  at  AFIT.  However,  no  conversion 
was  necessary.  The  EEG  file  was  transferred  from  the  Bio-logic  PC  to  h  network 
PC  in  the  AFIT  Motion  Sickness  Lab  using  a  data  diskette.  Frcm  thore  the  Silicon 
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Graphics  workstation  was  accessed  using  the  File  Transfer  Protocol  (FTP)  applica¬ 
tion  and  the  file  transferred  in  binary  mode. 

A  3.2  The  Data  Preprocessing  Program.  The  goal  of  this  programming  effort 
was  to  calculate  the  EEG  channel  voltages  from  the  data  contained  in  the  file.  The 
first  step  was  to  produce  a  table  of  the  values  from  the  file  header  which  are  required 
in  the  voltage  conversion  equation.  These  values,  along  with  the  channel  samples 
from  the  body  of  the  file,  could  then  be  used  to  calculate  the  voltages  of  each  channel 
for  each  sample  in  the  input  file.  Tnese  tables  were  created  early  in  the  programming 
effort  to  use  in  ensuring  that  the  data  was  being  handled  correctly.  A  portion  of 
the  results  are  in  Appendix  B,  but  the  code  was  not  preserved  in  the  form  used 
to  produce  the  tables.  Instead,  it  was  converted  into  functions  which  output  their 
results  for  graphical  interpretation  on  the  display. 

Since  each  data  point  was  stored  in  a  l-byte  location  in  memory,  they  were 
read  from  the  Brain  Atlas  file  as  C  char  (character)  data  types  since  they  also 
require  one  byte  each  (40:A-4).  The  data  points  were  then  converted  to  integers 
prior  to  use  in  calculations.  Required  values  from  the  file  header  (e.g.,  calibration 
voltages)  were  read  into  arrays  and  then  used  to  convert  each  EEG  sample  from 
A-D  counts  into  microvolts  prior  to  writing  it  to  the  appropriate  output  file.  No 
correction  was  made  for  the  offset  between  the  calculated  value  and  the  Brain  Atlas 
value  (see  Table  2).  Fourteen  separate  output  files  were  created,  each  holding  data 
from  one  EEG  channel.  Data  from  each  channel  were  also  analyzed  to  determine  the 
approximate  distribution,  with  a  maximum,  minimum,  and  standard  deviation  for 
each  channel  written  to  a  single  header  file.  Within  the  toj  oscope  code,  the  function 
prep-data()  performs  data  preprocessing. 
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A.  4  Summary 

The  conversion  between  analog  EEG  signals,  binary  Brain  Atlas  data  files,  and 
C-program  output  quickly  became  routine  and  dependable.  This  was  possible  due  to 
the  discovery  of  the  Brain  Atlas  file  format  and  the  flexibility  of  the  C  programming 
language. 
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Appendix  B.  Data  Obtained  from  Brain  Atlas  Data  File 


CV  Values  by  Channel: 


BYTE 

CH 

INTEGER 

HEX 

428 

1 

200 

c8 

429 

2 

200 

c8 

431 

3 

200 

c8 

432 

4 

200 

c8 

433 

5 

200 

c8 

434 

6 

200 

c  8 

436 

7 

198 

c6 

437 

8 

200 

c8 

438 

9 

202 

ca 

439 

10 

200 

c8 

441 

11 

202 
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Appendix  C.  Brain  Atlas  Data  Collection  Calibration 


PURPOSE:  To  prepare  the  Brain  Atlas  (BA)  program  for  data  collection.  Not 
required  if  you  have  not  edited  the  program  since  last  calibration. 

STEP  1:  At  the  C:>  prompt,  type  BA  ai.d  press  [ENTER], 

STEP  2:  Select  the  CPU  turbo  mode. 

ACTION:  Press  [ALT]-[CNTL]-[-]( “minus"  on  keypad). 

STATUS:  Lights  the  “Fast  Clock"  LED  on  front  of  CPU. 

STEP  3:  Select  “EEG  Collect"  on  Main  menu. 

STEP  4  '  Select  “Calib”  on  Collection  menu. 

STEP  5:  Create  400 mV  1-Hz  square  wave  using  a  function  generator  and  o-scope. 
STEP  6:  Set  GEN  MODE  on  function  generator  to  TRIG. 

STATUS:  Oscilloscope  shows  flat  trace. 


CAUTION:  Make  A-D  port  connection  carefully  to  avoid  shorting  the  A-D 
board  in  the  CPU. 

STEP  7:  Route  calibration  signal  to  A-D  port  on  CPU  using  25-pin  connector. 

STEP  8:  Select  "Start  Calibration"  and  wait  for  the  beep. 

STEP  9:  Approximated  0.5  secs  after  the  beep,  turn  GEN  MODE  switch  on  function 
generator  to  CONT. 

STEP  10:  Count  4  wavelengths  of  square  wave  and  then  return  GEN  MODE  to  TRIG. 

If  BA  indicates  "Calibration  out  of  limits,”  do  Steps  11- 13. 

STEP  11:  Press  [ESC]  and  note  calibration  wave  forms  collected. 

STEP  12:  Adjust  timing  of  GEN  MODE  switch  positioning  on  subsequent  attempts  to  fill 
the  time  window  displayed. 

STEP  13:  Repeat  Steps  7-10  as  needed. 
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If  BA  indicates  successful  calibration,  continue. 


STEP  14-'  Remove  25-pin  connector  from  A-D  port  on  CPU. 


CAUTION:  Exiting  the  DA  program  deletes  this  calibration! 
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Appendix  D.  Brain  Atlas  Data  Collection 


PURPOSE:  To  digitize  analog  EEG  data  or  test  stgnals  and  create  DOS  binary 
files  (“Exxxxxxx.dat")  in  the  Drain  Atlas  (BA)  format. 

STEP  1:  Ensure  the  BA  program  is  calibrated  (see  Calibration  Procedure). 

STEP  2:  Select  the  following  in  sequence: 

"EEG  Collect”  on  Main  menu. 

"Hardware”  on  C’olection  menu. 

"Sample  Rate”  on  pull-down  menu. 

"64  Hz”  on  option  menu. 

STEP  3:  Turn  on  Disking  function  and  open  new  data  file. 

ACTION:  Select  “File”  on  Collection  menu. 

Select  "Disking”  on  pull-down  menu. 

Enter  7-digit  file  name  (without  "E"  prefix 
or  “.dat"  extension)  and  select  "OK.” 

STATUS:  DSK  ON  flag  on  right  side  of  screen. 

STEP  4'  Connect  Bio-logic  BNC  channel  input  box  to  A-D  port  on  CPU. 


NOTE:  Data  collection  does  not  begin  until  approximately  2  seconds  after 
“Start"  is  selected. 

STEP  5:  Initiate  input  signal  (e.g.,  turn  on  recorder). 

STEP  6:  Select  "Start”  on  Collection  menu  (2  second  delay). 

STEP  7:  To  pause  data  collection,  press  [ESC]  or  right  mouse  button  once. 

STEP  8:  To  stop  collection,  press  [ESC]  or  right  mouse  button  twice. 

STEP  9:  Disconnect  input  cable  from  A-D  port  on  CPU. 
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Appendix  E.  Brain  Atlas  ASCII  Conversion 


PURPOSE:  Converts  the  first  256  samples  of  digitized  Brain  Atlas  (BA)  data 
into  ASCII  format.  Data  must  first  be  stored  in  a  BA  "bank.  ” 

NOTE:  The  Cursor  Analysis  bank  data  option  is  not  included  here  (see  BA 
manual,  page  EEG-191). 

STEP  1:  At  the  C:>  prompt,  type  BA  and  press  [ENTER]/ 

STEP  2:  Select  the  CPU  turbo  mode. 

ACTION:  Press  [ALT]  [CNTLj-[-j  ( “minus''  on  keypad). 

STATUS:  Lights  the  “Fast  Clock’’  LED  on  front  of  CPU. 

STEP  3:  Select  the  following  in  sequence: 

“EEG  Analysis”  on  Main  menu. 

“File”  on  Analysis  menu. 

“Disking”  on  pull-down  menu. 

STEP  4:  Select  desired  drive  and  data  file. 

STEP  5:  Select  “Banks”  on  Analysis  menu. 


NOTE:  Attempting  to  allocate  more  than  one  bank  can  result  in  insufficient 
memory  for  later  operations. 


STEP  6:  Select  “Bank  Data'  on  pulldown  menu. 

STEP  1:  Select  “1”  on  top-line  option  menu  to  bank  data. 
STATUS:  None  provided. 

STEP  8:  Press  [ESC]  or  right  mouse  button  to  exit  to  Main  menu. 
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STEP  9:  Select  the  following  in  sequence. 

“Map  Analysis’’  on  Main  menu. 

“Process”  on  Analysis  menu. 

“Dank  Math”  on  pull-down  menu. 

STEP  10:  Tab  to  the  following  fields,  make  entries,  and  press  [ENTER]: 

Source  bank:  1 

Destination  bank:  A 

STEP  11:  Press  [SHIFT]— [F3]  to  transfer  data  into  a  Bank  Math  letter  bank. 

STATUS'  Lines  at  top  of  screen  indicate  “A"  bank  is  allocated. 

STEP  12:  Back  tab  to  the  following  fields,  make  entries  and  press  [ENTER]: 

Source  bank:  A 

NOTE:  Failing  to  deallocate  the  letter  bank  will  result  in  insufficient  memory  for 
later  operations. 

STEP  lp.  Press  [F8]  to  deallocate  the  letter  bank. 

STEP  15:  Press  [ESC]  or  right  mouse  button  to  exit  Bank  Math. 

To  view  the  ASCII  file  contents,  continue. 

STEP  16:  Select  “File"  on  the  Analysis  menu. 

STEP  17:  Select  “Temp  Exit”  on  the  pull-down  menu. 

STEP  18:  Type  type  txxxxxxx.asc  L  more  and  press  [ENTER]. 

STEP  19:  Press  [ENTER]  to  view  subsequent  pages. 

STEP  20:  At  C:>  prompt,  type  exit  and  press  [ENTER]  to  return  to  BA. 
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Appendix  F.  Brain  Atlas  Cursor  Analysis  of  Data 


PURPOSE:  To  display  instantaneous  voltages  of  all  channels  in  a  Brain  Atlas 
(BA)  data  file.  Also  provides  an  interpolated  brain  map  based  on  those  voltages. 

STEP  1:  At  the  C:>  prompt,  type  BA  and  press  [ENTER]. 

STEP  2:  Select  the  CPU  turbo  mode. 

ACTION:  Press  [ALT]-[CNTL]-[-]  (“minus”  on  keypad). 

STATUS:  Lights  the  “Fast  Clock"  LED  on  front  of  CPU. 

STEP  3:  Select  the  following  in  sequence: 

"EEG  Analysis”  on  Main  menu. 

“File”  on  Analysis  menu. 

“Disking"  on  pull-down  menu. 

STEP  4:  Select  desired  drive  and  data  file. 

STEP  5:  Select  the  following  in  sequence: 

“Display”  on  Analysis  menu. 

“Speed”  on  pull-down  menu. 

"Double”  on  option  menu  to  display  all  data  pts. 

“Montage"  on  Analysis  menu. 

"Define”  on  pull-down  menu. 

“Select  Group”  on  User  Montage  menu. 

“I”  on  pull-down  menu. 

“Select  Montage"  on  User  Montage  menu. 

“0”  on  pull-down  menu. 

STEP  6:  Backout  to  Analysis  menu  using  [ESC]  or  right  mouse  button. 

STEP  7:  Select  “Cursor”  on  Analysis  menu. 

To  vary  the  scale  of  displayed  plots,  cont inue. 

STEP  8:  Select  '  uisplay  on  the  Cursor  menu. 

STEP  select  "Raise”  or  “Lower”  on  the  pull-down  menu,  as  needed. 


NOTE:  Start  time  of  displayed  data  is  located  at  the  bottom-left  comer  of  the 
screen. 

To  move  to  earlier/later  portion  of  data  by  entering  start  time,  continue. 

STEP  10:  Select  “Tinier”  on  Cursor  menu. 

STEP  11:  Select  “Set  to  hh:mm:ss”  on  Set  Timer  menu. 

STEP  12:  Enter  desired  start  time  and  select  “OK." 

NOTE:  Fine  adjustment  of  cursor  position  is  via  the  straight-arrow  icons  at 
bottom-nght  of  the  screen.  Time  of  selected  data  point  is  displayed  at  bottom-left 
comer. 

To  move  to  earlier/later  portions  of  data  one  screen  at  a  time,  continue. 

STEP  13:  Select  Z-shaped  arrow  icons  a*  bottom-right  of  the  screen  for  the  desired  direction. 
To  view  a  brain  inap  of  interpolated  voltages,  continue. 

STEP  14-'  Position  cursor  by  clicking  the  right  mouse  button  at  the 
desired  point  in  the  displayed  plot. 

To  view  instantaneous  voltages  in  each  data  channel,  continue. 

STEP  15:  Position  cursor  by  clicking  the  right  mouse  button  at  the 
desired  point  in  the  displayed  plot. 

STEP  16:  Select  “Voltage”  on  Cursor  menu. 

NOTE:  Screen  cannot  be  printed  while  “ Voltage  "  display  is  present. 

To  print  the  screen,  continue. 

STEP  17.  Select  the  printer  icon  at  the  bottom-right  corner  of  ine  scrcu: 
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Appendix  G.  TOPOS  User’s  and  Programmer’s  Guide 


Gl  Introduction 

This  guide  is  designed  for  two  classes  of  readers:  the  user  and  the  student 
or  researcher  who  wishes  to  become  familiar  with  the  code’s  design  in  order  to 
make  modifications.  Design  notes  will  be  mostly  high  level  and  intended  to  be 
supplemented  by  the  body  of  the  thesis  report  it  accompanies,  as  well  as  the  TOPOS 
code  included  in  Appendix  I.  While  the  program  behind  the  tool  was  written  with 
great  pride  and  care,  it  is  not  the  product  of  an  experienced  programmer,  so  leniency 
and  empathy  is  requested  from  all  future  reviewers. 

G.2  Why  TOPOS? 

This  tool  belongs  in  a  class  of  spatial/temporal  signal  analyzers  called  topo- 
scopes  and  is  designed  for  the  study  of  human  brain  electrical  activity.  The  first 
tGposcope  used  for  similar  purposes  was  conceived  and  constructed  by  Grey  Walter 
and  Harold  Shipton  in  the  late  1940s.  TOPOS  is  an  application  of  the  principles  of 
that  original  device.  Prior  to  the  advent  of  high-speed  computers  and  digital  process¬ 
ing,  Walter  and  Shipton 's  Toposcope  offered  researchers  a  picture  of  brain  activity 
based  on  frequency  commonalities  in  spatially  separate  locations  of  the  brain.  Those 
commonalities  led  to  their  early  reports  of  functional  links  between  neuron  groupings 
in  the  brain  as  they  worked  together  to  process  a  visual  stimulus. 

Almost  45  years  later,  AFIT  redeveloped  this  tool  on  a  computer  workstation 
to  study  brain  activity  patterns  resulting  not  from  a  visual  stimulus,  but  from  the 
stimuli  from  multiple  sources  accompanying  motion  sickness.  The  frequency  band 
below  4  Hz  was  the  region  of  interest,  and  so  several  of  the  options  are  specifically 
designed  to  focus  the  tool  on  those  frequencies. 
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As  a  toposeope.  TOPOS  offers  an  instantaneous  viewpoint  of  frequency  rela¬ 
tionships  between  separate  EEG  electrode  locations.  While  Fast  Fourier  Transforms 
(FFTs)  and  contour  maps  arc  popular  methods  of  processing  and  presenting  the 
same  information,  a  toposeope  does  so  without  the  time  averaging  required  by  FFTs 
or  the  interpolation  of  contour  maps. 

G.3  Using  TOPOS 

G.3.1  Prerequisites.  The  program  requires  the  following  files  to  run: 
topos 

sinetest  ,dat  .halfhz 
sinetcst.dat.  lhz 
sinetest.  dat.l_5hz 
sinetest. dat. 2hz 

an  input  data  file  in  either  Brain  Atlas  format 
or 

the  14  "datafile. chXX"  files  with  accompanying  “ datafile.h"  file. 

Attempting  to  start  the  tool  without  all  of  these  files  will  result  in  an  error  message 
to  the  console. 

G.3. 2  Starting  the  Tool.  TOPOS  is  started  by  typing  its  name,  followed  by 
the  data  file  name  as  shown: 


>  topos  inputfile 

If  the  input  file  is  a  raw  Brain  Atlas  data  file,  its  full  name  should  be  used.  If  the 
TOPOS  Prepare  Data  option  has  already  been  run  on  the  file,  the  name  of  the  raw 
data  file  should  still  be  used,  since  TOPOS  will  append  the  necessary  file  extensions 
before  accessing  the  files. 
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G.3.3  Help  Screen.  Upon  successful  startup,  a  help/instruction  screen  is  dis¬ 
played  which  provides  the  program’s  default  settings  of  display  options  and  overviews 
the  following  mouse  button  functions. 

Left  Mouse:  Halts  display 

Middle  Morse:  Continuous  Drawing  -  Press  pauses,  release  restarts 
Step  Drawing  -  Press/release  both  draw  next  step 
Batch  Drawing  -  Press  draw  next  batch  (redraws  all  bins) 

Right  Mouse:  Main  menu 

The  meaning  of  each  type  of  drawing  mode  is  explained  below. 

G.3.4  TOPOS  Control.  User  menus  greatly  simplify  the  task  of  viewing  a 
data  record  repeatedly  while  varying  display  parameters.  The  one  control  option 
not  located  on  a  menu  is  Display  Pause.  To  momentarily  pause  the  display  when 
Step-Draw  and  Batch-Draw  are  not  active,  the  user  should  depress  the  middle  mouse 
button  and  hold  it.  Releasing  this  button  restarts  the  display  at  the  point  it.  was 
paused. 

All  of  the  following  menu  options  are  available  on  the  main  menu.  Many  of 
the  descriptions  below  also  include  comments  on  usefulness. 

•  Prepare  Data  —  Selected  only  the  first  time  that  a  new  EEG  data  record  is  to 
be  displayed.  Data  preprocessing,  and  the  delay  it  requires,  can  be  bypassed  for 
that  data  record  from  then  on.  Subsequent  sessions  with  a  previously  processed 
data  file  uses  the  separate  channel  output  files  as  input  streams. 

t  Refercnct  Channel  —  Permits  the  online  selection  of  any  of  the  14  channels 
or  one  of  four  sinusoids  as  the  reference  signal.  Choosing  a  channel  will  cause 
TOPOS  to  display  relationships  of  all  channels  to  the  instantaneous  frequency 
(constantly  vat'ying)  of  that  reference  channel.  Selecting  a  sinusoid  as  a  refer¬ 
ence  results  in  a  display  of  all  channels  as  they  relate  to  a  fixed  frequency. 
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•  Bar  Width  —  Permits  selection  of  a  narrow  to  broad  bar  pattern  in  the  display. 
User  should  vary  to  their  taste.  Droader  bars  are  helpful  in  simpler  patterns 
and  narrower  bars  in  complex  displays. 

•  Bar  Intensity  —  Offers  three  options:  scale  to  reference  signal  value,  clamp  to 
maximum  value,  or  threshold  at  one  of  three  reference  signal  values.  Amplitude 
information  provided  bv  sealed  option  is  of  limited  value  and  can  be  confusing 
when  trying  to  ascertain  a  pattern  Thresholding  can  be  helpful  in  simplifying 
.  ittern  identification. 

•  Display  Mode  —  Offers  various  combinations  of  single-bar/time-shift  displays 
and  continuous/peak/axis-crossing  draw  event  triggers.  Single-bar  option  is 
entertaining,  but  rarely  useful.  Peak  drawing  is  fine  with  constant-amplitude 
test  signals,  but  complicates  the  pattern  for  EEG  signals  (peaks  of  which  vary 
widely  in  amplitude).  Axis-crossing  drawing  is  best  choice. 

•  Draw  Interval  —  Permits  the  modification  of  draw  event  intervals  by  skipping 
every  nth  event.  Useful  in  getting  down  to  lower  frequencies  of  both  smoothed 
and  unsmoothed  EEG  data.  Watch  "Reference  Signal  Frequency"  status  line 
on  screen  to  determine  best  selection. 

•  Bin  Number  —  Sets  the  number  of  bins  used  for  time-shifted  display.  A  low- 
bin  number  is  the  kev  to  displaying  recognizable  patterns  in  complex  signals. 
Helps  the  user  to  look  at  the  trees  instead  of  the  forest.  Also  note  that  the 
program  is  currently  set  to  slow  the  display  as  the  bin  number  is  decreased. 
This  permits  the  user  to  better  view  the  displays  presented. 

•  Start  Display  -  Allows  the  program  to  be  started  or  restarted  at  the  beginning 
of  the  data,  at  a  specified  time  mark,  or  at  the  point  where  the  display  was  last 
paused  during  that  session.  This  option  is  require.d  anytime  the  screen  shows 
‘'DISPLAY  HALTED"  over  the  status  box. 
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•  Test  —  Controls  the  tost  modes  of  single-buffer  drawing,  step-drawing  (draw 
one  time-shift  bin  and  wait),  and  batch-drawing  (draw  one  set  of  bins  and  wait). 
Single-buffer  drawing  is  the  way  around  the  "flicker"  plaguing  the  TOPOS 
display.  It  is  a  "must"  when  using  the  lower  bin  numbers  for  time-shift  drawing. 
These  displays  update  much  more  slowly  and  double-buffer  drawing  is  not 
suitable.  Step  drawing  is  useful  when  the  user  wants  to  control  the  display  of 
each  draw-event  with  the  middle  mouse  button.  Batch  drawing  prevents  the 
overwrite  of  any  bins  prior  to  the  user's  signal  via  the  middle  mouse  button. 
This  is  especially  helpful  when  capturing  images  from  the  screen. 

•  Help  —  Brings  up  the  holp/instruction  screen. 

•  Swap  Duffers  —  Continuously  swaps  the  two  display  buffers  to  enable  the 
observer  to  see  the  points  drawn  on  both  planes  whenever  the  display  is  paused. 

G.3.5  On-Screen  Information  User  information  appears  in  three  different 
locations  on  the  screen. 

•  Status  Display  —  In  the  lower-right  corner,  shows  the  data  file  being  displayed, 
the  sampling  rate  it  was  recorded  at,  the  elapsed  time1  since  the  start  of  the 
data  display,  and  the  instantaneous  reference  signal  frequency,  as  calculated 
from  the  number  of  samples  occurring  between  draw  events.  Other  messages, 
such  as  “Searching,"  “Display  Paused.”  and  “Display  Halted"  appear  in  this 
region  as  well. 

•  Configuration  Display  —  In  the  lower  left  corner,  shows  the  data  channel  being 
used  as  the  reference  signal  for  all  14  channel  displays,  the  draw  event  interval 
trigger  (axis  crossing  or  signal  peak),  any  adjustments  to  that  interval  (skip¬ 
ping),  and  whether  the  intensity  of  the  bars  is  scaled  to  the  reference  signal. 

‘A  calculated  value  based  on  the  original  sampling  rate  and  the  number  of  samples  displayed 
up  to  that  point. 
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damped  to  a  maximum  value,  or  thresholded  at  a  displayed  reference  signal 
value. 

•  Data  File  Comments  —  In  the  center  of  the  screen,  shows  the  comments  man¬ 
ually  added  to  the  data.h  file  at  the  times  specified.  Sec  Section  G.4.2.3  for 
further  information. 

G  4  Programming  TOPOS 

G.4.1  Prerequisites.  The  program  requires  the  following  files  to  compile: 

topos.e 

topos.h 

tdraw.c 

tmenu.e 

teonf  role 

tcvent.c 

t  event,  h 

The  makefile  included  in  Appendix  1.7  turns  these  files  into  the  executable  code 
topos.  Each  program  file  is  a  functional  grouping  of  functions  coded  in  Silicon 
Graphic's  implementation2  of  C  and  was  executed  on  their  1RIS-4D  workstation 
without  parallel  processing  optimization. 

G-4-2  Program  Flow.  While  the  functional  relationships  in  the  TOPOS  pro¬ 
gram  could  be  derived  from  the  code  in  Appendix  I,  this  section  and  the  flow  chart  in 
Appendix  H  are  intended  to  provide  a  familiarization  of  how  the  program  operates. 
Figure  28  is  a  very  high-level  diagram  of  the  relationships  within  TOPOS. 

G. 4-2.1  Initialization.  Invocation  calls  the  function  main( )  which  ini¬ 
tializes  the  TOPOS  window  and  graphics  and  displays  the  initial  screen  before  hand- 

2The  C  language  used  by  Silicon  Graphics  based  ci,  the  language  as  defined  by  Kerncy  and 
Ritchie's  1978  text,  The  C  Programming  Language.  (40:  A.- 1). 
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Figure  28.  Functional  Relationships  Within  TOPOS.  Primary  function  names  ap¬ 
pear  in  italics.  A  detaiied  How  chart  is  in  Appendix  H. 

ing  over  control  to  the  function  event ().  This  function  constantly  “listens"  for 
mouse-clicks  and  other  events  and  then  performs  the  task  assigned  to  that  event. 
An  example  is  the  pop-up  main  menu  that  appears  when  the  right  mouse  button  is 
depressed. 

The  main  menu  is  drawn  by  do  jnenus  ( 5  and  is  the  gateway  to  the  remainder  of 
the  program.  Each  menu  option  calls  a  new  function,  sometimes  passing  a  parameter 
as  well. 


G. 4-2.2  Display  Options.  All  of  the  display  parameters  listed  in  Sec¬ 
tion  G.3.4  are  controlled  by  separate  functions,  as  shown  in  the  non-exhaust ive  list 
below.  The  program  default  settings  are  determined  in  tcontrol .  c.  These  functions 
assign  numerical  and/or  boolean  states  to  specific  variable  which  are  then  used  in 
the  draw  function  of  TOPOS. 
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•  set_b()  —  Assigns  the  reference  channel  number  to  the  variable  b  (brightness 
driver). 

•  set.lvQ  —  Assigns  the  selected  number  of  pixels  to  the  variable  lw  (fine 
width). 

•  set_bins()  —  Assigns  the  selected  number  of  bins  to  use  in  the  time-shifted 
display. 

•  set.intensityO  —  Determines  if  the  displayed  bar  intensity  will  be  constant 
or  scaled  to  the  reference  signal’s  amplitude. 

•  change_display( )  —  Determines  whether  the  display  will  be  time-shifted,  if 
reference  channel  events  will  trigger  draw  function  calls,  and  whether  those 
events  will  be  peaks  or  axis  crossings. 

•  zero.ptrsO  /  leave_ptrs()  /  positions-ptraO  —  Controls  the  point  in 
the  input  data  file  where  the  program  will  begin  to  read  in  data 

•  set_iuterval_skip()  —  Controls  the  interval  between  draw  events  by  assign¬ 
ing  the  number  of  reference  channel  peaks  or  axis  crossings  to  skip  and  by- 
assigning  whether  the  midpoints  between  such  signal  events  will  also  trigger 
draw  events,  effectively  halving  the  time  between  events. 

•  set.testO  —  Turns  on  and  off  the  test  modes  of  single  buffer  drawing,  step 
drawing,  and  batch  drawing. 

•  svapO  —  Causes  the  back  and  front  display  buffers  to  swap  continuously  until 
an  event  occurs,  such  as  a  mouse  click. 

G.4.2.3  Converting  EEG  Data.  Selecting  the  menu  option  “Prepare 
Data "  starts  the  function  prep.dataO,  which  in  turn  calls  other  functions  such  as 
read_hdr()  and  read.chO  to  implement  the  techniques  discussed  in  Appendix  A 
for  converting  the  Drain  Atlas  EEG  files  into  a  form  usable  by  TOPOS.  The  output 
is  written  to  separate  files,  one  for  each  EEG  channel  as  well  as  one  containing  the 


96 


maximum,  minimum  and  standard  deviation  of  amplitude  values  of  each  channel. 
This  latter  file  has  the  same  name  as  the  input  file,  but  with  an  additional  “.h" 
extension.3  The  remaining  output  files  are  marked  similarly  with  “.chXX"  extensions, 
corresponding  to  their  channel  number. 

G.4.2.4  Displaying  Data.  After  initially  drawing  the  display's  back¬ 
ground.  the  function  drav_display()  is  a  loop  with  multiple  conditional  statements. 
The  execution  of  this  function  is  controlled  by  the  software  switches  provided  by  the 
functions  in  Section  G.4.2.2  as  draw.displayO  works  through  the  data  files  one 
sample  at  a  time.  Because  there  are  numerous  possible  combination  of  switch  set¬ 
tings,  ttiere  are  many  possible  paths  through  this  function.  However,  only  a  single 
software  path  which  was  used  often  during  testing  is  presented  below.  The  descrip¬ 
tion  provided  will  cover  the  more  basic  components  common  to  most  paths  within 
draw.displayO ,  as  well  as  the  functions  of  the  specific  switch  options  involved. 

TOPOS  Configuration:  Reference  CH:  Ch  2 

Bar  Width:  5  pixels 

Bar  Intensity:  Scaled 

Display:  Time-shifted  ‘Of  axis  crossings 

Bin  Number:  10 

Test  Modes:  Off 

Start:  Beginning  of  Data 

1.  Loop  from  first  sample  in  the  data  files. 

2.  Break  loop  if  mouse  click  event  has  occurred. 

3.  Find  successive  data  samples  in  Channel  2  which  change  signs.  Remainder  of 
the  function  will  use  the  first  value  of  this  data  pair,  corresponding  to  sample 
s-1  in  the  program. 

3Data  file  comments,  if  they  are  desired  on-screen,  must  be  manually  added  to  this  file  following 
the  channel  parameters  and  in  three  fields:  1 )  minute  time- tag  (integer);  2)  seconds  time- tag  (float): 
3)  comment  (alpha-numeric,  no  spaces,  less  than  21  characters). 
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4.  Swap  display  buffers  to  bring  current  display  to  the  back  for  editing  and  to 
show  the  most  recently  drawn  display. 

5.  erase_bars()  —  Erase  the  bar  segment  previously  drawn  in  the  time  bin 
corresponding  to  the  next  draw  event. 

C.  draw_boxes()  /  draw.statusO  —  Redraw  the  channel  boxes  and  configura¬ 
tion  and  status  boxes. 

7.  draw_bars()  —  Set  line  width  to  5  pixels. 

8.  draw-bars ()  —  As  required,  clip  any  channel  value  at  sample  s-1  to  three 
standard  deviations  above  the  channel's  mean. 

9.  draw_bars()  —  Calculate  the  vertical  displacement  of  the  bars  to  be  drawn 
from  each  channel's  value  at  sample  s-1,  scaling  the  maximum  possible  dis¬ 
placement  to  three  times  the  standard  deviations  of  each  channel. 

10.  draw.baraO  —  Calculate  the  start  and  stop  coordinates  of  the  bar  being 
drawn,  time-shifted  to  the  right  from  the  previously  drawn  bar. 

11.  draw-bars  ()  —  Set  the  draw  color  to  the  gray-scale  value  calculated  from  Ch 
2’s  amplitude  at  sample  s-1,  scaling  the  maximum  possible  displacement  to 
three  times  the  standard  deviation  of  each  channel. 

12.  drav_bar8()  —  Draw  a  5-pixel-wide  line  (the  ;bar’’). 

13.  draw.baraO  —  Repeat  loop  until  interrupted  or  at  end  of  data.  Reset  bin 
number  to  0  every  time  Bin  9  is  drawn. 

G.4-2.5  Obtaining  Screen  Prints.  Screen  images  of  TOPOS  displays 
were  captured  using  the  Silicon  Graphics  tool  snapshot  in  a  raster  scan  format. 
These  were  first  converted  to  run-length  encoded  format  using  the  program  tourt 
and  then  to  enca  mlated  postscript  files  using  rle2eps. 
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Table  6.  TOPOS  Test  File  Contents  by  Channel. 


CH 

nevtest . dat 

vartest . dat 

1 

0.5  Hz 

1  Hz 

2 

1  Hz 

2  Hz  | 

3 

2  Hz 

3  Hz 

4 

1  Hz 

4  Hz 

5 

2  Hz 

1-4  Hz  (varying) 

6 

3  Hz 

Sum  (CHs  1-4) 

mm 

4  Hz 

Sum  (CHs  1,2,4) 

8 

Sum  (CHs  1-3) 

Sum  (CHs  1,2, 4, 5) 

9 

Sum  (CHs  4-6) 

Sum  (EEG,  CH  5) 

10 

CH  4  with  45  deg  phase-shift 

Sum  (EEG,  SHI 

11 

CH  4  with  90  deg  phase-shift 

Sum  (EEG,  CH  1) 

12 

CH  6  with  40  deg  phase- shift 

Sum  (EEG, 

13 

CH  6  with  90  deg  phase-shift 

KiBSI!: 

14 

1.5  Hz 

EEG 

G.4-3  Test  and  Diagnostic  Files.  Table  6  shows  the  content  of  each  test  file 
used  during  program  development  and  included  with  the  executable  and  source  code 
upon  thesis  completion. 

In  addition,  a  diagnostic  file  listing  is  generated  for  every  TOPOS  session.  The 
file’s  name  is  the  data  file  with  a  .diag  extension.  It  includes  information  regarding 
each  draw  event:  sample  number,  displacements  by  channel,  and  other  parameters. 

G.4.4  Program  Bugs.  The  following  are  known  problems  which  the  creator 
either  did  not  know  how  or  had  insufficient  time  to  overcome. 


•  User  loses  control  over  program  execution  if  the  window  is  repositioned,  resized, 
or  collapsed  to  an  icon  while  the  display  is  not  showing  “Display  Halted”  over 
the  Status  box.  Requires  user  to  kill  process  from  console. 

♦  All  data  display  is  erased  when  window  is  resized  or  redrawn  (i.e.,  if  user 
displays  another  application  over  it  and  then  quits  or  moves  that  new  appli- 
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cation's  window).  However,  data  file  pointers  have  not  moved  and  display  can 
be  restarted  using  any  of  the  three  “Start  Display’’  menu  options. 

•  No  more  than  10  comments  can  be  read  from  the  data.h  file  for  display  during 
TOPOS  execution.  If  fewer  comments  are  desired,  then  large  times  and  empty 
comment  fields  should  pad  out  the  required  10  lines. 
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Appendix  H.  TOPOS  Flow  Chart 


Figure  29.  TOPOS  Flow  Diagram  (Part  1  of  4) 
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Prepare  Data 
prep.dataQ 

Reference  Ch 
set_b() 


Bar  Width 
set_lw() 


Bar  Intensity 
set_intensity() 

Display  Mode 
change_display() 


- 

Draw  Interval 
sotJnterval skip() 

- 

Help 

help() 

- 

Bin  Number 
set_bins() 

- 

Test 

set_test() 

Start  Display  @ 

Beginning 

zero_ptrs() 

Current  Pos’n 
leave_ptrs() 

Selected  Time 
position_ptrs() 

- 

Swap  Buffers 
swap() 

-*► 

Quit 

quit() 

Figure  30.  TOPOS  Flow  Diagram  (Part  2  of  4) 
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draw_display() 


draw_title() 


prep_data() 
been  run? 


open  input  files  & 
read  .h  file 


open  fixed-frea  reference  files 


read  data  file  comments 


BEGIN  DRAW  LOOP 


desired  start  time  reached? 


read  &  disreaard  data 


time  for  next  data  comrnentT^>  display  next  comment  1 


stop  display 


rt  or  left  mouse  pushed? 


continuous  or  step  draw  & 
jriiddle  mouse  held  down Z 


pause  display 


read  sample  from  each  file 


mutton  released 'i 


Figure  31.  TOPOS  Flow  Diagram  (Part  3  of  4) 
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swapbuffersQ 


erase_bars() 


dravy_boxes() 


draw_status() 


draw_bars()  Z>-0  (Return  to  start  of  Draw  Loop ) 


Figure  32.  TOFOS  Flow  Diagram  (Part  4  of  4) 


Appendix  I.  TOPOS  Code 


I.l  Main  File:  topos.c 


/* 

*  topos.c 

• 

*  The  main  file  for  the  AFIT  toposcope  TOPOS. 

* 

*  Capt  Dwight  Roblyer,  AFIT/GS0-92D ,  November  1992 

* 

*  Built  over  the  program  "CURVE  DEMO" 

*  by  Howard  Look  for  Silicon  Graphics,  June  1989 

* 

*/ 


•include 

•include 

•include 

•include 

•include 

•include 

•include 

•include 


<stdio . h> 
<string.h> 
<gl.h> 
<device.h> 
<math.h> 
<fmclient .h> 
"tevent.h" 
"topos . h" 


/**•*  GLOBAL  VARIABLE  DECLARATIONS  •***/ 


FILE 

•ifp; 

/* 

pointer  to  ".dat"  input  file 

*/ 

FILE  *ofp[19]  ; 

/* 

".dat.chXX"  files 

*/ 

FILE 

•hfp; 

/• 

".dat.h"  file 

*/ 

char 

AD, 

/* 

holding  var  for  read-in  hdr  data 

*/ 

caluv, 

/* 

calibration  volt  in  microvolts  (from 

i  hdr) 

*/ 

sens , 

/* 

sensitivity  (from  hdr) 

*/ 

bb, 

/* 

bit  bucket  var  for  unwanted  data  (from  hdr) 

*/ 

cv[15] , 

/* 

cal  volts  in  A-D  cr mts  by  ch  (from 

hdr) 

*/ 

dc [15] , 

/* 

DC  offsets  by  channel  (from  hdr) 

*/ 

v [15] ; 

/* 

holding  var  for  read-in  EEG  data 

*/ 

int  i=0, 

/* 

indicator  for  byte  number  of  a  value 

*/ 

ch“l , 

/* 

channel  counter 

•/ 

k, 

!* 

counter  for  •  time  thru  main  drawing 

loop  */ 

8, 

/* 

counter  for  •  of  EEG  samples 

*/ 

C  , 

/* 

variable  used  to  chk  for  EOF 

•/ 
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minutes=0,  /*  *  minutes  elapsed  (for  display)  */ 
hardware,  /*  type  of  hardware  running  */ 
origin.*,  origin.y,  /*  Lower  left  corner  of  the  window  */ 
size.*,  size.y;  /*  Size  of  the  window  */ 


float  seconds=0., 

data[18] ,  /*  "volt"  is  a  15-element  real  number  array  */ 

/*  used  to  hold  the  calculated  voltage  for*/ 

/*  each  channel  (Oth  element  is  not  used)  */ 

max[19],  /*  maximum  value  in  each  channel  */ 

mini [19],  /*  minimum  value  in  each  channel  */ 

sum.of _sqs  [15] ,  /•  sum  of  squared  values  in  each  channel  */ 

std_dev[l9],  /•  std  deviation  in  each  channel  */ 

t=0.,  /*  tracks  elapsed  time  of  display  */ 

aspect;  /*  Window’s  aspect  ratio,  x/y  */ 


char  sample_rate_char=0 ,  /*  header  value  indicating  record  sample  rate  */ 

*newfilel,  /*  filename  transfer  variable  */ 

*newfile2,  /*  filename  transfer  variable  */ 

xfer[25] ,  /*  this  worked  (*xfer  didn’t)  for  sprint  (see  draw.status)*/ 

timel[50];  /*  elapsed  time  status  line  */ 


Boolean  noprep ,  /•  tracks  whether  prep.dataO  has  been  run  •/ 

draw.active ;  /•  Is  drawing  active  ?  (a  holdover  from  CURVE  DEMO)  */ 


Coord  half,  offset,  text.offset;  /*  length  offsets  for  text  •/ 


/*•**  FUNCTION  PROTOTYPES  •*•*/ 


void  prep.data(void) ; 

void  readhdr(void) ; 

void  readch(char  cv[],char  dc[] 

void  time.print(void) ; 

void  init_window(char  *argv[]); 

void  init.events(void) ; 

void  redraw(void) ; 

void  process.lef t.down(void) ; 

void  process.lef t.up(void) ; 

void  process.middle. down (void) ; 

void  process. middle. up(void) ; 

void  draw.instructions(void) ; 

void  quit(void) ; 

void  init.fonts(void) ; 

void  help(void) ; 

void  end.help(void) ; 

Boolean  helping(void) ; 


/•  calls  readhdrO  and  readchO  (see  below)  */ 
/*  reads  in  data  from  record  header  */ 

,char  sens, char  caluv) ; 

/*  reads  data  in  counts  and  converts  to  volts*/ 
/*  prints  elapsed  time  to  status  box  */ 

/*  initializes  window  */ 

/*  defines  valid  actions  for  event  queue  */ 

/*  redraws  display  when  window  is  redrawn  */ 

/*  links  event  to  action  »/ 

/*  links  event  to  action  */ 

/•  links  event  to  action  •/ 

/*  links  event  to  action  */ 

/*  draws  help  screen  */ 

/*  quits  program  gracefully  */ 

/*  defines  text  font  */ 

/*  displays  help  screen  */ 

/*  removes  help  screen  */ 


/ 


•  tracks  presence  of  help  screen  */ 


extern  void  init.menus(void) ;  /*  creates  menus  */ 
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extern  void  draw.display (void) ;  /*  draws  toposcope  display  */ 


/•*•••••• . . . . .  MAINO  . **************************,*./ 

void  main(int  argc ,  char  *argv[]) 

{ 

if  (argc  ==  1) 

{ 

printf("\n  USAGE:  topos  inf ile . dat\n") ; 
quit() ; 


strcpy(xfer ,argv[l] ) ; 
drav_active  «  FALSE; 


/•*•*••**•  initialize  everything  ••*****•**/ 
init_window(argv) ; 

frontbuffer(TRUE) ;  /*  Enables  writing  to  front  buffer  */ 

/*  (Default  is  only  back  is  enabled)*/ 

czclear(Oxcccccccc ,0x7FFFFF) ;  /*  Clears  the  viewport  to  It  gray  •  / 

/*  and  sets  the  z-axis  to  minimum  */ 

/*  (0x7FFFFF)  in  enabled  buffers  */ 

frontbuffer (FALSE) ;  /*  Disables  writing  to  front  buffer*/ 

/*  (back  buffer  still  enabled)  */ 

init.fontsO  ; 
init_menus( )  ; 
init_events() ; 
draw, instruct ions () ; 

while(TRUE)  /*  Continuously  calls  function  which  */ 

eventO;  /*  monitors  for  mouse  actions,  etc.  */ 


> 

. . *** . ******  PREP.DATAO  ************************** . .*/ 


void  prep.data(void) 

{ 
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noprep  *  FALSE; 

if  (fopec(xfer,  "r")  ='  MULL) 

{ 

printf("\n  Input  file  */.s  does  not  exist  in  this  directory\n" , xfer ) ; 
quitO ; 

> 

else 

ifp  *  fopen(xfer,  “r"); 


/******••*••*•  create  ".dat.cbXX"  files  *••••********•*/ 

for  (ch=l ;ch<=l4;++ch) 

newf lie  1  =  "MAXIMUWF1LENAME . EXTENSION" ; 
nevf ile2  “  "MAXIMUMFILENAME. EXTENSION"; 
strcpy(newfilel ,xfer) ; 
strcpy(newf ile2 ,xfer) ; 

if  (ch==l)  strcat (newf ilel chi") ; 
if  (ch==2)  strcat(newf ilel , " . ch2") ; 
if  (ch«3)  strcat(newf ilel , " . ch3") ; 
if  (ch*«=4)  8ticat(netffilel,".ch4") , 
if  (ch==5)  8trcat(newf ilel , " . ch5"; ; 
if  (ch==6)  atrcat (newf ilel, ".ch6"); 
if  (ch==7)  etrcat(nevf ilel , " . ch7") ; 
if  (ch==8)  8trcat(neufilel,".ch8") ; 
if  (ch==9)  dtrcat(newf ilel . ch9") ; 
if  (ch==10)  strcat(newfilel,".chlO") ; 
if  (ch==ll>  atrcat(nevf ilel . chll") ; 
if  (ch==12)  strcatCnewf ilel ,chl2") ; 
if  (ch==13)  strcat(newf ilel .chl3") ; 
if  (ch==14)  strcat (newf ilel chl4") ; 

ofp[ch]=fopen(newfiiel ,  "v+">; 

max[ch]=0. ; 
mini  [ch] =0. ; 


atrcat(ne¥file2,'!.h") ; 
hfp=f open (newf  ile2,  '«+"); 


readhdr(,  ; 

f sec* (  i  .1310,0;,  /*  Sets  the  fiie  pointer  to  start  of  */ 

/*  the  EEG  data  (1920th  byte)  from  the*/ 
/*  beginning  of  the  file  (byte  0)]  */ 
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for  (ch=l ;ch<=14 , ++ch) 
sua.of _sqs [ch]  =d .  0 ; 


/*  Resets  array  used  in  Std  Dev  calc  */ 


k=l; 


while  ((c  =  getc(ifp)1  EOF) 

{ 

s=k-l ; 

readch(cv,  ,  sens,  caluv) ; 

♦+k; 

\ 

for  (ch=l ;ch<=14 ; ++ch) 

atd_dev[ch]  *  sqrt (sum.of _sqs [ch] /k) ; 

for(ch=l ;ch<=14 ;**ch) 

fprintf (hfp,"%f  Xf  ‘/.i  ‘/.f\n"  ,max  [ch]  .mini  [ch]  ,  sample. rate. char , 

std.dev [ch] ) ; 


/•••••••••a*********************  READHDRO  *••**•*••••*•****••**•*••*•****•*/ 

void  readhdr(void) 

{ 


/♦read  in  sampling  rate  of  AD  conversion*/ 
f seek(iip,416 ,0) ; 

fscanf (ifp/'Xc" ,ksattple_rate_char) ; 

if  (aajnple_rate.char=B7)  sanple_rate_char=128; 
if  (saniple_rate_char3S15)  s ample. rat e_chax=64 ; 


/•read  in  calibration  voltages*/ 

for  (ch=l;ch  <=  14;**ch) 

{ 


if  (ch  ==  1) 

feeek(ifp, 428,0) ; 

if  (ch  ==  3) 
f seek ( ifp ,1,1); 
if  (ch  ==  7) 

f seek(ifp, 1 , 1) ; 

if  (ch  ==  11) 


/•  Sets  "ifp"  to  byte  428  from  file  s^-art  (0)*/ 

/*  byte  431  (skipping  430)  •/ 

/*  byte  436  (skipping  435)  */ 

/*  byte  441  (skipping  440)  */ 
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fseek(ifp, 1 , 1) ; 

if  (ch  **  14)  /*  byte  445  (skipping  444)  •/ 

fs8ek(ifp, 1 , 1) ; 

fscanf  (ifp,  kAD) ; 

cv[ch] =AD; 


/•read  in  dc  offsets*/ 
for  (ch*l;ch  <°  14;*-*-ch)  { 


if  (ch  ==  1) 

fseek(ifp, 449,0) ; 
if  (ch  ==  3) 
fseek(ifp, 1 , 1) ; 

if  (ch  **  7) 

fseek(ifp, 1 , 1) ; 


if  (ch  ==  11) 
fseek(ifp, 1 , 1) ; 

if  (ch  ==  14) 
fseek(ifp, 1 , 1) ; 


/*  Sets  "ifp"  to  byte  449  from  file  start  (0)  */ 
/*  byte  452  (skipping  451)  */ 

/*  byte  457  (skipping  456)  •/ 

/*  byte  462  (skipping  461)  */ 


/•  byte  466  (skipping  465)  */ 


fscanf (ifp,  "%c",  kAD); 
dc  [ch]=AD; 


> 

fscanf  (ifp,  "7,c",  ksens) ; 

fseek(ifp, 1 , 1) ;  /•  byte  469  (skipping  468)  •/ 

fscanf (ifp.  "%c",  kcaluv) ; 


READCHO  •**•••••*•*•*********••*«>*****•**/ 

void  readch(char  cv[],  char  dc[],  char  sens,  char  caluv) 

{  /*  readin  one  sample  of  all  14  channels  */ 

f  scanf  ( if  p ,  "%cXclclcXc%cy,c%c%aac%cXcUlcy,cXcy.Mc  " , 
kbb,kbb,kv[l] ,kv[2] ,kbb,kv[3] ,kvt4] ,kvl5] ,kv[6] ,kbb, 
kv [7] ,kv [8] ,iv [9] ,kv[lO] ,kbb,kv[ll] ,kv[12] ,kv[13] ,kv[14] ,kbb) ; 

/*  Reads  in  blocs  of  20  CHAR  (1  byte)  A-D  values  into  */ 
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/•  either  the  bit  bucket  variable  or  the  CHAR  array  */ 
/*  "v"  according  to  the  format  found  in  the  lab  */ 

/*  convert  to  voltages  •/ 
f  or(ch=l ;  ch<=14  ;-*-+ch) 

{ 

data[ch]  =  ((float)v[ch]-(float)dc[ch])*(2.0*caluv/cv[ch]  )* (sens/ 10.0)  ; 

if  (data[ch]  >  max[ch]) 

max[ch]  =  dataCch] ; 

if  (data[ch]  <  minifchl) 

ninitch]  =  data[ch]; 

/*  write  each  channel’s  data  to  separate  files  •/ 
fprintf  (ofp[ch]  ,  "y.f\n"  ,data[ch]  ) ; 

su«_of _3qe [ch]  *  sun. of _sqs [ch]  +  pow(data[ch] ,2) ; 


> 


> 


/•*•*•*••«»••••••■••••••••••*•**  INIT.WINDOWO  ****♦***♦**»♦***•♦*•**•»*•*****/ 

/* 

•  Opens  the  window  and  sets  up  the  graphic  node. 

*/ 

void  init_windov(char  »nrgv[]) 

{ 

char  *t,  estrrchrO; 
int  planes ; 

keepaspect( .8b,l) ;  /*  x/y  ratio  of  window  */ 

/*  open  the  window  with  its  invocation  as  the  title  (less  the  path)  */ 
winopen((t=etrrchr(argv[0]  ,  ’/’))  !=  MULL  ?  t«-i  :  argv[0]); 

/• 

*  If  we  run  on  an  8-bit  eclipse,  we  don’t  use  RGB  node  or  depth 

*  cuing.  Speed  things  up  on  all  Eclipses. 

•/ 

planes  *  getplanesO  ; 

if  (planes  <■  8) 

hardware  ■  ECLIPSE8; 
else  if  (planes'  <=  24) 

hardware  *  ECLIPSE24; 

else 

hardware  »  GT; 
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doublebuffer ()  ; 


if  (hardware  !■  ECLIPSE8) 

RGBmode  0 ; 

zbuf fer(TRUE) ; 
gconf  igO ; 

getorigin(fcorigin_x,  Aorigin.y) ; 
get8ize(tsize_x ,  tsize.y); 

aspect  =  (f loat)aize_x/(f loat)aize_y ; 

> 


. . . . .  INIT.EVENTSO  **.•*••••*•****•*•**********•**/ 

/* 

•  Tells  event  manager  vhat  to  pay  attention  to. 

*/ 

void  init.events(void) 

{ 

/*  ESC  key  and  WINQUIT  quit  the  program  */ 
add.e vent (ANY,  ESCKEY ,  UP,  quit,  0); 
qdevice(ESCKEY) ; 

add_event(ANY,  WINQUIT,  ANY,  quit,  0); 
qdevice(WINQUIT) ; 

/*  vinduw  redraw  events  •/ 
add.eventCANY,  REDRAW,  ANY,  redraw,  0); 
qdev ice (REDRAW) ; 

add_event(ANY,  LEFTMOUSE,  DOWN,  process.left.down,  NULL); 
add_event(ANY,  LEFTMOUSE,  UP,  proceos.left.up,  NULL); 
qdevice (LEFTMOUSE) ; 

add_event(ANY,  MIDDLEMOUSE,  DOWN,  process.middle.down,  NULL); 
add.e vent (ANY,  MIDDLEMOUSE,  UP,  proceBs.middle.up,  NULL); 
qdevice (MIDDLEMOUSE) ; 

/•  Draw  the  display  whenever  the  user  does  not  request  anything  */ 

/•  i.e.,  there  are  no  events  happening  */ 
add_update(Adr aw. active,  draw. display ,  NULL); 

> 


/•*•****•***•*••**••••••••*•*•*•**»  REDRAW ()  **.*********•*****.•***********/ 

/• 

*  Called  by  the  event  manager  whenever  a  redraw  event  occurs, 
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•  e.g.  when  the  window  ia  resized. 

*/ 

void  redraw(void) 

{ 

getorigindtorigin.x,  korigin.y) ; 
getsize(*size_x ,  fcaize.y) ; 
aspect  =  (float)aize_x/{float)Bize_y ; 
viewport(0,  size_x-l,  0,  size.y-1); 

if  (helpingO) 
helpO  ; 

else 

{ 

frontbuffer(TRUE) ; 

czclear (Oxcccccecc ,0x7FFFFF) ; 

irontbuffer(FALSE) ; 

drav..title() ; 
drav_boxes() ; 
draw_status( ) ; 


/*«,4.*******.**.*****.*****.******0THER  FUNCTIONS*************1  *  ',***r*: . . 

void  process. left. down(void) 

{ 

int  marker; 

if  (helpingOH 

end.LelpO ; 
return; 

> 

> 


void  proceS8_left_up(void) 

< 

onqdevice  (MOUSED  ; 
unqdevico (MOUSEY) ; 

> 


void  proce88_niddle_down(void) 

{ 

int  marker; 
end.helpO  ; 

> 
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void  process_middle_up(void) 

{ 

unqdevice(HOUSEX) ; 
unqdevice (MOUSEY) ; 

} 


/*  limits  of  the  view  volume  */ 
Coord  box.limit  [3]  [2]  = 

{  {  -1.0,  1.0>, 

{  -1.0,  1.0), 

<  -1.0,  1.0> 

> ; 


void  quit(void) 

{ 

gexitO  ; 
fclose(ifp) ; 

for  (ch=l;ch<=14;++ch) 
fclo8«(ofp[ch] ) ; 

fcloee(hfp) ; 

exit (0) ; 

> 


/* 

*  Routines  for  writing  the  intro  messages  to  the  screen. 

*/ 


double  min(double  x,  double  y) 

{ 

if  (x  <  y) 

return  x; 

else 


> 


return  y; 


static  fmfonthandle  base. font,  scaled.font; 


void  init.fonte(void) 

{ 

fminitO  ; 

base.font  ■  fmf indfont("Times-Romac") ; 

> 
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void  setup_string(char  str[100],  int  size,  Coord  *half) 

{ 

double  h.size,  v.size,  f.size; 

Coord  length; 

h.size  =  (double)size_x/(double)size ; 
v.size  =  (double)size_y/(double)aize; 
f.size  =  min(h_size,  v.size); 
scaled.font  =  fmscalefont(base_font ,  f.size); 
fmsetfont (scaled.font) ; 

length  *  ftngetstrwidth(sealed_font ,  atr) ; 

•half  =  (Coord) (length/aize.x) ; 

> 

static  int  looking.at.help; 
static  int  active.state ; 

void  help(void) 

{ 

if  ( ! looking.at.help) 

( 

looking.at.help  »  TRUE; 

active.state  ®  drav.active; 
draw.active  *  FALSE; 

> 

drav. instruct ions () ; 

> 

void  end.help(void) 

{ 

if  (looking.at.help) 

{ 

looking.at.help  *  FALSE; 
drav.active  »  active.state; 

> 

> 

Boolean  helping(void) 

{ 

return  looking.at.help; 

> 

void  time.print(void) 

{ 

ortho2(-l .0,  1.0,  -1.0,  1.0); 
depthcue (FALSE) ; 
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cpack(OxaaOOOO) ; 


■inutes“(int) (t/60) ; 
seconda=t-( (float)minutes*60) ; 

sprintf  (t.imel ,  "Elapsed  Time:  7,i:X.3f  sec",  minutes  .seconds) 
aetup_string(timel ,  70, khalf); 
cmov2( .58,-0.80) ; 
fmprstr(timel) ; 


void  draw.inatructionsivoid) 

{ 

Coord  half,  offset,  text.offset; 
char  str[l00]; 

ortho2(-1.0,  1.0,  -1.0,  1.0); 
czcle&r  (Oxaaaaaaaa ,0x7FFFFF) ; 
depthcue (FALSE) ; 

sprintf (atr ."Topos") ; 

cpack (0x00655565) ; 

8«tup_Btring(str , 10 .khalf ) ; 
caov2(0 .006  -  half,  0.706); 
fmprstr(atr) ; 

/*  set  draw  color  to  BLUE  */ 
if  (hardware  **  ECLIPSE8) 
color(BLUE) ; 
elae  cpack(OxbbOOOO) ; 

setup_string(str , 10 .khalf ) ; 
cmov2(0.0  -  half,  0.7); 
fmprstr(8tr) ; 

/•  set  draw  color  to  BLUE  */ 
if  (hardware  «  ECLIPSE8) 
color(BLUE) ; 

else 

cpack(0x880000) ; 

8printf(str, "Interactive  AFIT  Toposcope"); 
setup_string(str, 30, khalf) ; 
cbov2(0.0  -  half,  0.5); 
faprstr(Btr) ; 

/*  set  draw  color  to  WHITE  */ 
if  (hardware  »■  ECLIPSE8) 
color(WHITE) ; 
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else 


cpack(0x222222) ; 


/*  just  sets  font  size  for  body  of  help  screen  •/ 
setup_string(st.r  ,50,thalf) ; 

offset  *  -0.6; 
text.offset  =  -0.30; 

aprintf (str, "Mouse  Left:"); 
cmov2(off3et,  0.3); 
fmprstr (atr) ; 

sprintf (str ."Halt  display."); 
cmov2(text. offset .  0.3); 
fmprstr(atr) ; 

aprintf (atr, "(Requires  V'Start  DisplayV  from  menu  to  restart.)"); 
cmov2(text_offset,  .22); 
fmpratr(str) ; 

aprintf (str, "House  Middle:"); 
cmov2(offset,  0.10); 
fmprstr (str) ; 

sprintf (str, "Continuous  Drawing  -  Press  pauses,  release  restarts."); 
cmov2(text_offset ,  0.10); 
fmprstr(str) ; 

aprintf (str, "Step  Drawing  -  Press/release  both  draw  next  step."); 
cmov2 (text.offset ,  .02); 
fmprstr (str) ; 

sprintf (str, "Batch  Drawing  -  Press  draws  next  batch  (redraws  all 
bins) .") ; 

cmov2(text_offset ,  -.06); 
fmprstr (str) ; 

aprintf (atr, "Mouse  Right:"); 
cmov2(offset,  -.18); 
fmpratr(str) ; 
sprintf (str, "Main  menu"); 
cmov2(text.off set ,  -.18); 
fmprstr(str) ; 

swapbuffersO  ; 
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1.2  Screen  Creation:  tdraw.c 


/• 

•  tdraw.c 

* 

*  Part  of  "TOPOS"  by  Capt  Dwight  Roblyer,  AFIT/GS0-92D 

*  November  1992 

* 

*  Baaed  on  "CURVE  DEMO" 

*  by  Howard  Look  for  Silicon  Graphics,  June  1989 

* 

*  This  file  contains  functions  which  draw  the  screen  displays. 

* 

•/ 


♦includa  <atdio.h> 
♦include  <gl.h> 
♦include  <device.h> 
♦include  <math.h> 
♦include  "topos.h" 
♦include  "tevent.h" 


/•  Prototypes  */ 
void  draw_display (void) ; 
void  drau_boxea(void) ; 
void  drav.bars(void) ; 
void  erase.bars(void) ; 
void  draw. atatus (void) ; 
void  draw.title(void) ; 
void  draw.f ilenaae(void) ; 
void  draw.update(void) ; 
void  draw.freq(void) ; 
void  draw.matrix(void) ; 
void  draw.line(void) ; 

extern  void  time.printO  ; 


/«  coordinates  overall  screen  display  */ 

/*  draws  blue  channel  boxes  */ 

/*  draws  bars  displayed  in  each  channel  box  */ 

/*  overwrites  bars  with  the  background  color  */ 

/*  displays  status  and  configuration  info  */ 

/•  draws  title,  electrode  locations/labels/lines  */ 
/*  fills  in  status  for  data  file  name  */ 

/*  fills  in  orig  digitization  rate  of  KG  file  */ 

/*  fills  in  instantaneous  frequency  of  ref  signal*/ 
/*  displays  correlation  matrix  •/ 

/*  line-drawing  function  */ 


int  b, 

bb_ integer , 
bins, 

samplo.rate , 

lv, 

n, 

s.old, 
s.diff , 
shift , 


/*  the  ch  selocted  as  reference  •/ 

/*  bit-bucket  for  integer  data  types  */ 

/*  number  of  time-shifted  bars  before  repeat  */ 

/*  record  rate  of  digitized  EEG  record  •/ 

/*  line  or  bar  width  (or  height)  */ 

/*  counter  used  for  reading  parameters  for  */ 

/*  hardcoded  reference  sinusoids  */ 

/•  sample  number  of  most  recent  event  */ 

/*  number  of  samples  between  two  most  recent  events  */ 
/*  tracks  which  time  bin  currently  in  use  •/ 
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shift.last,  /»  tine  bin  used  last  */ 

pack=0,  /*  used  to  define  color  of  bare  */ 

s=0,  /*  data  sample  currently  being  processed  */ 

cmt.minClO] ,  /«  minute  portion  of  data  comment  time-tag  */ 

cmt_sec[10],  /*  seconds  portion  of  data  comment  time-tag  •/ 

cmt_sample[10]={0},  /*  sample  #  at  which  to  display  each  of  10  comments  */ 

cmt_cnt=0,  /*  counter  to  track  data  file  comments  */ 

flag,  /*  counter  for  data  file  comment  processing  */ 

buffer,  /*  diagnostic:  tracks  display  buffers  */ 

buffer_now;  /*  diagnostic:  holds  current  buffer  */ 

float  bb_float,  /*  bit-bucket  for  decimal  data  types  */ 

packd=0.,  /*  used  to  define  color  of  bars  */ 

data_old[l9] ={0 . > ,  /*  data  pt  in  each  ch  read  last  iteration  »/ 

data_oldold[19]={0.},  /*  data  pt  in  each  ch  read  two  iterations  ago  */ 

startClS] [76] [3]={0>,  /*  XYZ  of  start  of  bar  by  channel  and  shift  */ 
start_old[15] [76] [3]={0>,  /*  last  value  of  "start"  */ 

stop[16] [76] [3]={0} ,  /*  XYZ  of  bar  stop  by  channel  and  shift  •/ 

stop_old[l5] [76] [3]={0>,  /*  last  value  of  "stop"  */ 

scale. factor[19]={0.},  /*  3  stand,  dev.  scale  for  each  channel  */ 

correl[15]={0.},  /*  channel  correlations  for  matrix  display  */ 

cursor_top[16] [3] ,  /*  coord  of  cursor’s  upper  pt  in  each  ch  */ 

cursor_bottom[15] [3] ,  /*  coord  of  cursor’s  lower  pt  in  each  ch  */ 

cursor_top_old[15] [3] ,  /•  last  value  of  "cursor.top"  */ 

cursor_bottom_old[15]  [3] ,  /*  last  value  of  "cursor.bottom  •/ 

ypos[19]={0> ,  /•  amplitude  of  each  channel  at  time  t  */ 

freq=0 .  ,  /•  calculated,  instant  freq  of  reference  signal  •/ 

diap[19] ,  /*  displacement  of  bar  from  x-axis  by  channel  */ 

ptl [3] ,pt2[3] ;  /*  used  for  electrode-to-box  lines  •/ 

FILE  *tfp ;  /*  file  ptr  to  datafile .DIAG  output  file  */ 

Boolean  emt_flag[20] ;  /•  holds  which  display  comment  has  been  displayed  */ 

unsigned  long  packcolor [15] ;  /*  special  variable  type  used  to  hold  */ 

/*  the  "packed  integer"  required  by  •/ 

/*  the  function  "cpack"  */ 

char  title [20],  statue [50],  electrode [20] ,  channel [10]; 

/•  strings  for  drawing  text  to  screen  */ 

char  cmt[20][20];  /*  strings  for  drawing  data  file  cmnts  */ 

extern  Boolean  time.shift, 
interval , 
axis. crossing, 
peak, 

constant.intensity , 
zero.ptrs , 
wait , 

8ingle_buff , 
synch. constrained , 
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batch; 


extern  int  hardware, 

total.skip , 

skip.seconds , 

ch, 

c, 

k, 

interval. skip; 

extern  float  t, 

max[] .  mini[]  , 
atd.dev  [] , 
data[] , 
volt  []  , 
synch.limit; 

extern  char  cv[] , 
dc[]  , 
sens , 
caluv, 
xfer[] , 

•newfilel, 

•neufile2; 

extern  FILE  *ofp[] , 

•hfp; 

extern  Boolean  noprep, 

zero.pointera, 
aynch.conatrained ; 

extern  Coord  half ; 


. . .  FUNCTION  "DfUW.DISPLAY"  . . . 

void  draw.diaplay(void)  /*  Function  used  to  draw  the  screen  */ 

{ 


qresetO  ; 
ahift=0; 

frontbuffer(TRUE) ;  /*  Enables  front  k  back  buffers  for  drawing  •/ 

czclear (Oxcccccccc ,0x7FFFFF) ;  /*  Re-clears  the  viewport  to  gray  and  */ 

/*  sets  the  z-axis  to  nininuo  0x7FFFFF  */ 


aetlinestyle (SOLID) ; 
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frontbuffer(FALSE) ; 


/*  Disables  front  buffer  for  drawing  */ 


dxaw_title() ; 


if  (noprep)  /*  If  prep.dataO  was  not  run,  program  looks  for  */ 

/*  the  -CHxx  and  .H  files.  Reads  in  ,H  file.  */ 

{ 

for  (ch=l ;ch<=14;+-*ch) 

{ 


newfilel  =  "MAXIMUMFILENAME . EXTENSION" ; 
nevf ile2  =  "MAXIMUMFILEHAHE . EXTENSION" ; 
strcpy (newf ilel ,xfer) ; 
strcpy(newf ile2,xfer) ; 


if  (ch==l)  8trcat(newfilel,‘  -hi"); 
if  (ch==2)  strcat(nevfilel,".ch2") ; 
if  (ch=»=3)  strcat  (newfilel ch3")  ; 
if  (ch==4)  strcat(newfilel," ,ch4") ; 
if  (ch==6)  strcat (newfilel ch5")  ; 
if  (ch==6)  8trcat(newfilel," ,ch6") ; 
if  (ch==7)  8trcat(newfilel,".ch7") ; 
if  (ch==8)  strcat(newf ilel,".ch8") ; 
if  (ch==9)  strcat (newf ile 1,". ch9" ) ; 
if  (ch==10)  strcat(newfilel,".chlO") 
if  (ch«ll)  8trcat(netffilel,".chll") 
if  (ch*=12)  8trcat(newfilel,".chl2") 
if  (ch==13)  strcat(newf ilel ,".chl3") 
if  (ch==14)  strcat (newfilel ,".chl4") 


if 

{ 


> 


(fopen(newf ilel ,  "r")  »»  HULL) 

piiulf ("\n  Input  file  Xs  does  not  exist  in  this 

directoryW , 

newfilel) ; 

quitO ; 


else 

ofp[ch]=fopen(aewfilel ,  "r"); 


max[ch]*0. ; 
mini [ch]=0. ; 


strcat (newf ile2 , " . h" ) ; 
bfp=fopen(newfile2,  "r"); 

f8eek(hfp,0,0) ; 

for  (ch*l ;ch<*14;++ch) 

facanf (hfp, "XfXfXiXf " ,tmax[ch] ,fcmini[ch] , Asample.rate , 

tstd_dev[ch] ) ; 
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> 


noprep*FALSE ;  /•  sets  flag  to  not  reop«n  files  during  session*/ 

ofp[0] sfopenC'sinetest .dat .h" , "r") ;  /*  opens  data  files  for  sine  */ 

ofp[15]*fopen("sinetest .dat . 1" , "r") ;  /*  reference  signals  */ 

ofpll6]=fopen("sineteet .dat. 2" ,"r") ; 
ofp[l7]=fopen("sinetest .dat . l.Shz" ,"r") ; 
of p [ 18] =f open (" sine test .dat .half hz" ,"r") ; 

for  (n=lS;n<=18;**n) 

fscanf  (ofp[0]  ."XfXfXiV.f"  , fcmar  [n]  ,kmini[n] , Abb. integer  .tstd.dev [n]  )  ; 
total_skip  **  akip_seconds*sample_rate ; 

if  (zero. pointers)  /*  if  starting  3  beginning  of  data*/ 

for  (ch*l ;ch<*14;*+ch) 

f seok(ofp[ch] ,0,0) ; 

k*l;  /*  tracks  •  of  times  thru  main  draw  loop  «/ 

s.diff^O;  /•  tracks  •  of  samples  between  draw  event  •/ 

s_old=0;  /*  sample  •  of  laBt  draw  event  */ 


newfilel*"MAXIMUMFILENAHE.  EXTENSION'1 ;  /•  opens  file  for  diagnostic  •/ 
strcpy (newf ilel ,xfer) ;  /*  output  and  prints  hdr  •/ 

strcat(newf ilel ," .diag") ; 
tfp=fopen(newf ilel ,"v") ; 

fprintf (tfp , "Xs  (reference  ch=Xi)\n\n",JCfer,o) ; 

fprintf (tfp,"  s  SHIFT  buff  data[r]  gray  1-data  l-disp  2-disp 

3-disp  4-disp  5-disp  6-disp  7-disp\n"); 
fprintf (tfp,"  8-disp  9-disp 

10-disp  11-disp  12-disp  13-disp  14-disp\n\n") ; 

buffer=2 ; 


while  ((c  =  gotc(hfp))  !«  EOF)  /*  reads  in  data  fiJe  comments  •/ 

{ 

fscanf (hfp,"XiXiXs" ,Acmt_min[cmt_cnt] ,kcmt.sec [cmt.cnt] ,cmt [cmt.cnt] ) ; 
cmt_sample[cmt_cnt]s( (cmt_min[cmt _cnt] *60) +c*t_soc [cmt.cnt] )*64; 
♦♦cmt.cnt ; 

> 

for  (cmt_cnt=0 ; cmt_cnt<=9 ; ++cmt_cnt) 
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cmt_f lag[cmt_cnt]=FALSE; 


while  ((c  =  getc(ofp[l]))  !*  EOF)  /*  begins  main  loop  for  screen  u/d  */ 

{ 


while  (s  <  total.skip)  /*  bypasses  data  samples  until  */ 

{  /•  selected  start  time  is  reached  */ 

for  (ch*l ; ch<=18 ; ♦♦eta) 

fscanf <ofp[ch] ,"Xf" ,kbb_float) ; 
frontbuffer(TRUE) ; 

if  (s==l) 

{ 

cpack(OxOOOOOOaa) ; 

sprintf (status ,  "SEARCHING"); 

setup_string(status,  SO.fchalf); 

cbov2( .625 , -.62) ; 

fmprstr(status) ; 

cpack (0x777777) ; 

rect ( .61 ,  -  .67 , .87 .74) ; 

cpackOxaa) ; 

rectf ( .62 , - -68 , . 625 .73) j 

> 

if  (s*=(total_akip/8) ) 

rectf(  .62, -.68, .65, -.73); 
if  (s*-(total_skip/4)) 

rectf (  .62,-. 68, .68,-. 73) ; 
if  (a”(total_Bkip/8*3) ) 

rectf ( .62,- .68, .71, -.73); 
if  (s»»(tctal_skip/2) ) 

rectf ( .62 ,-. 68, . 74 , - .73) ; 
if  (s==(total_skip/8*5) ) 

rectf ( .62,-. 6b, .77 ,- .73) ; 
if  (s»"(total_skip/4*3) ) 

rectf ( .62, -.68. .80, -.73) ; 
if  (e**(total_skip/8»7) ) 

rectf ( .62, -.68, ,83, -.73) ; 
if  (s==(total_skip-4) ) 

rectf (.62, -.68, .86,-. 73) ; 
frontbuffer(FALSE) ; 


for  (cmt_cnt=9 ; cmt_cat>=0 ; — cmt.cnt)  /*  chks  if  time  for  new  data  */ 
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/*  comment 


if  (s>=cst_sample[cmt_cnt] ) 

{ 

if  (cmt_f lag(cmt_cnt] ““FALSE) 

{ 

frontbuf fer (TKUE) ; 
cpack(Oxcccccccc) ; 
rectf (-.17, .4, .17, .26) ; 


*/ 


break; 

> 

> 

if  (cmt_flag[cmt_cnt] —FALSE)  /*  draws  data  cmnt  to  screen  */ 

{ 

cpack(OxOOffOOOO) ; 

sprintf (statue,  "Xs",cmt[cmt_cnt]) ; 

setup_string(cmt [cmt.cnt] ,  50,Ahalf); 

cmov2(0.0  -  half ,.3); 

fmprstr (statue) ; 

frontbuf fer (FALSE) ; 


> 


for  (f lag=cmt_cnt ;f lag>«0; — flag) 
cmt_f lag [cmt.cnt] “TRUE ; 


frontbuf fer (TRUE) ;  /‘erases  nsg  area  over  STATUS  box*/ 

cpack(Oxcccccccc) ; 

rectf ( .61 .626 , . 95 , - .525) ; 

frontbuffer(FALSE) ; 


if 

((qtestO)  =*  101  II  qtestO  ==  103) 

/• 

displays  msg  if  rt  or 

*/ 

{ 

/* 

1ft  mouse  pushed  and 

•/ 

frontbuf fer (TRUE) ; 

/• 

stops  draw  loop 

*/ 

cpack(OxOOOOOOaa) ; 

sprintf (status,  "DISPLAY  HALTED"); 

setup_string(status ,  50,fthalf); 

cmov2( . 60,- . 62) ; 

fmprstr (status) ; 

frontbuf fer (FALSE) ; 

break; 

> 

if  (qtest()==102  Aft  Ibatch)  /*  pauses  display  if  middle  mouse  */ 

{ 

qresetO ; 
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f rontbuffer (TRUE) ; 
cpack(OxOOOOOOaa) ; 
sprintf (status,  "DISPLAY  PAUSED"); 
setup_string(status,  50,thalf), 
cnov2( .58, -.62) ; 
foprstr(status) ; 
f rontbuffer (FALSE) ; 

while(qtest() • =102) 

qresetO  ; 

f rontbuffer (TRUE) ; 
cpack(Oxcccccccc) ; 
rectf ( . 57 , - .625 , . 92,- . 525) ; 
frontbuffer(FALSE) ; 


} 

t=s*l .0/sampIe_rate ; 
♦+s; 


/*  time  of  sample  from  start  of  record*/ 
/•  increments  "s"  for  next  pass  */ 


for  (ch*l ; ch<*18 ;**ch)  /•  saves  most  recent  data  values  and  */ 

{  /*  reads  in  new  •/ 

data_oldold[ch]  *  data_old[ch] ; 
data_old[chl  *  data[ch] ; 


fscanf (ofp[ch] ,  "It",  tdata[ch]); 

> 


if  (interval)  /*  enter  if  peak  or  axis-crossing  only  */ 

{ 

if  (peak)  /*  enter  if  peak-driven  drawing  */ 

{ 

if  (  (data_old[b]  >=  data_oldold[b] )  At  (data[b]  <  data_old[b])) 
{  /•  chka  for  two  increases  then  one  drop  */ 


if  ((k  X  interval.skip)  ==  0) 

( 


/*  enter  every  Nth  event  as  */ 
/*  set  by  Draw  Interval  menu  */ 


if  (!wait  At  tbatch)  /*  scales  display  spd  to  #  bins*/ 

for  (ch*l ;ch<=((76-binB)*60000) ;++ch) 
ch"ch ; 

> 


u.diff  *  s-s_old; 
s.old  =*  s; 
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evapbuffersO ; 

♦♦buffer;  /*  diagnostic*  for  .diag  */ 

buffer. now=(buffer  X  2) ; 

erase.barsO ; 

draw.boxesO ; 

draw.statusO ; 

draw.baraO  ; 

} 


♦+k; 


> 

) 

if  (axis. crossing)  /*  enter  if  axis-xings  drive  draw  events  */ 

{ 

if  (  (data_old[b]  <=  0.)  kk  (data[b]  >0.)  II 

(data_old[b]  >0.)  kk  (datafb]  <■  0.)) 

{  /*  chks  for  up  or  downward  crossings  »/ 


if  ((k  X  interval.skip)  »■  0)  /*  enter  every  Nth  event  as  */ 

/•  set  by  Draw  Interval  menu  */ 

{ 

s.diff  ■  s-s_old; 
s.old  *  s; 

if  Ovait  kk  !batch)  /»  scales  display  spd  to  #  bins*/ 

{ 

for  (ch=l ;ch<*( (75-bins ) *60000 ) ;+*ch) 
ch*ch; 

> 

swapbuffersO ; 

♦♦buffer;  /•  diagnostics  for  -diag  */ 

buffer_now»(buffer  X  2) ; 

erase.barsO ; 

draw.boxesO ; 
draw.statusO  ; 

draw.barsO  ; 
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> 


**k; 


> 

> 

} 

if  ((interval)  /*  enter  if  continuous  drawing  •/ 

{ 

swapbuf f ers( ) ; 
erase.barsO  ; 
draw_boxes() ; 
draw_Btatus() ; 
drav_bars() ; 

} 


>  /*  ends  the  main  loop  */ 

while  (!qtest()  left  time_shift) 
swapbuffersO  ; 


> 


/*  ends  the  function  •/ 


/ 


FUNCTION  "DRAW.BARS" 


/ 


void  draw_bars()  /*  used  to  drew  the  bar(s)  in  each  box  */ 

{ 

/•  left  and  right  edges  of  14  channel  boxes  */ 


static  float 

left.edge[16]*{ 

0. , 

-•4, 

-.7,  -.45,  - 

.65, 

-.6, 

-.35, 

-.1, 

.15, 

•  4, 

.35, 

.25, 

■  5, 

2, 

-.1); 

static  float 

right_edg#[15]  ■{ 

0.  . 

-.2, 

-.5,  -.25, 

-.35, 

-.4, 

-.15, 

•  1, 

.35, 

.6, 

.55, 

.45, 

•7, 

•  4, 

-l); 

/*  distance  from  window’s  x-axis  to  center  of  each  channel  box  */ 
static  float  vert.offset[16]={  0.,  .625,  .326,  .175,  -.125, 

-.425,  -.575,  -.725,  -.675,  -.426, 
-.125,  .175,  .325,  .625,  .575); 


linewidth(lw) ;  /*  expands  the  line  segment  to  "lv"  pixels  wide  */ 
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if  (single. buff )  /*  diagnostic  for  .diag  print  */ 

f rontbuffer(TRUE) ; 

for  (ch*l;ch<=18;<"*-ch)  /*  clips  and  scales  bar’s  displace  in  each  ch  •/ 

{ 

if  ( (data_old[ch] >0 . )  tt  (data_old[ch]>(3.*std_dev[ch]))) 
data_old[ch]=(3.*  std_dev[ch] ) ; 

if  ((data_old[ch]<0.)  tt  (data_old[ch]<(-3. *std_dev[ch] ) )) 
data_old[ch]=(-3.*  std_dev[ch] ) ; 

scale. factor [ch]= . 12/ (3 . *  std.devfch] ) ; 

disp[ch]=(  data_old[ch]  *  scale.factor [ch]  ); 

ypos[ch]  =  disp[ch]+vert_offset[ch] ; 


for  (ch=l;ch<=14;**ch)  /*  calcs  start  and  stop  of  each  "bar1',  keeping  */ 
{  /*  track  of  which  bins  they  are  in,  if  time-shftd  */ 


if  (tine. shift) 

{ 

start [ch] [shift] [0]«  (left_edge[ch]*0.0025*((float>shift* 

0 .2/ (f loat)bins) )  ; 

start [ch] [shift] [1]»  ypoe[ch]; 
start [ch] [shift] [2]"  0.; 

stop[ch] [shift] [0]-  (start[ch] [shift] [0] -0 . 0026+C0 . 2/ (f loat)bins) ) ; 
stop[ch] [shift] Cl] *  ypostch]; 
stop [ch] [shift]  [2]*  0.; 

> 


if  Otime.shift) 

{ 

start_old[ch] [shift] [0]=start[cb] [shift] [0] ;  /*  xfers  new  to  old  */ 

start_old[ch] [shift] [1] «start [ch] [shift] [1] ;  /*  (old  needed  to  erase)*/ 

start_old[ch] [shift] [2] “start [ch] [shift] [2] ; 

stop_old[ch] [shift] [0]=stop[ch] [shift] [0] ; 

stop_old[ch] [shift]  [l]*stop[eh] [shift]  [1] ; 

stop_old[ch] [shift] [2]*stop[ch] [shift] [2] ; 

start [ch] [shift] [C]»  (left .edge  [ch]+ .005) ; 
start[ch] [shift] [1]=  ypos[ch]; 
start [ch] [shift] [2]*  0.; 

stop[ch]  [shift]  [0]  =  (right .. edge  [ch]  +  .  005) ; 
stop[ch] [shift] [1]*  ypos[ch]; 
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stop[ch] [shift]  [2]  =  0.; 

} 

if  ( !ccnstant_intensity  kb  !synch_conatrained)  /*  set  gradient  bar  color*/ 

< 

packd*(data_old[b]*(  107 .  /(3 .  •  std.devfb] )  )-*-147. ) ; 

pack=packd; 
packcolor  [ch]*»  pack; 

packcolor [ch] =  (packcolor [ch]  <<  3)  I  pack; 
packcolor [ch] =  (packcolorfch]  <<  8)  I  pack; 
packcolor [ch]*  (packcolor [ch]  «  8)  I  pack; 

/«  creates  a  "packed  integer"  conaiating  of  4  8-bit  */ 
/•  integers  which  are  the  truncated  results  of  */ 

/*  scaling  the  EEG  Ch  b  voltage  to  a  255-step  gray  */ 

/*  scale.  All  4  integers  are  the  sane  value  to  •/ 

/*  create  the  gray  scale  by  assigning  identical  vala*/ 

/*  to  the  alpha,  red,  green,  and  blue  components  */ 

cpack(packcolor [ch] ) ;  /*  seta  the  draw  color  */ 


if  (synch.constrained)  /•  aet  thresholded  bar  color  */ 

{ 

if  ((diaptb]  <»  synch. linit)  kk  (diap[b]  >«  -synch.lioit) ) 
cpack(OxO) ; 

else 

cpack(Oxaaaaaaaa)  ; 


if  (constant.intensity)  /*  set  constant  bar  color  */ 

cpack(OxO) ; 


bgnlineO;  /•  draws  the  bar  using  the  arrays  "startl"  and  "atopl"*/ 

v3f (start  [ch]  [shift] ) ; 
v3f (atop[ch] [shift] ) ; 
endline  0 ; 

if  (tine.shift)  /•  draws  cursor  if  display  is  tine-shifted  */ 

{ 

’  cursor_top_old[ch]  [0]  =cursor_top[cli]  [0]  ; 
cursor. top_old[ch] [1] *cursor_top[ch]  [l] ; 
curacr_top.old(ch] [2] *cursor.top[ch] [2] ; 
cursor„bottom_old[ch]  [0]*cursor_botton[cb] [0] ; 
curaor_bottom_old[ch]  [l]*cursor.botton[ch] [1] ; 
cursor_botton_old[ch]  [2]*cursor_botton[ch] [2] ; 
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cureor_top[ch] [0]*atop[ch] [shift] [0] ; 
cursor .top [ch] [1] =*  vert  .of  f  set [ch] ♦ . 125 ; 
cursor.top [ch] [2] *0 . ; 
cursor_bottoa[ch] [0]=cursor.top[ch] [0] ; 
cur sor .bottom  [ch]  [1] -vert .of fast [ch] + . 11 ; 
curaor_bottoa[ch] [2]*0. ; 

cpack (OxOOf f 0000) ; 
linewidthO ) ; 

bgnlineO  ; 

v3f (curaor_top[ch] ) ; 
v3f ( cur sor.bottom [ch] ) ; 
endline () ; 

linewidth(lv) ;  /*reset  linewidth  to  bar  value*/ 


> 


/*  prints  .diag  file  */ 

fprintf (tfp,"X4i  X2i  Xi  X7.3f  X3i  X7.3f  X7.4f  X7.4f  X7.4f  X7.4f 

X7.4f  X7.4f  X7.4f  Xi:Xf\n" , (s-1) .shift ,buffer_nov,data_old[b] .pack, 

data.oldtl] ,disp[l] ,disp[2] ,disp[3] ,disp[4] ,disp[5] ,disp[6] ,diap[7] , 
( (s-D/64/60) , ( ( (((float )s-l . )/64 . /60 . )- ( (b-1)/64/60) )*60 . ) ) ; 
fprintf  (tfp, "  X7.4f  V  At  V  At 

X7.4f  X7.4f  X7.4f  X7 .4f\a\n'‘ ,disp[8]  ,disp[9]  , 

diap[10] , dispCll] ,diap[12] ,diap[13]  ,disp[l4]) ; 


if  (interval  tt  bina<»24)  /*  draws  correllation  grid  */ 

drau.matrixO ; 

if  (time.shift) 

{ 

if  (shift  <*  (bins-2))  /*  increments  shift  #  •/ 

♦♦shift ; 

else 

{ 

shift  ■  0; 

if  (batch  Jkt  qteat()s=102)  /*  handles  middle  mouse  in  */ 

qresetQ;  /*  step  and  batch  modes  */ 

while  (batch  tt  qtestO  !*102) 

{ 

if  (qtest()«101  II  qtest()*“i03) 
break ; 

> 
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if  (batch  kk  qtest()==102) 
qreeetO ; 


> 

> 


if  (wait)  /*  purges  unwanted  mouse  events  */ 

{ 

while  (qtestO  !=101  kk  qteatO  •  =  102) 

( 

if  (qtest()==103) 
qreaetO ; 

> 

if  (qtestO  !*101) 
qresetC ) ; 

> 

>  /*  ends  the  function  */ 


. . . . 

void  erase_bars() 

{ 

if  (single.buff) 

frontbuffer(TRUE) ; 


FUNCTION  "ERASE.BARS"  *« 

/•  used  to  erase  the 
/*  have  been  noved 

/*  diagnostic  for  .diag 


•••,••••**««*•*.*********/ 

bars  after  they  */ 
to  the  back  buffer  */ 

print  •/ 


cpack(Oxcccccccc) ;  /*  sets  drawing  color  equal  to  backgrnd  •/ 

linevidth(lw) ;  /*  expands  the  line  segment  into  a  lv-pixel-wide  bar  */ 

if  (tine.shift) 

{ 


for  (ch=l  ;ch<=14;*+ch) 

{ 

bgnlineO;  /*  overwrites  bar  #ch  in  backgrnd  color  •/ 

v3f (start[cb] [(shift)] ) ; 
v3f (stop[ch] [(shift)] ) ; 
endline () ; 


frontbuf far (TRUE) ; 
linevidth(l) ; 

if  (! single.buff)  /*  overwrites  last  cursor  if  time-shifted  •/ 
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bgnlineO ; 

v3f (cursor_top_old[ch] ) ; 
v3f (cursor_bottom_old[ch] ) ; 
endlineO ; 

> 

else 

{ 

bgnlineO  ; 

v3f (cursor_top[ch] ) ; 
v3f (cursor_bottom[ch] ) ; 
endline () ; 

> 

linewidthdu )  ; 
if  ( ! single.buff ) 
frontbuffer (FALSE) ; 


> 

> 

if  ( !time. shift) 

{ 

for  (eh=l ;ch<=14;++ch) 

i 

bgnlineO;  /*  overwrites  bar  ich  in  backg  rod  color  */ 

v3f (start. old [ch] [shift] ) ;  v3f (stop_old[ch] [shift] ) ; 

endline () ; 

> 

} 

>  /*  ends  the  function  */ 

. . .  FUNCTION  "DRAW.HATRIX"  . . . . . 

void  draw.matrix(void)  /*  draws  correlation  grid  fi  screen’s  center  */ 

{ 

static  float  xedge[6]={  -.15,  -.09,  -.03,  .03,  .09,  .15); 

static  float  yedge[6]a{  .025,  -.05,  -.125,  -.2,  -.276,  -.35}; 

static  float  xoff».004; 

static  float  yoff*.005; 

float  xlshift*0.; 

float  x2ahift=0.; 

frontbuffer (TRUE) ; 


linewidth(l) ; 

cpack(0x00ff0000) ; 

rect(xedge[l] ,yedge[0] ,xedge[2] ,yedge[l]) 
rect(xedge [2] ,yedge[0] ,xedge[3] ,yedge[l]) 
rect(xedge[3] ,yedge[0] ,xedge[4] ,yedge[l]) 
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rect (xedge [0] ,yedge[l] ,xedge[l] ,yedge[2]) ; 
rect (xedge [1] .yedge  [1] .xedge [2] .yedge [2]) ; 
rect (xedge [3] .yedge  [1] .xedge [4] .yedge [2]) ; 
rect(xedge[4]  ,yedge[l]  ,xedge[5?.  ,yedge[2] ) ; 
rect  (xedge  [1]  , yedge  [2]  , xedge  12]  , yedge  [3] )  ; 
rect (xedge [3] .yedge  [2] , xedge [4] ,yedge[3] ) ; 
rect  (xedge  [0]  .yedge  [3]  , xedge  [1]  .yedge  [4] )  ; 
rect  (xedge  [1]  .yedge  [3]  .xedge  [2]  .yedge  [4]) ; 
rect  (xedge  [3]  .yedge  [3]  .xedge  [4]  .yedge  [4]) ; 
rect  (xedge  [4]  .yedge  [3]  .xedge  [5]  .yedge  [4]) ; 
rect (xedge [2] ,yedge[4] ,xedge[3] ,yedge[5] ) ; 

xlshift-shift* ( .052/(f loat)bins) ; 
x2shift=(shift+l)*( .052/(float)bine) ; 

tor  (ch=l ;ch<=14;*+ch) 

correl[ch]=fabs(diep[ch]-di8p[b] ) ; 

if  (correl[ch]<=( . 12*. 2)) 
cpack(OxO) ; 

if  (correl[ch]>(.12*.2)  tk 

correl[ch] <=( (3. *atd_dev[ch] )* .3) ) 
cpack(0x00777777) ; 

if  (correl[ch]>(.12*.3)) 
cpack(OxOOeeeeee) ; 

if  (ch— 1) 

rectf ( (xedge [l]+xoff ♦xlahift) , (yedge [0]-yoff ) , 
(xedge[l]+xoft+x2shift) , (yedge [l]+yoff)) ; 
if  (ch— 14) 

rectf ((xedge [2] ♦xoff ♦xlahift) , (yedge [0]-yoff) , 
(xedge [2] ♦xoff +x2shift) , (yedge [l]*yoff)) ; 
if  (ch==13) 

rectf ((xedge [3] +xoff+xlahift) , (yedge [0]-yoff) , 
(xedga[3]+xoff+x2shift) , (yodge[l]+yoff)) ; 
if  (ch— 2) 

rectf ( (xedge [0] +xoff +xlahif t) , (yedge [1] -yoff ) , 
(xedge [0]+xoff+x2ahift) , (yedge [2]+yoff)) ; 
if  (ch— 3) 

rectf ((xedge [1] +xoff ♦xlahift) , (yedge [1] -yoff ) , 
(xedge [l]+xoff+x2ahift) , (yedge [2] ♦yoff ) ) ; 
if  (ch— 11) 

rectf ((xedge [3] +xoff+xlahift) , (yedge [1] -yoff) , 
(xedge [3]«-xoff*x2Bhift)  .(yedge [2] ♦yoff)) ; 
if  (ch— 12) 

rectf ((xedge [4] ♦xoff ♦xlahift) , (yedge [1] -yoff) , 
(xedge [4] ♦xof f ♦x2shift) , (yedge [2] +yof f ) ) ; 
if  (ch— 4) 
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rectf ((xedge[l]+xoff+xlahift) , (yedge [2] -yoff ) , 
(xedge [1]  +xoff+x2shif t) ,  (yedge  [3]  +yoff ) )  ; 
if  (ch*»10) 

rectf ((xedge[3]+xoff+xlshift) , (yedge [2] -yoff ) , 
(xedge [3] +xof f+x2ahif t) , (yedge [3J  *yoff ) ) ; 
if  (c.V==5) 

rectf ((xedgerjj+xoff+xlahift) , (yedge [3] -yoff ) , 
(xedge 'd] +xoff+x2shift) , (yedge[4]+yoff)) ; 
if  (ch==6) 

rectf ((xedge  Ll]+xof f+xlshift) , (yedge [3] -yoff) , 
(xedge [1] +xoff*x2shif t) , (yedge [4] +yof f ) ) ; 
if  (ch-=8) 

rectf((xedge[3]+xoff+xlshift) , (yedge [3]-yoff) , 
(xedge[3]+xoff+x2ahift) , (yedge[4]+yoff)) ; 
if  (ch=-9) 

rectf ((xedge [4] +xof f+xlshift) , (yedge [3] -yoff) , 
(xedge [4] +xoff +x2shift) , (yedge[4]+yoff)) ; 
if  (ch*=7) 

rectf ((xedge[2]+xoff+xlahift) , (yedge [4] -yoff ) , 
(xedge[2]+xoff+x2shift) , (yedge [5] +yoff ) ) ; 


> 


frontbuffer (FALSE) ; 

> 


. . .  FUNCTION  "DRAW.BOXES"  **+••*•••**•••**** . . 

void  draw_boxea(void)  /*  draws  boxes  on  screen  for  ch  displays  */ 

{ 

linewidth(l) ; 

cpack(OxOOffOOOO) ;  /*  aeta  line  color  to  blue  */ 

rect(-.4, .75,-. 2, .625) ;  /*  CH  il  •/ 

rect(-.4, .625, -.2, .5) ; 

sprintf (channel ,  "CH  1"); 
setup. stting(channel ,  60,*half); 
caov2(-.39, .51) ; 
faprstr( channel) ; 

rect(-.7, .45, -.5, .325) ;  /*  CH  #2  */ 

rect(- .7 , .325, -.5, .2) ; 
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aprintf (channel ,  "CH  2"); 
setup. string(channel ,  60,Jthalf); 
cmov2(- . 69 , .21) ; 
fmprstr (channel) ; 

rect (- . 45 , .3,- .25 , . 175) ;  /•  CH  #3  */ 

rect(-.4S, .175, -.25,  .05) ; 

8printf (channel ,  "CH  3"); 
setup. string(channel ,  60,4half); 
cmov2(- . 44 , .06) : 
fmprstr (channel) ; 

rect (-.55,0. ,-.35,-. 125) ;  /*  CH  *4  */ 

rect (-.55, -.125, -.35, -.25) ; 

sprintf (channel ,  "CH  4"); 
setup_string(channel ,  60,thalf); 
cmov2(-.54,-.24) ; 
fmprstr(channel) ; 

rect(- .6 .- . 3,- .4 ,- . 425) ;  /•  CH  #6  */ 

rect (-.6, -.425, -.4, -.55) ; 

sprintf (channel,  "CH  5"); 
setup_string(channel ,  60,khalf); 
cmov2(-.59,-.54) ; 
fmprstr (channel) ; 

rect (- . 35 , - .45 ,- . 15 , - .575) ;  /•  CH  #6  */ 

rect (- . 35 , - .575 , - . 15 , -.7) ; 

sprintf (channel ,  "CH  6"); 
setup. string (channel ,  60,thalf); 
cmov2(-.34,-.69) ; 
fmprstr (channel) ; 

rect (-. 1 ,- .6, . 1 .725) ;  /*  CH  #7  */ 

rect (-. 1 ,- .725 , . 1 ,- . 85) ; 

sprintf (channel ,  "CH  7"); 
setup. string(channel ,  60,thalf); 
cmov2(-.09,-.84) ; 
fmprstr (channel) ; 

rect ( . 15 ,- . 45, .35,- . 575) ;  /•  CH  #8  */ 

rect( . 15,- .575, .35 ,- .7) ; 

sprintf (channel ,  "CH  8"); 
setup_string( channel ,  60,thalf); 
caov2(. 16, -.69) ; 
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fmprstr( channel) ; 


rect (.4, -.3, .6,- .425) ;  /*  CH  #9  */ 

rect (.4, -.425, .6, -.55) ; 

sprintf (channel ,  "CH  9"); 
eetup_8tring(channel ,  60,thalf); 
cmov2( .41 54) ; 
fmprstr (channel) ; 

rect (. 35,0 ... 55 125) ;  /*  CH  #10  */ 

rect( .35,- . 125, .55,-. 25); 

sprintf (channel ,  "CH  10"); 
setup_string(channel ,  60,4half); 
cmov2( .36 ,-.24) ; 
fmprstr(channel) ; 

rect (.25, .3,-45, .175) ;  /*  CH  #11  */ 

rect( .25, .175, .45, .05) ; 

sprint! (channel,  "CH  11"); 
setup_8tring (channel ,  GO.fchalf); 
cmov2( .26, .06) ; 
fmprstr (channel) ; 

rect (.5, .45, .7. .325) ;  /*  CH  #12  •/ 

rect( .5, .325, .7, .2) ; 

sprintf (channel,  "CH  12"); 
setup_string(channel,  60,#half); 
cmov2(  .51 , . 21) ; 
fmprstr (channel) ; 

rect(.2,.75, .4, .625) ;  /*  CH  #13  */ 

rect( .2, .625, .4, .5) ; 

sprintf (channel ,  "CH  13"); 
setup_string(channel,  60,thalf); 
cmov2(  .21 , .51) ; 
fmprstr (channel) ; 

rect(-.l, .7, .1..575);  /*  CH  #14  */ 

rect(-.l, .575, .1,,45) ; 

sprintf (channel ,  "CH  14"); 
setup. string(channel ,  60,fchalf); 
caov2(-09, .46) ; 
fmprstr (channel) ; 


> 


/*  end  of  function  */ 
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. . MM  FUNCTION  "DRAW.STATUS"  ..***.**•***.*•*** 

void  drau_statu8(void) 

{ 

cpack(OxOOeeeeee) ;  /•  draws  Status  box  in  near-white  */ 

rectf ( . 55 , - . 65 , . 95 , - . 95) ; 

cpack(OxOOeeeeee) ;  /*  draws  Configuration  box  in  same  */ 

rectf (-.96,- .65, -.55, -.95) ; 


time.printO  ; 
cpack(OxOOaaOOOO) ; 


/•  draws  elapsed  time  «/ 

/*  sets  color  to  dk  red  */ 


/ 


/*  remainder  draws  information  into  display  boxes  */ 

sprintf (status ,  "STATUS"); 
setup_3t.nng(8tatus ,  SO.Jthalf); 
cmov2( .68 , - . 70) ; 
fmprstr(status) ; 

draw_f  ilenameO  ; 
draw.updateO ; 
draw_freq() ; 

cpack(OxOOaaOOOO) ; 

aprintf (status,  "CONFIGURATION") ; 
setup. string(8tatus ,  50,thalf); 
cnov2(-.925,~.70) ; 
fmprstr(status) ; 

cpack(OxOOaaOOOO) ; 


if 

(b<*14) 

sprintf (statue , 

"Reference 

Ch: 

Xi",  b); 

if 

(b==l5) 

sprintf (status , 

"Reference 

Ch: 

1-Hz  SINE"); 

if 

(b==16) 

sprintf (status , 

"Reference 

Ch: 

2-Hz  SINE"); 

if 

(b==l7) 

sprintf (status , 

"Reference 

Ch: 

1.5-Hz  SINE"); 

if 

(b==18) 

sprintf (status , 

"Reference 

Ch: 

0.5-Hz  SINE") ; 

setup. 8tring(status,  70,fchalf) 
cmov2(-.93,-.75) ; 
fmprstr(status) ; 
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if  ((interval)  tt  (peak)) 

sprintf (status,  "Draw  Event :  Peak"); 
if  ((interval)  kk  (axis.crossing)) 

sprintf (status,  "Draw  Event:  Axis  Xing"); 

if  (! interval) 

sprintf (status,  "Draw  Event:  Constant"); 
set up. st ring (status,  70,fchalf); 
cmov2(- .93 , - .80) ; 
fmprstr(status) ; 

if  (interval) 

sprintf  (status,  "  ->  Display  Interval:  */,i"  .interval  sk)p) 

setup_8tring(status ,  70,4half); 
cmov2(-.93,-.85)  ; 
fmprstr(status) ; 

if  (wait) 

{ 

cpack(Cxaa) ; 

sprintf (status,  "STEP-DRAW"); 
setup_string( status ,  50,khalf); 
cmov2(- . 86 , - . 58) ; 
fmprstr (status) ; 

sprintf (status,  "(Middle  Mouse  *  NEXT)"); 
setup.st  ring  (status,  70,lthalf); 
cmov2(-.89,-.62) ; 
fmprstr (status) ; 

cpack(OxOOffOOOO) ; 

> 

if  (batch) 

{ 

cpack(Oxaa) ; 

sprintf (status,  "BATCH-DRAW"); 
setup..string(status,  SO.thalf); 
cmov2(-.89,-,58) ; 
fmprstr(status) ; 

sprintf (status,  "(Middle  Mouse  =  NEXT)"); 
setup. string(status,  70,Ahalf); 
cmov2(-,91  ,-.62) ; 

ffflprstr(status) ; 

cpack(OxOOffOOOO) ; 

> 

if  (constant.intensity) 
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sprintf (status,  "Ref  Ch  Freq:  X4.2f  Hz",  freq) ; 

else 

sprintf (status,  "Ref  Ch  Freq:  N/A"); 

setup„string(status ,  70,Ahalf); 
cmov2( .58 ,- .90) ; 
fmprstr (status) ; 


> 

. . ****.....*****  FUNCTION  "DRAV.TITLE"  *.**•*•.**•***«*»*****•***/ 

void  drav.title (void)  /*  draws  title,  electrodes,  etc  •/ 

{ 

frontbuff er (TRUE) ; 

ortho2(-1.0,  1.0,  -1.0,  1.0); 
depthcue(FALSE) ; 

cpack (0x00777777) ; 

sprintf (title,  "TOPOS") ; 
setup_string(title ,  15,Ahalf); 
cnov2(- .925 , .805) ; 
fmprstr (title) ; 

cpack(Oxaa) ; 
cmov2(-.93, .80) ; 
fmprstr(title) ; 


cpack(0x77) ; 

sprintf (title ,  "The  AFIT"); 
setup_string(title ,  40,fchalf); 
cnov2(-.93,.70); 
fmprstr (title) ; 

sprintf (title ,  "Interactive"); 
e#tup_string(titl« ,  40,thalf); 
c»ov2(- .93, .63) ; 
fmprr,r(title) ; 

sprintf (title,  "Toposcope") ; 
setup. string(title ,  40,Ahalf); 
cmov2(- . 93 , . 56) ; 
fmprstr(title) ; 

cpack(0x777777); 
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sprintf (title ,  "Face"); 
aetup_atring(title ,  45,fchalf); 
cnov2 (0.0-half, .95)  ; 
fBpratr(title) ; 

aprintf (title,  "Right"); 
aetup_string(title ,  45,fchalf); 
cmov2(  .86,0.;; 
fnpratr (title) ; 

aprintf (title ,  “Ear"); 
aetup_atring(title ,  45,4half); 
cbov2(.88,-.07); 
fnpratr(title) ; 

aprintf (title ,  "Left"); 
83tup_atring(title ,  45,fchalf); 
cmov2(-.95,0. ) ; 
fBpratr(title) ; 

sprintf (title ,  "Ear"); 
aetup_atring(title ,  45,thalf); 
cbov2(-.945,-.07)  ; 
fnprBtr(title) ; 


circf (-.15, .625, .01) ; 

/*  electrode:  f3  */ 

circf (-.5, .6, .01) ; 

/*  electrode:  f7  */ 

circf <>.76, .125, .01); 

/*  electrode:  t3  */ 

circf (-.15, . 125, .01) ; 

/*  electrode:  c3  */ 

circf (-.7, -.35, .01)  ; 

/*  electrode:  t5  */ 

circf (-.25, -.35. .01) ; 

/*  electrode:  p3  */ 

circf (-.25, -.775, .01) ; 

/•  electrode:  ol  */ 

circf (  .25,- .775, .01) ; 

/*  electrode:  o2  */ 

circf ( .25,- .35, .01)  ; 

/*  electrode:  p4  */ 

circf ( .7, -.35,  .01) ; 

/*  electrode:  t6  */ 

circf ( . 15, . 125, .01) ; 

/•  electrode:  c4  */ 

circf ( .75, . 125, .01) ; 

/*  electrode:  t4  •/ 

cir ;f ( .6, .6, .01) ; 

/*  electrode:  f8  */ 

circf ( . 15, .625, .01)  ; 

/*  electrode:  f4  */ 

aprintf (electrode ,  "F3"); 
eetup_atring(electrode ,  60,thalf) ; 
cbov2(-. 165, .645) ; 
fnpr8tr(electrode) ; 

Bprintf (electrode ,  "F7"); 
aetup_8tring(electrode ,  60,thalf ) ; 
cbov2(-.515, .62) ; 
fBpratr(electrode) ; 
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sprintf (electrode ,  "T3"); 
aetup_string(electrode ,  60,fchalf) ; 
cmov2(- .765, .075) ; 
fmprstr (electrode) ; 

sprintf  (electrode ,  ''C3"); 
setup_string(electrodo ,  50 .fehalf ) ; 
cmov2(- . 165 , . 145) ; 
fmprstr (electrode) ; 

sprintf  (electrode  ,  *'TS  “ ) ; 
setup_string(electrodn ,  60,fchalf ) ; 
cmov2(-.715,-.4) ; 
fmprstr (electrode) ; 

sprintf (electrode ,  "P3" ) ; 
setup_string(electrode ,  60,fchalf) ; 
cmov2(- . 265 ,- .33) ; 
fmprstr (electrode) ; 

sprintf (electrode ,  "01"); 
8«tup_string(electrode ,  60 , fchalf ) ; 
cmov2(- . 265 .825) ; 
fmprstr(electrode) ; 

sprintf (electrode ,  "02"); 
setup_string(electrode ,  60,*half ) ; 
cmov2( . 235 .825) ; 
fmprstr (electrode) ; 

sprintf (electrode ,  "P4"); 
setup_string(electrode ,  60 .fchalf ) ; 
cmov2( .235, -.33) ; 
fmprstr(electrode) ; 

sprintf (electrode ,  "T6"); 
setup„string(electrode ,  60,thalf) ; 
cmov2( .686, -.4) ; 
fmprstr(electrode) ; 

sprintf (electrode ,  "C4"); 
setup_string(electrode  ,  60  .Itltalf ) ; 
cmov2( . 135, . 145) ; 
fmprstr (electrode) ; 

sprintf (electrode ,  "T4"); 
setup_string(electrode ,  60,fthalf) ; 
cnov2( .735, .075) ; 
fmprstr(electrode) ; 

sprintf (electrode ,  "F8"); 
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setup_string(electrode ,  60 ,ihalf ) ; 
cmov2( .485, .62) ; 
faprstr(electrode) ; 

sprintf (electrode ,  "F4"); 
setup.string(electrode ,  60 , Jthali ) ; 
cmov2( . 136, .645) ; 
fmprstr (electrode) ; 

linewidth(2) ; 


arc (0 . , .86. .08,-100,1900) ; 
arc( .25 , .77, .06,200, 1600) ; 
arc (-.25,  .77, .06,200,1600) ; 

deflinestyled  ,0xcccc) ; 
aetlinestyle (DOTTED) ; 

ptl  [0] »  -.1; 
ptl[l]  =  .626; 
pt2[0]  =  -.2; 
pt2[l]«  .625; 
draw.lineO ; 

ptl[0]«  -.4; 
ptl  [13 *  .625; 
pt2[0]=  -.5; 
pt2[l]  =  .6; 
draw.lineO ; 

ptl[0]»  -.5; 

ptl  Cl]  ■  -6; 

pt2[0]»  -.6; 

pt2[l]»  .45; 
draw.lineO ; 

ptl[0]=  -.7; 
ptl[l]-  .25; 
pt2[0]*  -.75; 
pt2  [1]  =  .125; 
draw.lineO ; 

ptl[0]x  -.75; 
ptl[l]=  .125; 
pt2  [0] =  -.45; 
pt2  [1] =  .125; 
draw.lineO  ; 


pt  1  [0]  ■  -.25 
ptl[l]  =  .125 
pt2[0]»  -.15 
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pt2[l]  = 

.125; 

draw.lineO ; 

ptl[0]* 

-.15; 

ptl [1]= 

.125; 

pt2[0]  = 

-.35; 

pt2[l]  = 

-.05; 

draw.lineO ; 

ptl [0]= 

-.55; 

ptl[l]  = 

-.2; 

pt2[0]  = 

-■7; 

pt2  [1]  = 

-.35; 

draw.lineO  ; 

ptl [0]= 

-.7; 

ptl[l]  = 

-.35; 

pt2  [0]  = 

-.6; 

pt2[l]  = 

-.35; 

draw.lineO  ; 

ptl[0]  = 

-.4; 

ptl Cl]- 

-.36; 

pt2[0]=> 

-.25; 

pt2 [1]* 

-.35; 

draw.lineO  ; 

ptl [0]* 

-.26; 

ptl[l]  = 

-.35; 

pt2[0]» 

-.26; 

pt2  [1]  = 

-.45; 

draw.lineO ; 

ptl [0]= 

-.25; 

ptl  Cl]  - 

-.7; 

pt2C0]« 

-.25; 

pt2[l]  = 

-.775 

draw.lineO  ; 

pti[o]  = 

-.25; 

ptl[l]» 

-.775 

pt2[0]« 

-.1; 

pt2[l]  = 

-.775 

draw.lineO  ; 

ptl [0]* 

.  l ; 

ptl[l]» 

-.775 

pt2[0]» 

.25; 

pt2[l]= 

-.775 

draw.lineO  ; 

ptl[0]«  .25; 
ptl[l]  =  -.775 
pt2[0]»  .25; 
pC2 [1] “  -.7; 
drau.lineO  ; 

ptl  [0]  =  .25; 
ptltl]-  -.45; 
pt2[0]  =  .25; 
pt2[l]=  -.35; 
drau.lineO  ; 

ptl[0]  =  .25; 
ptl[l]  =  -.35; 
pt2f0]»  .4; 
pt2[l]  =  -.35; 
drau.lineO  ; 

ptl[0]=  .6; 
ptl[l]  =  -.35; 
pt2[0)=  .7; 
pt2[l]*  -.35; 
drau.lineO ; 

ptl  [0]=  .7; 
ptl[l]  =  -.35; 
pt2  [0]  -  .55; 
pt2[l]=  -.2; 
drau.lineO  ; 

ptl  [0]  =  .35; 
ptl[l]=  -.06; 
pt2[0]=  .15; 
pt2[l]3  .125; 
drau.lineO  ; 

ptl[0]=  .15; 
pti[l]=  .126; 
pt2[0]s  .25; 
pt2[l]=  .125; 
drau.lineO ; 

ptl[0]  =  .45; 
ptl[l]«  .125; 
pt2[0]=  .75; 
pt2[l]=  ,126; 
drau.lineO ; 

ptl  [0] =  .76; 
ptl[l]=  .125; 
pt2[0]“  .7; 


pt2[l]=  .25; 
drav.lineO ; 

ptl[0]=  .6; 
ptl[l]=  .45; 
pt2[0]=  .5; 
pt2[l]=  .6; 
drav.lineO ; 

ptl [0] =  .5; 
pt 1 [ 1 ] =  . 6 ; 
pt2[0]*  .4; 
pt2[l]=  .625; 
drav.lineO  ; 

ptl[0]=  .2; 
ptl[l]=  .625; 

pt2[0]=  .1; 

pt2[l]=  .625; 
drav.lineO  ; 

setlinestyle (SOLID) ; 

frontbuf fer(FALSE) ; 


/•***.•••**•****»*******••****  FUNCTION  "DRAW. LINE"  ***♦***•***•••**•*•***•***/ 

void  drav_line(void)  /*  used  immediately  above  by  draw.titleO  */ 

{ 

bgnlineC ) ; 
v3f (ptl) ; 
v3f (pt2) ; 
endlineO  ; 
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1.3  User  Menus:  tmenu.c 


/• 

* 

*  tmenu.c 

* 

*  Part  of  "TOPOS"  by  Capt  Dwight  Roblyer',  AFIT/GS0-92D 

*  November  1992 

* 

*  Based  on  "CURVE  DEMO" 

*  by  Howard  Look  for  Silicon  Graphics,  June  1989 

* 

*  This  file  contains  the  routines  used  for  initializing  and 

*  maintaining  the  menus  that  make  up  the  profile  program. 

* 

*/ 


•include  <Btdio.h> 

•include  <gl.h> 

•include  <device.h> 

•include  <math.h> 

•include  "tevent.h" 

•include  "topos.h" 

/•  Prototypes  */ 

void  init.menus(void) ;  /*  initializes  menu  */ 

void  do_menus(void) ;  /*  calls  menus  */ 

void  remake_curve_menu(void) :  /*  builds  a  menu  */ 

void  remake.display .menu (void) ; 

void  remake. interval. menu(void) ; 

void  remake. bin. menu(void) ; 

void  remake.start.menuCvoid) ; 

void  remake. timer.menu(void) ; 

void  remake. synch.menu(void) ; 

void  remake.line.menu(void) ; 

void  remake. intensity.monu(void) ; 

void  remake. all.menus(void) ; 

extern  void  prep.data(void) ; 
extern  void  set.b(int); 
extern  void  set.lw(int) ; 
extern  void  set.intensity(int) ; 
extern  void  set.bins(int) ; 
extern  void  eet.defaultn(void) ; 
extern  void  change.display(int) ; 
extern  void  set.interval.skip(int) ; 
extern  void  leave.ptrs(void) ; 
extern  void  zero.ptrs(void) ; 
extern  void  position. ptrs(int) ; 


147 


extern  void  save.quit(void) ; 
extern  void  quit(void); 
extern  void  invert.order(void) ; 
extern  void  help(void) ; 
extern  void  end_help(void) ; 
extern  void  svap(void) ; 
extern  void  set_teat(int) ; 
extern  void  drav_display(void) ; 

extern  char  xfer[] ; 

/*  Declare  menus  */ 

static  long  curve.menu;  /*  Top  level  Menu  */ 

static  long  display.menu; 

static  long  interval. menu; 

static  long  bin.menu; 

static  long  start. menu; 

static  long  timer.menu; 

static  long  synch. menu; 

static  long  line.menu; 

static  long  intensity.menu; 

static  long  test.menu; 


/************************  FUNCTION;  INIT.MENUSO  **.*************************/ 

/*  Allocate  and  initialize  the  menus,  then  tell  the  event  manager 
*  to  watch  for  right  mouse  button  clicks. 

*/ 

void  init_menus(void) 

{ 

curve.menu  *  nevpupO; 
display.menu  =  newpupC); 
interval.menu  *  nevpupO; 
bin.menu  =  nevpupO; 
start. menu  =  newpupO; 
timer.menu  =  nevpupO; 
scale. menu  =  nevpupO; 
synch.menu  =  nevpupO; 
line.menu  =  nevpupO; 
intensity.menu  =  nevpupO; 
test.menu  =  nevpupO; 
aet.initial.conditionsO ; 
remake.all.menusO ; 

add. event (ANY,  RIGHTHOUSE,  DOWN,  do.menus,  NULL); 
qde vice (RIGHTHOUSE) ; 

> 

/*•****»*• . * . .  FUNCTION;  DO.MENUSO  *********** . . . **/ 

/« 
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•  Called  by  the  event  manager  whenever  the  right  mouse  button  is 

•  pressed . 

*/ 

void  do.menus(void) 

{ 

end_help( ) ; 
dopup(curve.menu) ; 

> 


. . . . .  FUNCTION:  REMAKE.CURVE.MENIK)  ***#*.**...*.*».****../ 

void  remake_curve_menu(void) 

{ 

f reepup(curve_menu) ; 
curve. menu  =  newpupO  ; 

addtopup(curve_menu,  "TOPOS  */,t");  /•  title  */ 
addtopup(curve_menu,  "Prepare  Data  'Li",  prep.data) ; 
addtopup (curve. menu.  "Reference  CB  Xm",  synch.menu) ; 
addtopup(curve.aenu,  "Bar  Width  Xm",  line.menu) ; 
addtopup (curve .menu,  "Bar  Intensity  Xm",  intensity.menu) ; 
addtopup (curve. menu,  "Display  Mode  Xm",  display.menu); 
addtopup (curve .menu,  "  Draw  Interval  Xm",  interval.menu) ; 
addtopup(curve_aenu,  "  Bin  Number  Xm",  bin. menu) ; 
addtopup (curve.menu,  "Start  Display  Xm",  start.menu) ; 
addtopup (curve.menu,  "Help  Xf " .  help); 
addtopup (curve.menu,  "Test  Xm" ,  test .menu) ; 
addtopup (curve.menu,  "Swap  Buffers  Xf",  swap); 
addtopup(curve_menu,  "Quit  Xf",  quit); 

> 

. . .  FUNCTION:  REMAKE.DISPLAY.MENUO  •*•***•***************/ 

void  remake_display_menu(void) 

char  menu_string[32] ; 

freepup(di splay .menu) ; 
display.menu  =  newpupO; 

addtopup(display_menu,  "Display  Mode  Xt  XF",  change.display) ; 
addtopup (display.menu,  "Single  Bar  k  Continuous  XxO"); 
addtopup(di8play_menu ,  "Single  Bar  k  Interval  (axis  xing)  Xxl"); 
addtopup (display.menu ,  "Single  Bar  k  Interval  (peak)  Xx2") ; 
addtopup (display.menu ,  "Time-Shift  k  Continuous  Xx3"); 
addtopup (display.menu,  "Time-Shift  k  Interval  (axis  xing)  Xx4”); 
addtopup(display_menu ,  "Time-Shift  k  Interval  (peak)  XxS"); 

> 
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/*•••***••*•*•«.•••••*•  FUNCTION:  REMAKE. INTERVAL.MENUO  •*,•**•••»•*.*.*••**•/ 

void  remake.interval.menu(void) 

{ 

char  menu_string[32] ; 

freepup(interval.menu) ; 
interval.menu  =  newpupO; 


> 


addtopup ( interval.menu , 

addtopup ( interval.menu , 
addtopup (interval .menu , 
addtopup(interval_menu , 
addtopup (interval .menu , 
addtopup ( interval.menu , 
addtopup ( interval.menu , 
addtopup ( interval.menu , 
addtopup (interval .menu , 
addtopup(interval_menu , 
addtopup (interval.menu , 


"Set  Display  Interval  */.t  V.F", 

set.interval.skip) ; 
"  Every  2nd  event  Xx2"); 

"  Every  3rd  event  Xx3"); 

"  Every  4th  event  *4x4"); 

"  Every  5th  event  *4x5"); 

"  Every  6th  event  Xx6"); 

"  Every  7th  event  *4x7"); 

"  Every  8th  event  Xx8"); 

"  Every  9th  event  *4x9"); 

"  Every  10th  event  XxlO"); 

"  All  events  (Reset)  Xxl"); 


/* 


FUNCTION :  REMAKE.BIN.MENUO 


void  remake.bin.menu(void) 

{ 

char  menu_string[32] ; 

freepup(bin.menu) ; 
bin.menu  s  nevpupO  ; 


addtopup (bin. menu, 

addtopup (bin.menu , 
addtopup ( bin.menu , 
addtopup(bin.menu , 
aadtopup(bin.menu , 
addtopup (bin.menu , 


Set  »  of  Time-Shift  Bins  *4t  *4F" , 

set.bins) 

6  Xx6") ; 

10  XxlO")  ; 

24  Xx24") ; 

50  y,x50")  ; 

74  *4x74"); 


FUNCTION:  REMAKE.START.MENU ( ) 


void  remake_start_aenu(void) 

{ 

char  menu.string [32]  ; 


/ 


/ 
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freepup(start.menu) ; 
start. menu  =  nevpupO; 

addtopup (start .menu,  "Start  Display  at...  Xt“); 
addtopup (atart.men  "Beginning  of  Data  Xf",  2ero_ptrs); 
addtopup(staLrt.menu,  "Current  Position  Xf",  leave.ptrs) ; 
addtopup (start.menu,  "Selected  Time  Xm",  timer.menu); 

> 


/*m***«*. FUNCTION :  REHAKE.TIMER.MENUO  ***•*•*..*.*****..*.*./ 

void  remake_timer_aenu( void) 

{ 

char  menu.string [32] ; 

freepup(timer.menu) ; 
timer.menu  =  nevpupO; 

addtopup (timer .menu ,  "Elapsed  Time  (in  mins: secs)  Xt  XF", 

position.ptrs) ; 

addtopup(timer_menu,  "  1  Xx60"); 

addtopup (timer.menu,  "  2  Xxl20"); 

addtopup(timer_menu,  "  3  XxlSO"); 

addtopup (timer.menu,  "  4  Xx240"); 

addtopup (timer.menu,  "  4:4SXx285"); 

addtopup (timer.menu,  "  5  Xx300"); 

addtopup(timer_menu,  "  6  Xx360"); 

addtopup (timer.menu,  "  8  Xx480"); 

addtopup(timer_menu,  "  10  Xx600"); 

addtopup (timer.menu,  "  12  Xx720"); 

} 

/*»•*•«•**•****•••*•»•••  FUNCTION:  REHAKE.SYNCH.MENUO  »••**.*.*******.«*****/ 

void  remake_8ynch_menu(void) 

{ 

freopup(3ynch.menu) ; 
aynch.menu  =  netfpupO; 

addtopup (synch.menu,  "Reference  Channel  Xt  XF" ,  set.b) ; 

addtopup (aynch.menu,  "  CH  1  Xxl"): 

addtopup (synch.menu,  "  CH  2  Xx2"); 

addtopup (aynch.menu,  "  CH  3  Xx3"); 

addtopup (aynch.menu,  "  CH  4  Xx4“); 

addtopup(aynch_menu,  "  CH  5  Xx5"); 

addtopup(aynch_menu,  "  CH  6  Xx6”); 

addtopup (aynch.menu,  "  CH  7  Xx7"); 

addtopup  (synch.menu,  "  CH  8  Xx8''); 

addtopup (aynch.menu,  "  CH  9  Xx9”); 

addtopup (synch.menu,  "  CH  10  XxlO") ; 
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> 


addtopup(synch_menu , 
addtopup(synch_menu, 
addtopup (synch. menu , 
addtopup(synch_menu , 
addtopup ( synch.menu , 
addtopup(aynch_menu, 
addtopup(synch_menu , 
addtopu” ( synch.menu , 


'  CH  11  Xxll") ; 

CH  12  2j12"); 

CH  13  y.xl3"); 

CH  14  Xr.14"); 

0.5  Hz  Sine  ’/.xlS") 

1  Hz  Sine  */,xl5") 
1.5  Hz  Sinn  y.xl7") 

2  Hz  Sine  '/,xl6") 


/ 


FUNCTION:  REMAKE. LINE. MENU ( ) 


/ 


void  remake.line.menu(void) 

{ 

freepup (line.menu) ; 
line.menu  =  newpupO; 

addtopup(line_menu ,  1 
addtopupdine.menu,  ' 
addtopupdine.menu ,  ' 
addtopupdine.menu,  ’ 
addtopupdine.menu,  ' 
addtopupdine.menu,  ' 
addtopupdine.menu,  ' 
addtopupdine.menu,  ' 

> 

/»**! 


Pixel  Width  */.t  'IS",  set.lw) ; 
2  '/.x2") ; 

5  '/.X5" ) ; 

10  y.xlO"); 

15  '/,xl5") ; 

20  '/,x20") ; 

25  y,x25") ; 

30  y,x30") ; 


REMAKE. INTENS ITY.MENU ( ) 


FUNCTION: 


void  remake.intensity.meniHvoid) 

{ 

freepup(intensity_menu) ; 
intensity.menu  =  newpupO; 


addtopup(intensity_menu , 
addtopupdntensit.y.menu , 
addtopup(int.ensity_menu , 
addtopup(intensity_menu , 
addtopup (intensity.menu, 
addtopup (intensity.menu , 


"Bar  Intensity  ’/.t  '/.F" ,  set.intenBity)  ; 
"Scale  A  Clip  3  3  Sigmas  '/,x0"); 

"Clamp  to  Maximum  ’/,xl"); 

"Threshold  ®  Reference  <  0.005  */,x 2"); 
"Threshold  ®  Reference  <  0.01  ’/,x3"); 
"Threshold  ®  Reference  <  0.02  Xx4"); 


/ 


> 

. . .  FUNCTION:  REMAKE.TEST.KENUO  **** 

void  remake.test.menu(void) 

{ 

f reepup(test.menu) ; 
test.menu  ~  newpupO; 

addtopup  (test. menu ,  "Test  Modes  '/.t  %F"  ,  set.test) ; 


152 


addtopup(teat_menu ,  "Step  Draw  ON  */,xl"); 
addtopup(teat_menu ,  "Step  Draw  OFF  */,x2") ; 
addtopup(test_menu ,  "Batch  Draw  ON  */,x5") ; 
addtopup(teat_menu ,  "Batch  Draw  OFF  %x6"); 
addtopup(t98t_menu ,  "Single  Buffer  ON  7iX3"); 
addtopup(tsat_menu ,  "Single  Buffar  OFF  %x4"); 
addtopup(test_menu ,  "TEST  MODE  OFF  %x0") ; 

> 

/•*•*****«»»*»•*«•»•••••  FUNCTION:  REMAKE.ALL.MENUSO  »*•*«•»****•*•»••****•/ 

void  ramake_all_menu8 ( ) 

{ 

remake_display_menu() ; 
remake_intei  v?.l_menu() ; 
remake_bin_mer.u() ; 
remake_timer.aenu(. )  ; 
reoake_start_menu( ) ; 
remako_synch_menu() ; 
remike_line.menu() ; 
reoake_iutenaity_menu( ) ; 
remake_to8t.aenu() ; 
renake_curve_monu() ; 

> 
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1-4  Display  Option  Control:  tcontrol.c 


/* 

» 

•  tcontrol.c 

* 

*  Part  of  "TOPOS"  by  Capt  Dwight  Roblyer,  AFIT/GS0-92D 

*  November  1992 

* 

«  Based  on  "CURVE  DEMO" 

»  by  Howard  Look  for  Silicon  Graphics,  June  1989 

* 

*  This  file  contains  the  routines  used  for  the  program’s  display  and 

*  environment  options,  as  defaults  and  as  selected  by  the  user.  Most 

»  of  these  routines  are  called  by  menu  selection.  See  the  file  tmenu.c 

*  for  menu  definitions. 

* 

*/ 


•include  <stdio.h> 
•include  <gl.h> 
•include  <device.h> 
•include  <math.h> 
•include  "topos.h" 
•include  "tevent.h" 


/*  Prototypes  •/ 

void  set_initial_conditions(void) ; 

void  8et_defaults(void) ; 

void  change_display(int) ; 

void  change_8cale(int) ; 

void  aet_interval_skip(int) ; 

void  set.bins(int) ; 

void  set_b(int); 

void  set.lw(int) ; 

void  set. intensity (int) ; 

void  leave.ptrs(void) ; 

void  zero.ptrs(void) ; 

void  position. ptrs(int) ; 

void  awap(void) ; 

void  dbl.synch(void) ; 

void  set.testdnt) ; 

extern  void  clear_vindow( void) ; 
extern  void  remake_curve_menu(void) ; 
extern  void  remake_display_menu(void) ; 
extern  void  remake.all.menus(void) ; 


int  8kip.ainutes=0 , 

8kip_seconds=0 , 
total_skip=0 , 
interval.skip; 
line.atyle ; 

float  synch.limit; 

Boolean  time.shift, 
interval , 
axis.crossing, 
peak , 

constant.intensity , 
zero.pointers , 
wait , 

single.buff , 
synch. constrained , 
batch ; 

extern  Boolean  drau.active,  noprep; 
extern  char  tool; 
extern  float  myBcale; 
extern  int  hardware, 

saaple.rate , 

8  , 

bins , 

b, 

lw, 

ch; 

extern  FILE  *ofp[] ; 

/*  About  the  window...  •  / 

extern  int  origin.x,  origin.y,  size.x,  size.y; 
extern  float  aspect; 

/*•**********«*****• *****FUNCTIOH:  SET.INITIAl.CONDITIONSO***** 

/* 

*  Sets  default  values  for  the  TOPOS  options.  THOSE  WHICH  USERS 

*  MIGHT  WANT  TO  CHANGE  ARE  MARKED. 

•/ 

void  set. initial. conditions (void) 

{ 

ortho(-l .0,  1.0,  -1.0,  1.0,  -2.0,  2.0); 
line. style  =  SOLID; 
setlineotyledine.style)  ; 
unqdevice(INPUTCHANGE) ; 


noprep  =  TRUE; 


/*  which  draw  events  to  display  */ 

/«  pixel  height  of  bar  */ 

/*  time. shifted  or  single  bar  */ 

/•  •  of  bins  in  time-shifted  display*/ 

/*  continuous  or  interval  drawing  */ 

/*  if  interval,  then  axis-xing  */ 

/*  if  interval,  then  peaks  */ 

/*  starts  at  beginning  of  data  file  */ 

/*  reference  channel  */ 

/*  chl5=lhz  chl6=2hz  chl7=1.5hz  chl8=.5hz  */ 

/*  single  or  double  screen  buffers  */ 

/*  step-draw  if  TRUE  */ 

/»  batch-draw  if  TRUE  */ 

constant.intensity  *  TRUE;  /*  bar’s  color  is  constant  or  scaled  */ 

synch_constrained=FALSF,;  /•  bar’s  color  is  thresholded  */ 

synch_limit=0.01 ;  /*  threshold  values  for  colors  */ 

> 

/********•*****.*•••*.•*• FUNCTION:  SET.DEF AULTS ( . . . 

void  set_defaults(void) 

< 

set_initial_conditions() ; 
remake.all.menusO ; 
draw.displayC) ; 

> 

/********************* ***FUNCTI0N :  CHANGE_DISPLAY( )**************************/ 

void  change_display(int  new. mode)  /*  implements  Display  Mode  menu  option  */ 

{ 

if  (new_mode==0) 

{ 

time.shiit  =  FALSE; 
interval  =  FALSE; 
axis_cros8ing  =  FALSE; 
peak  =  FALSE; 

> 

if  (new_mode==l) 

{ 

tioe.shift  =  FALSE; 
interval  =  TRUE; 


interval_skip=l ; 
lw=25; 

time.shift  =  TRUE; 
bins=6; 

interval  =  TRUE; 
axis.crossing  =  TRUE; 
peak  «  FALSE; 

skip.minutes  =  0; 

b=17 ; 

s ingle. buf f=TRUE; 

wait=FALSE; 

batch=TRUE; 
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axia.crossing  =  TRUE; 
peak  =  FALSE; 

> 

if  (new_mode==2) 

{ 

time. shift  =  FALSE; 
interval  =  TRUE; 
axis.crosaing  =  FALSE; 
peak  =  TRUE; 

> 

if  (new_mode==3) 

< 

time. shift  -  TRUE; 
interval  =  FALSE; 
axis. crossing  =  FALSE; 
peak  =  FALSE; 

> 

if  (new_mode==4) 

{ 

time.shift  =  TRUE; 
interval  ■  TRUE; 
axis.crosaing  =  TRUE; 
peak  =  FALSE; 

> 

if  (nev_mode==5) 

{ 

time.shift  =  TRUE; 
interval  =  TRUE; 
axis.crossing  =  FALSE; 
peak  =  TRUE; 

> 


reaake.display.menuO  ; 
remake. curve.menuO ; 


/****♦ . * . ‘•♦••♦••FUNCTION :  SET. INTERVAL.SKIPO ••••••••••••••••• 

void  set_interval_8kip(int  skip)  /*  implements  Draw  Interval  menu  opt 

{ 


if  (skip) 

intervul.skip  =  skip; 

else 

interval.Bkip  =  1; 


•••••/ 

ion  •/ 


157 


THIS 

PAGE 

IS 

MISSING 

IN 

ORIGINAL 

DOCUMENT 


if  (intensity==4) 

{ 

constant_intensity=FALSE; 
synch_constrained=TRUE ; 
synch.limit=0 . 02 ; 

> 

> 


/••****•*«•***•**••*•• •••FUNCTION:  ZERO.PTRS ()***• . **•****•***,.*/ 

void  zero.ptrs(void)  /*  implements  Start  Display  (Beginning)  menu  option  */ 

{ 

s=0; 

zero.pointere=TRUE ; 
skip.seconds  =  0; 
draw.displayO  ; 

> 


/*••*.**•••*.*••*****• ••♦FUNCTION:  LEAVE.PTRS ()•***•*»***•***** ************/ 

void  leave.ptrs ( void)  /*  implements  Start  Display  (Current)  menu  option  */ 

{ 

zero_pointers=FALSE; 
skip.seconds  =  0; 
draw.displayO  ; 

> 


/***,*******.***********FUNCTI0N:  POSITION.PTRS ()•*•*♦**••**•*************•/ 

void  position.ptrsdnt  skip)  /*  implements  Start  Display  ($)  menu  option  */ 

{ 

zero_pointers=TRUE; 

s=0; 

skip.seconds  =  skip; 
draw.displayO ; 

> 
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/ 


FUNCTION:  SWAPO 


/ 


void  svap(void) 

{ 


> 


while( !qtest() ) 

svapbuff ars( ) 


/* 


implements  Swap  Buffers  menu  option  */ 


******  *  ****  ***************  FUNCTI ON :  SET.TEST ()****•*******♦**************•/ 

void  8«t_test(int  mode)  /*  implements  Test  Modes  menu  option  */ 

{ 

if  (mode==0) 

{ 

wait=FALSE; 
single.buff “FALSE; 
batch=FALSE; 

> 

if  (mode*“l) 

{ 

vait=TRUE; 

batch=FALSE; 

> 

if  (mode==2) 

{ 

vait=FALSE; 

bat.ch=FALSE; 

> 

if  (node*“3) 

single  'o'jff=TRUE; 

if  (mode==4) 

si!>£'.le_buff  “FALSE; 

if  (mode®,=5) 

{ 

batch=TRUE; 
wait “FALSE; 

> 

if  (aode==6) 

{ 

batch=FALSE ; 
uait=FALSE; 

> 

} 


ICO 


1.5  Event  Queue  Control:  tevent.c 


NOTE:  The  following  code  was  modified  very  little  during  this  research. 


/* 

«  tevent.c 

* 

*  Part  of  "TOPOS"  by  Capt  Dwight  Roblyer,  AFIT/GS0-92D 

*  November  1992 

* 

*  Used  with  few  changes  from  "CURVE  DENQ" 

*  by  Howard  Look  for  Silicon  Graphics,  June  1989 

* 

•  A  more  rational  way  to  handle  reading  the  event  queue 

•  Written  by  Wade  Olsen 

• 

•/ 

•include  <stdio.h> 

•include  <device.h> 

•include  "tevent.h" 


typedef  struct  eveDt.s 

{ 

int  window,  device,  state; 

void  (*func) (void  *,  int); 
char  *arg; 

struct  event.s  *next; 

>  event.t,  *event_p; 

typedef  struct  update.s 

{ 

int  *flag; 

void  (*ufunc) (void  *) ; 
char  *arg; 

struct  update. s  *next; 

>  update.t,  *update_p; 


static  event.p  ovent.list; 

static  update.p  update.list; 

/* 

*  This  routine  adds  an  event  to  be  checked  for  to  the  event  queue. 
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*  window  should  be  the  vid  of  the  window  to  respond  in,  or  ANY  if 

*  this  event  applies  to  all  windows,  device  is  the  device,  and  state 

*  is  the  device’s  value  (e.g.  ANY,  UP,  DOWN,  the  window  id  (for 

*  REDRAW),  etc).  Func  is  the  function  that  will  be  called,  with 

*  arguaents  ’arg’  and  the  device’s  value, 
e 

*  NOTE:  the  device  Bust  be  queued  for  it  to  be  found  by  the  event () 

*  routine —  aad.event  DOES  NOT  qdevice(device) . 

*/ 

void 

add_event (window ,  device,  state,  func,  arg) 
int  window,  device,  state; 
void  (*func)(void  *,  int); 
char  *irg; 

{ 

event.p  new_guy; 


new.guy 
new.guy ->window 
nev_guy->device 
new_guy->state 
new_guy->func 
new_guy->arg 
uew_guy->next 


( event _p) malloc (sizeof (event.t) ) ; 

window ; 

device ; 

state ; 

func ; 

arg; 

event. list ; 


event  ..list  «  new_guy; 


/• 

*  Specify  a  function  to  be  called  if  there  is  nothing  in  the  queue 

*  and  (*flag)  is  non-zero.  If  no  update  function  is  active,  or 

*  (*flag)  is  0,  then  eventO  will  block  on  a  qread.  If  there  is  an 

*  active  update  function,  eventO  will  continuously  call  the  update 

*  function,  hogging  the  cpu. 

*/ 

void 

add_update(flag ,  ufunc,  arg) 
int  *flag; 

void  (*ufunc) (void  *) ; 
char  *arg; 

{ 

update.p  new.guy; 


new.guy  *  (update. p)nalloc (aizeof (update.t) ) ; 
new_guy->f lag  *  flag; 
new_guy->ufunc  =  ufunc; 

new.guy->arg  *  arg; 
new.guy->next  *  update.list; 
update.list  =  new.guy; 


/* 
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*  The  aain  Event.  Call  this  repeatedly  to  automagically  handle 

*  reading  the  queue,  and  calling  your  functions  to  handle  what 

*  appears  there . 

•/ 

void 

eventO 

{ 

void  f ind.event(void) ,  event. inputchange(void) ; 
int  f ind.update(void) ; 
static  int  initialized  =  0; 

if  (initialized  ==  0) 

{ 

add_event(ANY ,  INPUTCHANGE,  ANY,  event.inputchange ,  NULL); 
initialized  =  1; 

> 


> 


/•  Something  in  the  queue?  Handle  it  */ 
if  (qtestO) 

f ind_event( ) ; 


int 

f  ind.updateO 

{ 

update. p  scan; 

int  updated  =  0; 

for  (scan  *  update.list;  scan  kk  updated  0;  scan  =  scan->next) 

if  (*8can->f lag) 

{ 

(*scan->ufunc) (acan->arg) : 
updated  *  1 ; 

} 

> 

return(updated) ; 

> 

int  context,  state,  device; 
void 

event.inputchange () 

{ 

context  =  vingetO; 

> 


void 

find.  eventO 

{ 
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event.p  scan; 
short  s ; 

device  =  qread(fcs); 
state  =  s; 

for  (scan  7  e'’ent_liat,  scan;  scan  =  scan->next) 

{ 

if  (  ((scan->window  =*  ANY)  i|  (context  *»  scan->windov)) 

tt  ((scan->device  =■  ANY)  1 1 (device  ■»  scan->device) ) 
tt  ((scan->state  ==  ANY)  II  (state  ==  scan->state) )) 

{ 

(*scan->func) (scan->arg,  state); 
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1.6  Header  Flits 


/* 

* 

*  topos.h 

«  This  fils  contains  the  defines  used  in  TOPOS. 

* 

/*  Linestyles  •/ 

•define  SOLID  0 

•define  DOTTED  1 

•define  DASHED  2 

/*  Patterns  for  the  linestyles...  */ 

/*  2  pixels  on,  2  pixels  off  •! 

•define  DOTTED. PATTERN  OxCCCC 
/*  8  pixels  on,  8  pixels  off  •/ 

•defit  DASHED.PATTERN  OxFFOO 


•define  ON  1 

•define  OFF  0 

/*  Hardware  types  */ 
•define  ECLIPSES  1 

•define  ECLIPSE24  2 

•define  GT  3 


/* 

*  event. h 

*  External  interface  and  defines  to  input-queue  event  handling 

*  routines. 

*  Written  by  Wade  Olsen  for  Silicon  Graphics,  Inc. 

*/ 

/* 

*  The  event  handler  understands  two  kinds  of  things;  events  and 

*  updates.  Events  are  reactions  to  things  occuring  in  the  input 

*  queue.  Updates  are  functions  that  should  be  called  whenever  there 

*  is  nothing  waiting  in  the  input  queue,  and  say  be  active  or 

*  inactive.  If  there  are  no  active  updates  and  nothing  in  the  input 

*  queue,  then  eventO  will  block,  using  up  no  CPU  time. 

* 

*  add. event  is  used  to  look  for  events.  The  first  three  arguments 

*  are  used  to  identify  which  event  to  look  for.  The  first  argument 

*  is  the  window  (gid)  the  event  must  happen  in;  if  this  value  is  ANY 

*  then  any  window  will  do.  Th6  second  argument  is  the  device  to  look 
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*  for  (e.g.  RIGHTMOUSE  or  REDRAW  or  KEYED,  etc).  Again,  if  it  is 

*  ANY  then  any  device  will  natch.  The  third  argument  is  the  value 

*  the  device  must  generate  (e.g.  DOWN  or  UP);  ANY  means  all  values 

*  match.  The  last  two  arguments  are  what  should  be  done  when  an 

*  event  is  generated.  The  fourth  argument  is  a  function  to  be 

*  called,  and  the  fifth  is  an  argument  that  should  be  supplied  to  the 

*  function.  In  addition,  the  value  generated  by  the  device  will  also 

*  be  passed  to  +he  function  when  it  is  called. 

* 

*  For  example, 

*  add_event(winget() ,  RIGHTMOUSE,  DOWN,  dopup,  my. menus) ; 

*  qdevice (RIGHTMOUSE) ; 

*  will  make  a  pop-up  menu  appear  when  the  right  mousebutton  goes 

*  down.  Note  that  you  must  do  the  qdeviceO  call  yourself. 

*/ 

void  add.evont(int ,  int,  int,  void  (*fn)(void  *,  inti,  char  *)  ; 


/• 

*  An  update  is  like  an  event,  only  simpler.  The  first  argument  is  a 

*  pointer  to  an  integer  flag  specifying  whether  or  not  thiB  update 

*  function  is  active.  The  second  is  a  function  to  be  called  when  it 

*  is  active,  and  the  last  is  an  argument  to  be  supplied  to  the 

*  function. 

*/ 

void  add_update(int  *,  void  (*fn)(void  *),  char  *)  ; 

/* 

*  Finally,  when  all  updates  and  events  have  been  added,  repeatedly 

*  call  event ()  to  handle  them  --  something  like 

* 

*  while  (quitflag  ==  FALSE)  eventO; 

* 

«  You  should  have  previously  added  an  event  that  sets  quitflag  to 

*  TRUE,  of  course. 

•/ 

void  event(void)  ; 

/* 

*  These  are  some  useful  defines  for  the  possible  values  buttons 

•  can  generate . 

•/ 

•define  ANY  -1 

•define  UP  0 

•define  DOWN  1 

/• 

•  And  a  few  external  variables  you  might  find  useful 

*/ 

extern  int  context,  state,  device; 
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1.7  Makefile 


*  Ttaifl  makefile  creates  the  test  AFIT  toposcopeo,  TOPOS 

* 

LLDLIBS  *  -Ifm.s  -lgl.s  -lm 

CFILES  =  topoe.o  tevent.o  tmenu.o  tdraw.o  tcontrol.o 
topos:  ${CFILES> 

<TAE  *  cc  -g  -o  topos  ${CFILES>  ${LLDLIBS> 
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